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Abstract 

by 


KELLY SCOTT CARNEY 


A method of dynamic substructuring is presented which utilizes a set of static 
Ritz vectors as a replacement for normal eigenvectors in component mode synthesis. 
This set of Ritz vectors is generated in a recurrence relationship, proposed by Wilson, 
which has the form of a block-Krylov subspace. The initial seed to the recurrence 
algorithm is based upon the boundary flexibility vectors of the component. 
Improvements have been made in the formulation of the inial seed to the Krylov 
sequence, through the use of block-filtering. A method to shift the Krylov sequence 
to create Ritz vectors that will represent the dynamic behavior of the component at 
target frequencies, the target frequency being determined by the applied forcing 
functions, has been developed. A method to terminate the Krylov sequence has also 
been developed. Various orthonormalization schemes have been developed and 
evaluated, including the Cholesky/QR method. Several auxiliary theorems and proofs 
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which illustrate issues in component mode synthesis and loss of orthogonality in the 
Krylov sequence have also been presented. 

The resulting methodology is applicable to both fixed and free-interface 
boundary components, and results in a general component model appropriate for any 
type of dynamic analysis. The accuracy is found to be comparable to that of 
component synthesis based upon normal modes, using fewer generalized coordinates. 
In addition, the block-Krylov recurrence algorithm is a series of static solutions and 
so requires significantly less computation than solving the normal eigenspace problem. 
The requirement for less vectors to form the component, coupled with the lower 
computational expense of calculating these Ritz vectors, combine to create a method 
more efficient than traditional component mode synthesis. 
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Chapter 1 


Introduction 


Component mode synthesis is a methodology for analyzing large structures by 
separating them into smaller components, reduced representations of which can then 
be recombined to analyze the entire system. This methodology has become well 
established and widely used in structural dynamic analysis. The advantages of 
component mode synthesis include lower computation costs associated with smaller 
components, and the flexibility of data management gained by working with the 
discrete components. 

The typical component mode synthesis algorithm is briefly described 8 . A large 
structure is broken into components, with each component having a set of boundary, 
or interface, points. At these interface points, fixed or free boundary conditions are 
assumed, and a corresponding set of component normal mode shapes, or eigenvectors, 
is determined. The eigenvectors are augmented by a set of modes which are 
associated with the component’s boundary flexibility. Depending on whether a fixed 
or free interface is selected, these modes are the constraint modes or the attachment 
modes, respectively. The combined set of component normal modes and boundary 
modes are used to represent the component in subsequent system analysis by using 
the following transformation process. The combined set of modes form a coordinate 
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transformation matrix which transforms the physical coordinates of the structural 
model into a combination of modal coordinates and boundary coordinates. The 
boundary coordinates are retained in the physical space, so they can be used to couple 
the components for subsequent system analysis. 

A component’s size, although smaller than that of the complete structural 
model, can still be large enough to be computationally expensive. The rapid reduction 
in cost per calculation in today’s digital computers has not necessarily led to a 
reduction in total computation cost. Instead, engineers have exploited the increased 
computational resources by creating larger structural and component models. The 
larger models have allowed for more structural details to be represented, as well as 
more refined data recovery, but they may be expensive to formulate and analyze. In 
order to reduce the computational cost associated with large component models, it is 
desirable to develop more efficient methods of formulation. Since the solution of the 
normal eigensystem problem requires the largest computational effort in component 
formulation, it is logical to develop alternate methods which circumvent the 
eigensystem solution entirely. 

A method, which does circumvent the eigensystem solution, has been defined 
in literature and is briefly described 110 ' 30 . The boundary flexibility modes, specifically 
either the same constraint modes or attachment modes that were mentioned previously, 
are multiplied by the component mass matrix to create a force matrix. Static analysis 
is then performed, using this force matrix and the component stiffness matrix, to 
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obtain a matrix, or block, of vector displacements. A recurrence relationship of matrix 
multiplications, which have been shown to be a Krylov sequence 10 19 , then defines a 
series of matrices, or blocks, of vector displacements. The calculated vectors are 
orthogonalized, using normalized Gram-Schmidt orthogonalization 18 . These vectors, 
which can also be thought of as static modes or static Ritz vectors, replace the normal 
modes in the component formulation methodology. Because the static Ritz vectors 
are calculated in blocks and are based on a Krylov sequence, the subspace defined by 
these vectors is called a block-Krylov subspace. 

The work described in this dissertation originated by identifying areas of 
potential improvement in the implementation of the existing static Ritz vector, block- 
Krylov, boundary flexibility methodology. Specific potential improvements in the 
form of equations and software were created, implemented and assessed. If useful, 
the improvement was adopted. These adopted improvements are briefly discussed 
in the following paragraphs. 

Mathematical theory predicts that the vectors obtained from the Krylov 
sequence, after orthogonalization with the previous two vectors in the sequence, are 
independent. However, in practical applications these vectors usually converge, and 
sometimes quickly, to a nearly dependent state. Using totally dependent or nearly 
dependent vectors, the Gram-Schmidt orthogonalization algorithm fails. As a result, 
a set of vectors sufficient to define the dynamics of the component may not be 
obtainable. A solution to this problem has been available, consisting of repeated 
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Gram-Schmidt orthogonalization with all previous vectors, whenever that procedure 
initially fails. This simple solution is usually successful. However, successive re- 
orthogonalization can be computationally expensive and there is no guarantee that the 
resulting vectors will have any physical significance, or that this process will not also 
eventually fail 22 . Several alternate re-orthogonalization procedures have been 
investigated, evaluated, and discussed in this work. 

In addition to the investigation of Gram-Schmidt re-orthogonalization 
procedures, an alternative procedure for orthogonalization in the block-Krylov 
sequence has been presented. Rather than the Gram-Schmidt vector by vector 
numerical scaling approach, a matrix transformation which orthogonalizes the vectors 
simultaneously can be created, without calculating eigenvectors. This matrix 
transformation utilizes the transpose of the Cholesky factor of a matrix product of the 
original dependent set of vectors. This process has been shown to be more efficient 
than the Gram-Schmidt orthogonalization procedure in this application. Cholesky 
orthogonalization has been integrated into the boundary flexibility method of 
component synthesis using generalized static Ritz vectors. 

Another enhancement in the boundary flexibility methodology was possible 
because of the initiation of the Krylov sequence with the boundary flexibility modes. 
Recall that the boundary flexibility modes are multiplied by the mass matrix to create 
a force matrix. The size of this force matrix is the same as the size of the component 
interface, which subsequently determines the block size in the Krylov sequence. 
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When working with finite element models which have simple interfaces the resulting 
force matrix and solutions have physical significance and the blocks are convenient 
to process. However, in typical aerospace structures where a more complicated 
interface and a large number of interface nodes exists, the Krylov block size will be 
very large. This can lead to several complications in convergence and truncation. 
The problem of large block size has been solved in this work by discarding nearly 
dependent vectors from the block, previous to orthogonalization. Vectors which are 
nearly dependent are, geometrically, nearly identical to each other. As a result, no 
particularly useful information is being lost when they are removed from the block. 
This reduction also reduces the size of subsequent blocks. This method has been 
identified in this work as block filtering. 

A problem which all static Ritz vector, Krylov sequence solutions have is the 
lack of a sound mathematical basis for judging when to terminate the sequence. The 
sequence should be terminated when the resulting set of Ritz vectors is sufficient for 
dynamic representation of the component. When using normal eigenvectors to 
represent a component, modal truncation, based upon an eigenvalue cutoff, is the 
most popular basis forjudging if the dynamic representation is sufficient. Static Ritz 
vectors do not have an eigenvalue with which to associate a truncation limit and so 
an alternate method of the termination of the Krylov sequence was required and 
developed. This method is an heuristic methodology based upon the density of the 
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modal space in the component and the demonstrated observation that the most 
dynamically significant static Ritz are generated early in the Krylov sequence. 

While investigating the Lanczos eigenvalue extraction method, done so because 
of it’s use of the Krylov sequence and it’s resulting similarity with static Ritz vectors, 
the tool of inverse operators was recognized. The use of inverse operators is 
commonly called sequence shifting. By shifting the sequence, an eigenvalue 
extraction routine can locate missing eigenvalues within a specific frequency range. 
The use made of shifting in this work differs from that of eigenvalue extraction in 
that, instead of targeting ranges of missing eigenvalues, the shift is targeted to the 
frequency of the applied forcing function. In this manner, vectors which can 
contribute to an accurate dynamic response prediction are generated. The use of 
shifting in this manner has been identified in this work as targeted shifting. 

Of the various enhancements discussed in the preceding paragraphs, only block 
filtering is specifically tied to component mode synthesis and the boundary flexibility 
method. All the other improvements, such as targeted shifting and modal density 
series truncation, can be used in the general static Ritz vector methodology. Even 
so, in this work these improvements have primarily been implemented and evaluated 
in the context of the boundary flexibility method. Several general math proofs, which 
are also not specific to the boundary flexibility method, have also been derived and 
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Chapter 2 


Literature Review 

2.1) Introduction: 

Wilson, Yuan, and Dickens 29 originally proposed the use of Ritz vectors, based 
upon external loading, for structural dynamic analysis. This formulation reduced an 
entire structure, not a component. The algorithm begins with a set of externally 
applied loads. The displacements from the static solution to the applied loads become 
the initial Ritz vector. That vector is then multiplied by the mass matrix to become 
the next force vector. This sequence is repeated to form a recurrence relationship and 
a series of special Ritz vectors, which are referred to in this work as static Ritz 
vectors. The proposed recurrence relationship is used in the papers discussed below 
and throughout this work. 

Wilson’s methodology was applied, using MSC/NASTRAN, to several large 
finite element models by Arnold, Citerley, Chargin, and Galant 3 . It was found that 
Wilson’s methodology was computationally more efficient than the standard normal 
modes procedure. A recent application using a simple model of the space station was 
presented by Escobedo-Torres and Rides 12 . This work compared the predicted 
transient response of the space station due to a docking force using load dependent 
Ritz vectors with predictions using eigenvectors. 
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Nour-Omid and Clough 19 investigated Wilson, et al’s methodology and found 
that the proposed recurrence relationship actually generated a Krylov sequence. A 
Krylov subspace of order j is a vector space defined by 

[<|>,A<| (21) 

where <j) is a column vector and A is a square, symmetric matrix. If A is n x n, 
nonsingular, and if j = n, the Krylov vectors span the n dimensional space 10 , and an 
exact solution can be produced. In structural dynamics, the Krylov subspace can be 
defined by the following Krylov sequence, 

[r,K~ l Mr t (K~ x r] ( 2 - 2 ) 

where K is the stiffness matrix, M is the mass matrix, and r is a starting vector (or in 
block-Krylov, a set of vectors). The matrix product K 'M is not symmetric and its use 
in the Krylov sequence yields subtle theoretical and practical differences which will 
be discussed below. The Krylov sequence is also the basis of the Lanczos eigenvalue 
extraction algorithm, and Nour-Omid and Clough refer to Wilson’s static Ritz vectors 
as Lanczos coordinates. The Lanczos eigenvalue extraction algorithm generates 
Krylov vectors, which are used to transform the system into a tridiagonal form 1522 . 
To extract the eigenvalues, this tridiagonal matrix is diagonalized using a QR, or 
related, algorithm. Nour-Omid and Clough utilize the tridiagonal system matrices 
directly to solve the dynamic response problem. 



9 


Nour-Omid and Clough extended their work to a more general dynamic 
loading represented by a linear combination of load vectors 20 , such as time dependent 
loading. The result was a structural dynamics application of the block-Krylov, or 
block-Lanczos method for the dynamic analysis of structures. A block is defined as 
the combination of a set of vectors, or modes, into matrices. They also presented the 
requirement of using the Gram-Schmidt procedure to orthogonalize the vectors within 
the Krylov block. 

The use of static Ritz vectors was shown to be applicable to component mode 
synthesis by Wilson and Bayo 28 . The static Ritz vectors calculated were based, once 
again, upon an external load. Only a formulation for components with fixed interface 
boundary conditions was presented. This work was also implemented and applied by 
Leger 16 to an example small beam. 

A similar development of static Ritz vectors in component mode synthesis was 
presented by Allen 2 . This paper provided the basis for the application work performed 
by Brunty 5 . The transient response of the Space Shuttle vehicle, during liftoff, was 
calculated using load-dependent static Ritz vectors and compared to the response 
predicted using eigenvectors in the component mode synthesis. Similar answers were 
obtained using less computer time. 

Abdallah and Huckelbridge 1 , and independently, Craig and Hale 10 , 
demonstrated a generalized methodology applicable to components with fixed or free 
interfaces, with or without rigid body modes, and with or without applied loading. 
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Components having no applied external loading were formed using the boundary 
flexibility matrix, multiplied by the mass matrix, to form a force matrix. This force 
matrix produces a set of static Ritz vectors. (Craig and Hale refer to these vectors as 
Krylov vectors). The boundary flexibility matrix is defined as either the constraint 
modes or the attachment modes, depending on whether fixed or free interface 
conditions are selected. The methodology contained in these two papers is reviewed 
in the next three parts of this section. Abdallah and Huckelbridge also quantified the 
advantages, in computational effort, that generalized static Ritz vectors have over 
normal eigenvectors. Carney, Abdallah, and Huckelbridge implemented this 
methodology in MSC/NASTRAN 6 . 

Yiu and Landess 30 also developed a similar methodology for forming a 
component which does not have an external applied load. However, their formulation 
is applicable to fixed interface components only. A criteria for concluding the 
recurrence sequence, based upon the rigid body mass and flexibility represented by 
the calculated static Ritz vectors, was proposed. 

Some applications of static Ritz (referred to in that article as Krylov) vectors, 
including unsymmetric, damped structural dynamics systems may be found in the 
work of Craig, Su, and Kim 11 . The focus of this effort is on the control of the flexible 
structure represented by the Krylov vectors. 

Both the Lanczos eigenvalue extraction algorithm and static Ritz vectors 
calculated using the boundary flexibility method are based upon the Krylov recurrence 
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sequence. As a result, some of the information and experience available in the 
published literature concerning Lanczos eigenvalue extraction has relevance to static 
Ritz vectors. Amongst the wide amount of available literature, the most complete and 
up to date source of information concerning the Krylov sequence and Lanczos 
eigenvalue extraction is Parlett’s book. The Symmetric Eigenvalue Problem 22 . This 
book also includes information on other eigenvalue extraction algorithms, as well as 
discussions on general linear algebra tools which are particularly useful in this type 
of analysis. 

In addition, two other excellent sources of information concerning the Lanczos 
eigenvalue extraction algorithm are the reports, A Shifted Block Lanczos Algorithm 
for Solving Sparse Symmetric Generalized Eigenproblems 13 , by Grimes, Lewis and 
Simon of Boeing Computer Services, and the MSC/NASTRAN Handbook for 
Numerical Methods 31 . These reports have a useful emphasis on the practical 
implementation of the Lanczos algorithm. The same implementation of the shifted 
block-Lanczos algorithm is presented in both of these documents. The first 
commercial implementation of Lanczos eigenvalue extraction was accomplished by 
Boeing Computer Services. MSC subsequently obtained this code from Boeing and 
implemented it in MSC/NASTRAN, and have since made a number of modifications 
and improvements to the algorithm. 
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2.2) Review of Boundary Flexibility Component Mode Synthesis: 

Almost all of the literature published on the variously titled, Krylov, Lanczos, 
or static Ritz, vectors focuses on the use of these vectors at the system level. These 
works deal with either Lanczos eigenvalue extraction, control dynamics, or external 
load derived Ritz vectors. The emphasis of this dissertation is the use of static Ritz 
vectors in component mode synthesis. As a result, a detailed review of the small 
amount of literature describing the static Ritz vector boundary flexibility method of 
component synthesis, which is the starting point of this work, is warranted. 

2.2.1) Fixed Interface Methodology: 

First, as is standard in component mode synthesis methods, the finite element 
component mass, m, and stiffness, k, matrices are partitioned into internal and external 
degrees of freedom, denoted by subscripts i and c, respectively. 


in 


cc 

Cl 

m 

771. 

1C 

u 

k 

k 

cc 

Cl 

k . 

k 

1C 

u 


The constraint modes are defined by 


= -t 

v ic K ii 


-1 


(2.3) 


(2.4) 


(2.5) 


which is the same definition used in standard component mode synthesis. 
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For Wilson’s method 29 , a set of externally applied loads is required to obtain 
the initial set of Ritz vectors. For the boundary flexibility method, this set of loads 
is created by multiplying the constraint modes by the mass matrix. (Craig 10 also 
included the off-diagonal mass matrix in his formulation.) Since the mass matrix is 
used to create the loads, they can be considered inertia loads. This set of inertia loads 
are then used to generate the initial set, or block, of Ritz vectors using the following 

ir-v'o*.®**" 1 *) (2 - 6) 


where the superscript ** indicates that the vectors in the matrix have not been 
normalized. The first block of vectors is normalized using the following equation, 
where the subscript r signifies that the block is normalized vector by vector. There 
are c vectors within each matrix, or block. 



r = l,2 v ..,c 


( 2 . 7 ) 


Note that the denominator of the right hand side of the equation is merely the inner 
product norm, ||^|| m , calculated with respect to the mass matrix. 

The subsequent sets of static Ritz vectors in the Krylov sequence are generated 
using the recurrence relationship which was defined in equation (2.2) 1929 , 

4j =k ti~ lm ti9j-i 


(2.8) 
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where the superscript * signifies that the vectors have not been orthogonalized or 
normalized. The additional sets of vectors are orthogonalized, with respect to the 
mass matrix, with all previous vectors. The process used to perform this 
orthogonalization is a normalized Gram-Schmidt procedure. 


where 


qj =q j -q lj -ic 


(2.9) 


c= 4ij-i m u9j 


( 2 . 10 ) 


and q, j.j is the concatenation of the previous sets of Ritz vectors. 


( 2 . 11 ) 


where all vectors have been normalized as follows. 


V 


q, 


l Jr 




r l,2,...,c 


( 2 . 12 ) 


The complete set of calculated Ritz vectors is included in the transformation 
matrix as Q t . (The resulting transformation matrix has the same form as that of 
"Craig-Bampton" component mode synthesis 9 , with the Ritz vectors replacing the 
normal modes.) 


Y = 


Icc 0 

K Qi 


(2.13) 
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The physical mass and stiffness matrices are transformed into the component modal 


matrices 

p = ¥ T m ¥ 


( 2 . 14 ) 


k = ¥ T k Y 


( 2 . 15 ) 


The resulting mass submatrices are 


where 



V-d 

V-u 


( 2 . 16 ) 


»cc = *ic T ( m ii®ic +m ic) +m ci*ic +m cc 

p fl = / fl = q] m ii Q l 


The resulting stiffness submatrices are 


where 



K cc = k ci*ic + k cc 


K fc = K w r = 0 


K a = Qi k »Qi 


( 2 . 17 ) 


( 2 . 18 ) 


( 2 . 19 ) 
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The use of constraint modes in the transformation matrix leads to the null off-diagonal 
partitions of the component stiffness matrix, just as in the approach based upon normal 
modes. 

2.2.2) Free Interface Methodology for Components with No Rigid Body Modes: 

When allowing the interface points of a component to be free to deflect while 
forming the component, a somewhat different basis for the initial vector of the Krylov 
sequence is required. The attachment modes, rather than the constraint modes are 
utilized in initial block definition. By definition, the attachment modes are the 
columns of the flexibility matrix which correspond to the interface degrees of freedom. 

g = r 1 (2.20) 
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( 2 . 21 ) 


The initial block of vectors in the free interface formulation is defined as 


q i** = k' l mg a 


(2.22) 


and is normalized as follows. 
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(2.23) 
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Note that the unpartitioned physical mass and stiffness matrices of the component are 
used in the free interface formulation. The recurrence algorithm then proceeds in the 
same manner as in the fixed interface methodology. 


% 


q/ = * m 9j. 1 

(2.24) 


(2.25) 

c= 9ij-i Tm 9/ 

(2.26) 

q ** 

' r r = 1 

(2.27) 




Formation of a static Ritz vector component then follows the normal 
component mode synthesis techniques which were presented by MacNeal 17 and 
Rubin 24 . To combine the "Rubin-MacNeal" method with the presented method, the 
normal eigenvectors are simply replaced with the static Ritz vectors, as in the fixed 
interface methodology. The free interface methodology uses residual flexibility terms, 
which fully define the stiffness missing from the modal space due to excluded modes. 
The flexibility contained in the calculated Ritz vectors is given by the following 
equation. 

g * = <?,((?/* C ?/)' 1 ^ 7 (2 - 28) 


The unrepresented flexibility, or residual flexibility, is defined as 


gd = S~g k 


(2.29) 
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The residual flexibility matrix is then partitioned in the same manner as the flexibility 
matrix was in equation (2.21), when the attachment modes were created for initial 
vector calculation. The result is the residual attachment modes. 


g 


a i 


Scc d 

Si Cj 


( 2 . 30 ) 


When the residual attachment modes, g^, are added to the Ritz vectors, Q h the 
complete flexibility of the component is represented. 

The residual attachment modes and the Ritz vectors are used to form the 
component transformation matrix. This matrix transforms the physical subspace, u, 
to the modal subspace, p, and is defined by the following equation. 
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( 2 . 31 ) 


In order to provide physical interface degrees of freedom, for use in component 
coupling, p c in the above equation is back-transformed to eliminate it from the right- 
hand side of the equation. This results in the following transformation matrix, 
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( 2 . 32 ) 
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where 

g. * = g. g ~ l 

* 1C O IC rf O cc rf 

<?«*=<?« -^^cc; 1 ^ (2 - 33) 

The transformation of the component mass and stiffness matrices then proceeds in a 
similar manner as shown in equations (2.14) through (2.19), with the following 
differences. The Q it * matrix partition replaces the Q, matrix partition. The g u * matrix 
partition replaces the <I> ic matrix partition. In the fixed-interface methodology, the 
definition of the constraint modes, <E> ic , leads to terms in the component stiffness 
matrix which cancel out. In the free-interface methodology, the definition of the 
transformation submatrices has changed and so this cancellation does not occur. 
Therefore, the corresponding equations in (2.19) are replaced by the following 
equations, respectively. 

\c=8ic' T ( k iigic' +k ic) +k ci8ic* + k cc 

K jc - = Qu r (Mic* + *ic) (234) 

2.2.3) Free Interface Methodology for Components with Rigid Body Modes: 

When a component has rigid body modes, the associated stiffness matrix is 
singular. The inverse of the stiffness matrix, the flexibility matrix, cannot be directly 
obtained, and therefore the attachment modes cannot be directly obtained. To 
circumvent this problem, Rubin 24 presented the following method for obtaining the 
residual elastic attachment modes of a component with rigid body modes. 
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First, the stiffness matrix is constrained from rigid body motion by partitioning 
out r degrees of freedom, where r is the number of rigid body modes. 



(2.35) 



(2.36) 


(2.37) 


A square projection matrix is defined by 

A=I m -m<t> r <S> T T (2.38) 

where <I> r is the rigid body modes matrix. The elastic flexibility matrix, g e , with rigid 
body motion removed, is shown in reference [24] to be 

g e = A T g c A ( 2 . 39 ) 


Now the analysis proceeds in a similar fashion to the previously discussed 
methodology of the free interface component with no rigid body motion. The major 
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difference between the two approaches is that the elastic flexibility matrix is used in 
place of the general flexibility matrix. The inertia relief attachment modes are 

(2.40) 



The initial block of vectors is calculated using the inertia relief attachment modes and 
the elastic flexibility matrix. 


(2.41) 


The subsequent static Ritz vectors are calculated, orthogonalized, and 
normalized as shown in equations (2.24) through (2.27). The residual elastic 
flexibility terms are also calculated as shown in the free interface with no rigid body 
modes discussion, equations (2.28) through (2.30). Creation of the transformation 
matrix, equations (2.31) through (2.34), is also similar to when no rigid body modes 
are present. The one exception is that the rigid body modes must be included in the 
transformation matrix. Therefore, equation (2.31) is replaced by 



e o . $ 

o cc d ^ cl cr 

2 O 7 ^ 

& ic d \r 



(2.42) 


Formation of the final transformation matrix, and subsequently the component mass 
and stiffness matrices, is then performed as described in the previous section. 
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2.3) Relationship to Lanczos Eigenvalue Extraction: 

In addition to being the result of the fundamental recurrence equation when 
calculating static Ritz vectors, a Krylov sequence is also the result of the fundamental 
recurrence relation used in the Lanczos eigenvalue extraction algorithm 15 22 . The 
Lanczos eigenvalue extraction algorithm uses a Krylov sequence to generate terms 
which can be used to transform the given system into a tridiagonal form. This 
tridiagonal matrix is then diagonalized using a QR, or related, eigenvalue extraction 
algorithm and the resulting diagonal terms are the eigenvalues of the original system. 
A detailed presentation of the Lanczos algorithm is beyond the scope of this work. 
Since Krylov vectors are the foundation of the Lanczos eigenvalue extraction 
algorithm, some of the existing literature which investigates orthogonalization, shifting, 
and practical implementation of the Lanczos algorithm is applicable to static Ritz 
vectors, since they are also derived from the Krylov sequence. 

There are also significant differences between the Lanczos eigenvalue 
extraction algorithm 1522 and the use of static Ritz vectors directly. These differences 
are summarized as follows. In the selection of the initial seed to start the Krylov 
sequence, Lanczos starts with a random vector. This is done to prevent the Krylov 
sequence from converging to a particular class of eigenvectors and skipping other 
eigenvectors entirely. Since another transformation will take place converting the 
tridiagonal matrix into the system eigenvalues, the vectors from the sequence need not 
have a particular physical significance. Wilson’s methodology initiates the sequence 
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with a matrix which has physical significance, i.e., a load vector. In the boundary 
flexibility method that load vector is the constraint or attachment modes multiplied by 
the component mass matrix, forming an inertia load vector block. That leads the 
vectors which result from subsequent iterations of the sequence to tend to also have 
physical significance. This difference alters the Krylov sequence from being a purely 
mathematical tool for eigenvalue extraction, into a mechanically based application in 
structural dynamics. 

These differences between the two methodologies stem from the differing goals 
of the two algorithms. As discussed above, the Lanczos algorithm searches for the 
eigensystem by first transforming the original system matrices into tridiagonal form, 
and then diagonalizing that tridiagonal form to obtain the eigenvalues. Some 
applications of Wilson’s methodology make use of the tridiagonal form, but the 
orthonormalized vectors obtained from the sequence are instead used to transform and 
reduce the original system directly. The eigenvalues and eigenvectors are not 
obtained, and both transformed mass and stiffness matrices do not assume a diagonal 
form. In all component mode synthesis the reduced matrices are not diagonal anyway, 
due to the coupling mass and stiffness partitions of the final reduced matrices 
(equations 2.17 and 2.34). Therefore, the non-diagonal form obtained from the use 
of static ritz vectors is not a disadvantage when used in component mode synthesis. 
In addition, the vector space in which a transformed system is tridiagonal is not the 
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same as in the component mode synthesis algorithm, and, as a result, is not useful as 
a final goal. 

Vector blocks are also utilized in the practical implementation of the Lanczos 
eigenvalue extraction algorithms 1 3,22-3 '. The original Lanczos algorithm had difficulty 
in determining a complete set of multiple eigenvalues. Using vector blocks in the 
Krylov sequence allows the algorithm to determine multiple roots, up to the dimension 
of the block. When blocks are used in the Krylov sequence, the tridiagonal form of 
the transformed system also assumes a blocked format. The bandwidth of the 
tridiagonal form is then determined by the dimension of the block. The block format 
complicates orthogonalization and sequence truncation schemes. The block format 
is a natural feature in the boundary flexibility method because of the multidimensional 
inertia load vector block. 

2.41 Orthogonalization: 

The most computationally expensive aspect of the formulation of static Ritz 
vectors is the process of orthogonalization (equations (2.9) through (2.12)). As a 
result, efficient orthogonalization is critical in determining this method’s efficiency 
when compared to the use of eigenvectors. Since Ritz vectors are not inherently 
independent, for them to be used in a similarity transformation, they must be 
orthogonalized. Obviously if this process is more computationally expensive than the 
calculation of eigenvectors, then it’s usefulness is limited. 
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2.4.1) Loss of Orthogonality: 

It has been shown in several works' 519 ' 22 , that if a Ritz vector obtained from 
the Krylov sequence is orthogonalized, using Gram-Schmidt, with the two previous 
Krylov vectors, it is theoretically orthogonal to all previously generated vectors. The 
orthonormalizing coefficients are assembled into a tridiagonal matrix. The assembled 
tridiagonal matrix is also the result of the Krylov coordinate transformation on the 
original system matrix (Q T AQ). In other words, use of the Krylov sequence and the 
properties of orthogonality, allow the system to also assume a tridiagonal form. 
Unfortunately, in practice the theoretical orthogonality that each new Krylov vector 
has with all vectors, after orthogonalization with the previous two, is usually lost due 
to either even minimal computational round-off error, or other factors which will be 
discussed in chapter 3. A mathematically rigorous explanation of this phenomena is 
given in references [21], [22] and [26]. As a result, additional orthogonalization and 
sometimes re-orthogonalization is required in order to perform a correct transformation 
and to maintain the tridiagonal form. 

A brief clarification of terms found in the existing literature and used in this 
work follows. When Gram-Schmidt procedures are required, they are sometimes 
referred to as orthogonalization and sometimes as re-orthogonalization, varying with 
author. In this work, orthogonalization refers to the initial Gram-Schmidt process, 
even if performed with all previous vectors. Re-orthogonalization refers to any 
additional orthogonalization steps following the initial orthogonalization. 
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The initial solution to the loss of orthogonality, which occurs in the Krylov 
sequence, was to simply explicitly orthogonalize with all previous vectors, rather than 
just the previous two. This orthogonalization scheme is computationally more 
expensive than orthogonalizing with the previous two vectors, especially as the 
number of previously determined vectors becomes large. The computational expense 
of orthogonalization is what limited Lanczos eigenvalue extraction use for many years. 
Many enhancements to the basic Lanczos eigenvalue extraction algorithm have 
appeared in the literature through the years, some of which were directed toward 
insuring orthogonality. Paige 21 established a theorem, using matrix norms and terms 
from the tridiagonal matrix, which yields a numerical criterion for determining when 
re-orthogonalization is required. The orthogonality of the obtained vectors is not 
explicitly calculated. When the requirement for orthogonalization does arise, the new 
vector would be orthogonalized with respect to all previous vectors. Parlett and 
Scott 23 also use a numerical criterion, derived using matrix norms and terms from the 
tridiagonal matrix, to determine when re-orthogonalization is required. They proposed 
the simple modification of orthogonalizing, using modified Gram-Schmidt 4 , with only 
those previous vectors with which the new vector is not orthogonal. Simon 26 clarified 
issues dealing with the loss of orthogonality, developed an additional re- 
orthogonalization scheme, and investigated the complete orthogonalization issue in 


depth. 
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2.4.2 Orthonormalization: 

In general, the existing literature describes the Lanczos eigenvalue algorithm 
as extracting the eigenvalues of a single symmetric matrix, A . 

Ax = Xx (2-43) 

The vectors obtained from the Krylov sequence are orthonormalized using the 
Euclidean vector norm, which is defined as ||x || 2 = ( |x, I 2 + |x 2 | 2 + ... + |x p | 2 ) 1/2 . 
Because of the spectral theorem 22 , the eigenvectors are also orthonormal with respect 
to the A matrix. However, as discussed above, in structural dynamics, the 
eigenproblem which is being solved is a system with both mass and stiffness matrices. 

Kx = XMx ( 2 -44) 

In the MSC/NASTRAN application of the Lanczos eigenvalue extraction algorithm, 
the mass and stiffness matrices are used directly in the Krylov sequence and the 
Krylov vectors are orthonormalized with respect to the mass matrix 31 . 

In the use of static Ritz vectors for structural dynamics applications, as 
described in detail in section 2.2, the mass and stiffness are used in the Krylov 
sequence and the vectors are orthonormalized with respect to the mass matrix. There 
is one exception to the use of the mass matrix. In the work of Su and Craig 27 , the 
static Ritz (referred to as Krylov) vectors are orthonormalized with respect to the 
stiffness matrix. The result of orthonormalizing with respect to the mass matrix will 
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be that the modal partition of the reduced system mass matrix is the identity matrix. 
One can orthonormalize with respect to the mass matrix, the stiffness matrix, or using 
the Euclidean vector norm, but not all three, with a single Gram-Schmidt 
orthogonalization procedure. Orthonormalizing with respect to mass (or stiffness) 
allows some flexibility in dealing with the reduced mass matrix in some applications 
and is required to generate the tridiagonal form when using the Krylov sequence of 
equation (2.2). This is a direct result of the matrix product, K ' l M, not being 
symmetric. As previously discussed, in component mode synthesis the resulting 
transformed mass and stiffness matrices are, by definition, not diagonal, and exist in 
a different space than that of the tridiagonal form. Therefore, the matrix form of the 
transformed component matrices is not an important issue, and it is not necessary to 
maintain the ability to generate the tridiagonal form. 

2.4.3) Choleskv/OR Decomposition: 

Parlett 22 presents the following discussion, relating to orthogonalization. Any 
non-null rectangular m by n matrix B can be written as B = QR with m by r Q 
satisfying Q T Q = I r , and r by n R upper triangular with non-negative diagonal 
elements. The QR factorization is the matrix formulation of the Gram-Schmidt 
process for orthonormalizing the columns of B. When B has full rank, then R is the 
upper Cholesky factor of B T B since 


R T R=R T Q T QR=B T B 


(2.45) 
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Use of the upper Cholesky factor of a matrix for orthogonalization will be presented 
in section 3.3.5. 

2.5) Shifting: 

A technique defined by Scott 25 shows how an inverse operator, either applied 
explicitly while using subspace iteration or applied implicitly by while using the 
Lanczos algorithm, can direct a solution to particular frequency range. The use of the 
inverse operator in the Lanczos eigenvalue extraction algorithm is commonly called 
shifting. Shifting has been implemented with success in the commercially available 
Lanczos eigenvalue extraction routine 1331 . The inverse operator is applied to the 
Krylov sequence within the Lanczos algorithm, and therefore can be applicable to any 
Krylov sequence based solution. 

In the Scott paper, the problem of computing some eigenpairs of the 
generalized eigenvalue problem is considered, 

(A-XM)x = 0 ( 2 - 46 ) 

with X being the eigenvalues and x the eigenvectors of the pencil (A,M). A system 
is created using the operator 


( A-oMY l M 


that has the same eigenvectors as (2.46). The shifted system eigenvalues are 
transformed to 1 /(k - a). This means that the eigenvalue nearest <7 becomes the 
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dominant eigenvalue and the sequence will converge to the corresponding eigenvector. 
In the commercially available Lanczos routines the algorithm is applied to the shifted 
and inverted eigenvalue problem of the following equation 

M(A - oM)~ 1 Mx = — —Mx (2.47) 

X-a 

The use of the shifted eigenvalue form allows the for good approximations to 
eigenvalues within specific ranges, even if they are clustered. The cost for having the 
advantage of shifting is the factorization of ( A - oM)' 1 . 

2.6) Spectrum Slicing: 

Parlett 22 presents the following theorem. When the triangular factorization of 
(A - gM) is calculated, if A is symmetric then 

(A-oM) =LDL t ( 2 - 48 ) 

where D is diagonal and M is positive definite. Then 

v(A-o/)=v(A-oM)=v(D) (2-49) 

where v is number of negative eigenvalues and A = diag (^, Xi , , K )• Th e 

number of negative elements of D equals the number of eigenvalues of the pencil 
(A,M) which are less than G. As a result, whenever a shift is undertaken the number 
of eigenvalues below the shift frequency is determined. 
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2.7) Krylov Sequence Termination Techniques: 

A large uncertainty in the use of static Ritz vectors and the boundary flexibility 
method is the lack of a criteria for terminating the Krylov sequence. When 
eigenvectors are used to form a component, typically a frequency range is defined and 
all eigenvectors with eigenvalues within that given range are determined. That option 
does not exist in the use of static Ritz vectors. A number of error criterions have been 
proposed which truncate the sequence when a somewhat arbitrary variable reaches a 
arbitrary value. 

Wilson, et. al. 29 , used a definition of the modal participation factor to define 
an error term. This factor is equivalent to the dependence coefficient in a Gram- 
Schmidt orthogonalization procedure. The error term is the linearly independent 
portion of the force vector, with respect to the Ritz vectors. The linearly independent 
portion of the force vector should be zero if the complete static solution is desired. 
No account is made for the dynamic response of the system in the error term. A 
similar error term was defined by Leger 16 . 

Nour-Omid and Clough also used a modal participation factor as a sequence 
truncation criterion 19 20 . However, they did not define or use an error term based upon 
the desired dependence of the force vector. They proposed a simple cut off when the 
mode participation factor (or dependence coefficient) reached an arbitrary numerical 
value. No account is made for the dynamic response of the system. 
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Su and Craig 27 stated that a modal participation factor, such as used by Nour- 
Omid and Cough, was not valid for a non-diagonal representations, which are the 
result of the representations not being based upon normal modes. Therefore, they 
proposed that the norm of the off-diagonal submatrices in the transformed mass and 
stiffness matrices should be used as sequence truncation criterion. Again, the 
sequence would be terminated when the norm reached an arbitrary numerical value. 

Yiu and Landess 30 proposed two similar sequence truncation criteria. The first 
used a flexibility convergence criteria, similar to the error term proposed by Wilson, 
et. al. A mass convergence criteria was also proposed. It was based upon the amount 
of the rigid body mass, rather than the force vector, which would be represented by 
Ritz vectors. The percentage of rigid body mass represented by Ritz vectors is 
commonly referred to as effective mass 14 . When the percentage of rigid body mass 
represented reaches an arbitrary numerical value the sequence is terminated. 



Chapter 3 


Theoretical Development 

3.11 Introduction: 

This short chapter presents two theorems which pertain to the use of static Ritz 
vectors in the boundary flexibility method of component mode synthesis. These 
theorems concern whether or not static Ritz vectors can be used to represent a 
component in a mathematically rigorous fashion. Proofs are presented which 
demonstrate the fact that static Ritz vectors can be used to correctly represent a 
component. 

3.21 The Exact Nature of the Methodology: 

Currently, most component mode synthesis applications use the normal 
eigenvectors of the substructure to form the component. If all the eigenvectors of a 
system are used to form the component, the complete dynamic properties of the 
component are represented and an "exact" finite element solution may be obtained. 
This is because, if all the modes are used, the component representation in not a Ritz 
vector approximation but is instead a complete linear coordinate transformation. The 
same principle holds true for components based upon a block-Krylov sequence. It 
was proven in references [10] and [27] that an n size Krylov subspace spans the entire 
n dimensional space. The following theorem is an alternative demonstration of the 
proposition that, if n Krylov vectors are used to form a component from a system of 
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size n, then the same complete solutions, as found directly from finite elements or 
normal modes, can be obtained. 

Theorem 3.1 - Given that * is an eigenvector of A, and that A is the 
corresponding eigenvalue, then for B, where B = F'AP, F'x is the associated 
eigenvector of B and A is the invariant eigenvalue of both B and A . 

Proof - The eigensystem of A is defined as 

Ax = X x (3-D 


and since, from the definition of B, 

A = PBP 1 


equation (3.1) can be re-written as 

PBP~ x x = Xx 


Premultiplying by P 1 yields 

BP 1 x = XP 1 x 


(3.2) 


(3.3) 


(3.4) 


If y is defined as P ! x, then 


By = Xy 


(3.5) 
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and it is evident that the eigenvalue of this system is A, and that the associated 
eigenvector is P'x. ■ 

It is important to note that the only constraint on P, implicit or explicit, is that 
it is an nonsingular matrix, with the same dimension as A. For P to be nonsingular 
it must be of full rank, or equivalently stated, its columns must be linearly 
independent. Therefore, in order for an exact coordinate transformation to be 
accomplished, the vectors from which the transformation matrix is assembled must be 
linearly independent. These vectors do not need to be the eigenvectors, and they do 
not need to be orthogonal. 

The fact that the vectors in a transformation matrix need only be linearly 
independent, and not orthogonal, is already utilized in component synthesis based 
upon normal eigenvectors. The transformation matrices which result from normal 
mode component synthesis consist of linearly independent vectors, not orthogonal 
ones. This can be demonstrated simply by inspecting the result of the matrix 
transformation, as shown in equations (2.14) through (2.19). The complete reduced 
mass and stiffness matrices are not diagonal. Only the modal partition of the matrices, 
which does result from an orthogonal transformation, is diagonal. Therefore, even 
when compared to normal eigenvectors, the use of static Ritz vectors contains no 
inherent disadvantages of matrix form or accuracy, since the resulting complete 
component mass and stiffness matrices in either case are not diagonal. 
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In the proof of theorem 3.1, the inverse of the transformation matrix was used 
in the pre-multiply position. Typically, as in equations (2.14) and (2.15), the transpose 
of the transformation matrix is used in this position. If a set of orthogonal vectors 
make up a transformation matrix, then the inverse and the transpose of the matrix are 
identical. If the vectors are only linearly independent, and not orthogonal, then the 
transpose and the inverse of the transformation matrix are not identical. However, in 
structural dynamics, this does not necessitate using the inverse of the transformation 
matrix. In structured dynamics, the eigensystem represents the space of the matrix 
multiplication M'K, which is derived from Kx = AMx. The following theorem shows 
that for the M'K space, use of the inverted transformation matrix and the transposed 
transformation matrix is interchangeable. 

Theorem 3.2 - Specifying that P 7 ?P' , so that P is not orthonormal, and if M a 
= P r MP and K a = P r KP, which represent a transformation using the transpose, and 
M h = P'MP and K h =P'KP, which represents a transformation using the inverse, then 
M;'K a = M h 'K h . 

Proof - This may be proven by substitution. The proposition is that 

M' l K = M b l K b (3-6) 

a a oo 

and since, by definition, 

M=P t MP M h = P l MP 

a D 

K=P t KP K h = P l KP 

a D 


(3.7) 
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These equations may be substituted into equation (3.6), yielding 

(P t MP) 1 P t KP = 0 P' 1 MP) 1 P 1 KP ( 3 8 ) 

Post-multiplying by P‘K ‘P , and canceling, results in 

(p t mp) 1 p t p = (p-'Mpy 1 ( 3 - 9 ) 

and then post-multiplying by P'MP yields, 

C p t mp) 1 p t mp=i ( 3 - 10 ) 

which reduces to 

/ = / B (3-11) 

In summary, it was proven in theorem 3.1 that orthogonal vectors are not 
required for an exact transformation, i.e., linearly independent vectors suffice. This 
proof used the inverse of the transformation matrix in the pre-multiply position. With 
a non-orthogonal transformation matrix, the transpose and the inverse of the 
transformation matrix are, by definition, not identical. As a result, a transformation 
which uses the transpose of a non-orthogonal matrix results in a mass and stiffness 
matrix different from the result obtained from a transformation using the inverse of 
the transformation matrix in the pre-multiply position. In theorem 3.2 it was shown 
that when the transformed mass and stiffness matrices, which result from the use of 



38 


the non-orthogonal transpose, are combined, the result is identical to the result 
obtained from using the inverse in the pre-multiply position. Obviously, the 
eigensystem of the two transformation results are then also identical. Therefore, if the 
dimension of the transformation matrix is equal to the dimension of the component 
matrix, then Ritz vectors, or any linearly independent set of vectors, can form an exact 
transformation. Also, if the system is a dual matrix system, such as one consisting 
of a mass and stiffness matrix, then the transpose of the linearly independent set of 
vectors may be used in the pre-multiply position of the transformation. 
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Chapter 4 

Orthogonalization, the Krylov Sequence and Static Ritz Vectors 
4.1) Introduction: 

This chapter examines several related issues pertaining to the orthogonalization 
of static Ritz vectors in the boundary flexibility method of component mode synthesis. 
Efficient orthogonalization is essential to the efficiency of the methodology because 
the bulk of the computing effort required to produce a static Ritz vector component 
is in vector orthogonalization. Section 4.2 examines general issues concerning 
orthogonalization of the vectors obtained from the Krylov sequence. The vectors 
derived from this sequence must be orthogonalized with respect to previously obtained 
vectors, usually using a Gram-Schmidt approach, to insure linear independence in the 
transformation matrix. Various orthogonalization schemes are proposed and examined 
for their accuracy, robustness, and efficiency. In section 4.3 issues concerning the use 
of blocks with the Krylov sequence are discussed. These topics include orthogonality 
within the block, and reducing the block to a manageable size. Algorithms containing 
the new orthogonalization schemes are presented in section 4.4 
4.21 Orthogonalization: 

As presented in chapter 2, the static Ritz boundary flexibility method of 
component synthesis is based upon the Krylov sequence. The vectors derived from 
this sequence must be orthogonalized with respect to previously obtained vectors to 
insure linear independence in the transformation matrix. If linear independence, or 
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orthogonality is not maintained, numerical errors in the transformation are a certainty. 
In addition, for a system of size n, the orthogonalization methodology should be 
robust enough to obtain n linearly independent vectors, to insure completeness in the 
representation. Orthogonalization routines such as Gram-Schmidt or modified Gram- 
Schmidt can be extremely expensive computationally if implemented inefficiently. 
Lanczos eigenvalue extraction was not used widely until efficient orthogonalization 
schemes were implemented within the algorithm. Consistent with this, for static Ritz 
vectors to be practical, the orthogonalization scheme must be accurate, robust, and 
efficient. 

Issues examined in this section include the reason why the loss of 
orthogonality occurs, checking of vector orthogonality, various Gram-Schmidt 
orthogonalization schemes, orthonormalization options, and alternate orthogonalization 
methodologies. An appropriate, workable orthogonalization scheme is suggested. 
This proposed scheme is contrasted with those used with Lanczos eigenvalue 
extraction. In Lanczos eigenvalue extraction, it is also important to maintain 
orthogonality to insure linearly independent vectors. In addition, maintaining the 
orthogonality also aids in determining multiple eigenvalues, as well as maintaining the 
tridiagonal form 22 . As discussed in chapter 2, maintaining the tridiagonal form is not 
important in the use of static Ritz vectors, when used in the boundary flexibility 
method of component synthesis. 



41 


4.2.1) Linear Independence and the Loss of Orthogonality: 

It was shown in chapter 3.2 that for an n size system, n independent vectors 
will form a transformation matrix which will allow an exact linear transformation. 
The boundary flexibility algorithms presented in chapter 2 would not be successful in 
obtaining n independent vectors without additional enhancements. Specifically, the 
normalized Gram-Schmidt orthogonalization technique utilized in these algorithms is 
not robust enough to obtain independent vectors which span the entire n space. It 
must be said that, in typical applications, there is no requirement for the entire space 
of size n to be represented. One major advantage of component mode synthesis is a 
reduction of system size. However, in some cases this loss of orthogonality begins 
quite early within the Krylov sequence. As a result, no guarantee can be made that 
a sufficient number of vectors, to adequately represent the component, is obtainable. 
In addition, whether the component is reduced by means of sequence truncation, 
vector selection, or any other approach, a correct reduction can not be guaranteed, 
unless the entire dynamic space is obtainable. Therefore, the static Ritz vector 
algorithm should be able to yield n independent vectors for an n size system. 

Theoretically, each Krylov vector is linearly independent of previous vectors. 
This can be demonstrated by inspection of equation (2.1). Since each vector is a 
product of multiplication of previous vectors, except for the special cases of null or 
unity spaces, each vector can not also be defined as a linear combination of previous 
vectors. However, in practice, numerical dependence does occur in vectors obtained 
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from a Krylov sequence. Due to this numerical dependence, n independent vectors, 
describing an n size system, can not be directly obtained from the Krylov sequence. 
Understanding of this problem can be aided by considering a Krylov sequence based 
upon a space described by the normalized vector, x 2 . Assume, in equation (2.1), that 
A, the iteration matrix, is x 2 and that <j>, the initialing vector, is e„ the unit column 
matrix of order one. The Krylov sequence becomes 


(4.1) 


Figure (4.1) shows a plot of the resulting vectors, to the ninth order, normalized to 
unity. In other words, the sequence of functions x 2 , x 4 , x 6 , x 8 , x’°, x 12 , x’ 4 and x 16 , 
have been plotted. The unit column matrix, e,, was not plotted. It can be seen that 
as this series continues to a higher order, the vectors become nearly dependent 
because of computational roundoff error, and eventually, numerically indistinguishable, 
despite theoretical independence. The Gram-Schmidt algorithm is not able to 
orthogonalize a vector which is numerically dependent on previous vectors. 

An analytical source of vector dependency, not dependent on computational 
roundoff, is the situation when the Krylov sequence converges to a normal 
eigenvector. This state might seem to be desirable, considering the traditional use of 
eigenvectors in dynamic analysis and component mode representation. (In this 
circumstance the Krylov sequence acts similar to a power method of eigenvalue 
extraction.) Unfortunately, when the Krylov vectors have converged and produced 
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an eigenvector, the next vector in the Krylov sequence is the same eigenvector, scaled, 
which is linearly dependent. Theorem 4. 1 demonstrates the convergence. 

Theorem 4.1 - If an eigenvector of the system Kx = AMx, x t , appears in the 
Krylov sequence, x i+ , = K 'Mx,, then the resulting vector, x M , is linearly dependent 
upon x t , and differs only by a scale factor of 7//L 
Proof - If 


and since the Krylov sequence is defined as 

K'Mx, 

then 

-x i =K- l Mx i 
k ' 

which leads to 

Kx i = XMx i 


(4.2) 


(4.3) 


(4-4) 


(4.5) 


which is true, from the definition of an eigenvector. ■ 

The potential convergence of the Krylov sequence to system eigenvectors has 
several implications. It is a reason, in addition to not prejudicing the sequence 
towards a certain eigenvector, that in Lanczos eigenvalue extraction, the Krylov 
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sequence should be initiated by a random vector. In Lanczos eigenvalue extraction, 
the system will be re-transformed to a diagonal form and so the eigenvectors are not 
required at that step. The potential convergence of the Krylov sequence to 
eigenvectors also restricts shifting of the sequence 23 31 . Shifting the sequence so that 
it would produce a vector close to an eigenvector would cause the dependency 
problem discussed above. 

In the use of generalized static Ritz vectors, where the vectors which initiate 
the sequence and the vectors resulting from the sequence have physical significance, 
the convergence of the Krylov sequence to eigenvectors can cause problems of linear 
dependency. Consider a cantilevered beam under a gravity load. The deflected shape 
of the beam is very close to the first normal eigenvector. (A cantilevered beam under 
a gravity load is used as an example and is illustrated in chapter 5.) If that vector 
initiates the sequence, subsequent vectors will be nearly linearly dependent on the 
previous vector. The greater the similarity a static Ritz vector has to an eigenvector, 
the greater will be the dependency of the subsequent vector in the Krylov sequence. 
In practice, total dependency does not occur on digital computers due to the same 
class of roundoff errors that lead to the requirement for re-orthogonalization discussed 
in section 2.4. 1 . When two Ritz vectors are nearly linearly dependent, the differences 
between the two vectors will tend to be random, and so physically significant 
information that the vector will contribute, after orthogonalization, would be minimal. 
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This random element resulting from orthogonalization does have some benefit. 
In Parlett 22 , the process of randomization is discussed. There, randomization is 
defined as a process of creating a set of random vectors, after several blocks of static 
Ritz vectors have been obtained, and then orthogonalizing these random vectors with 
the previous static Ritz vectors. A Krylov sequence initiated with a vector orthogonal 
to a space requiring representation may have difficulty in producing n independent 
vectors for a system of size n. For example, if the initial vector does not contain a 
displacement in one of the three ordinate directions, then a pure Krylov sequence 
would not generate a vector representing the system in that direction. Introducing a 
random element into all degrees of freedom contained in a vector, and then 
orthogonalizing, would allow all possible directions and shapes to be represented. 

It has been found that due to computational roundoff error an explicit 
randomization routine is not necessarily required. As discussed above, when a vector 
is orthogonal ized with a vector upon which it is nearly dependent, the purified vector 
will contain a random element. As a result, randomization occurs to some degree in 
all Krylov processes implemented on a digital computer. In this manner, vectors may 
be obtained which are orthogonal to the initiating vector in the Krylov sequence and 
the complete n size component space may be spanned. 

4.2.2) Use of The Euclidean Norm for Normalization: 

In section 2.4.2, it was discussed that in the creation of static Ritz vectors, 
where the Krylov sequence of equation (2.2) is used to generate the vectors, previous 
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Gram-Schmidt algorithms orthonormalized the vectors with respect to the mass matrix. 
This is required if obtaining a tridiagonal form is desired. If the Krylov sequence is 
based upon equation (2.1), then the Euclidean norm is used. In the implementation 
of the Gram-Schmidt orthogonalization procedure, it was found that the use of the 
mass matrix in orthonormalization becomes comparatively and extremely costly for 
anything other than a small problem size. As a result, alternatives were examined. 
The tridiagonal form is not a particular advantage using Ritz vectors, as opposed to 
eigenvectors, in component mode synthesis. In component mode synthesis, by 
definition, the transformed matrices have large off-diagonal components. Therefore, 
it is possible to orthonormalize using the Euclidean norm, even though doing this will 
not produce a diagonal modal mass matrix and a tridiagonal modal stiffness matrix. 
Again, this feature is not a disadvantage in component mode synthesis and very large 
savings in computational cost are achievable, as documented in chapter 5. 

4.2.3) Gram-Schmidt Failure and Reorthoeonalization: 

The simple normalized Gram-Schmidt orthonormalization procedure outlined 
in section 2.2 is inadequate to guarantee a set of orthogonal and linearly independent 
set of Ritz vectors. In the previous section, it was discussed how linear dependence 
can arise in vectors generated by a Krylov sequence. It is well documented that the 
Gram-Schmidt orthogonalization algorithm is not successful at producing an 
orthogonal vector from a nearly dependent vector 18 . The Gram-Schmidt procedure 
will fail on occasions when vectors, although theoretically independent, are dependent 
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within the numerical constraints of current digital computers. These vectors cannot 
then be made orthogonal, using the single normalized Gram-Schmidt step described 
in the chapter 2. 

One option for orthogonalizing numerically nearly dependent vectors is the 
modified Gram-Schmidt algorithm. Modified Gram-Schmidt is a computationally 
expensive procedure which is very successful at orthogonalizing nearly dependent 
vectors 418 . Implementation of modified Gram-Schmidt will be discussed in the next 
section and its computational expense will be discussed in the next chapter. A 
potentially less expensive option is to repeat the Gram-Schmidt orthonormalization, 
if it has been unsuccessful in the first attempt. Reorthogonalization will work, even 
though the first Gram-Schmidt attempt has been unsuccessful at producing an 
orthogonal vector, because it modifies the vector enough so that it is no longer 
numerically nearly dependent. The second normalized Gram-Schmidt step, or 
reorthogonalization, is then usually successful at producing a orthogonal vector. Two 
normalized Gram-Schmidt orthogonalizations will typically be less expensive than one 
modified Gram-Schmidt orthogonalization, especially considering that the 
reorthogonalization step will not be required for every Krylov vector. Computational 
costs comparisons for various models will be presented and discussed in chapter 5. 

If Gram-Schmidt orthonormalization is to be repeated when unsuccessful, a 
assessment of the Ritz vectors’ orthogonality is required. Upon completion of each 
Krylov vector calculation, and the associated normalized Gram-Schmidt 
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orthogonalization, each vector is checked for orthogonality. The orthogonality of the 
block-Krylov vectors is checked using one of the following equations; either 




(4.6) 


or 




= «U-1 



(4.7) 


depending on whether the vectors have been orthogonalized with respect to the mass 
matrix or using the Euclidean vector norm. (The mass matrix used in equation (4.6) 
is appropriate for fixed interface modes. For the free interface approach, the complete 
physical mass matrix is used.) If the new vector, q p is orthogonal to all previously 
calculated vectors, L m or L, will be a j by 1 size null vector. The infinity norm of the 
L vector, which is defined as 

||L|L = max{|/ 1 |,|/ 2 |,...,|/ 7 .|} (4-8) 


is then calculated and compared to a specified value, e. If L„ > e, then the associated 
vector is judged to be non-orthogonal and the Gram-Schmidt algorithm is repeated. 
4.2.41 Various Gram-Schmidt Orthogonalization Strategies: 

A number of combinations of Gram-Schmidt orthonormalization and 
reorthogonalization strategies are possible. These strategies, for the initial Gram- 
Schmidt step, include complete orthogonalization and orthogonalization with the 
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previous two blocks only, a possibility which was discussed in section 2.4. It is 
possible to use modified Gram-Schmidt exclusively, however that option is not 
competitive computationally. Complete reorthogonalization, selective 
reorthogonalization, and selective reorthogonalization using modified Gram-Schmidt 
were the reorthogonalization options investigated. Orthonormalizing with respect to 
the mass matrix or using the Euclidean vector norm can be performed with any of the 
above possibilities. The total number of possibilities investigated, amongst the 
different combinations possible, is twelve. Of these twelve, the options which were 
examined and presented in the next chapter are, total initial orthogonalizations with 
all three reorthogonalization options, and initial orthogonalization with the previous 
two blocks and with selective reorthogonalization, for a subtotal of four cases. 
Orthonormalization with respect to the mass matrix and the Euclidean norm for the 
above four cases was also performed for a total of eight Gram-Schmidt options 
considered. The computational time required for creating component mode models 
from various finite element models, using the various Gram-Schmidt options discussed 
above, is presented in chapter 5. 

4,2.5) Choleskv/OR Orthogonalization: 

Alternatives to Gram-Schmidt orthonormalization exist which perform the 
orthogonalization in a matrix format, rather than a vector by vector format. Use of 
explicit matrix orthogonalization can be an advantage in certain programming 
applications. Gram-Schmidt is classified as a method of performing a B = QR 
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decomposition, where B is a set of vectors (such as unorthogonalized Krylov vectors), 
Q is the orthonormalized vector set, and R is an upper triangular matrix, which can 
be assembled from the coefficients in the Gram-Schmidt algorithm. In typical Gram- 
Schmidt the R matrix is not assembled or used explicitly. The QR decomposition can 
be performed explicitly by using a Householder technique 18 or by recognizing that the 
R matrix is also the upper Cholesky factor of the B T B subspace 22 . It is instructive to 
note the similarity of the B T B matrix multiplication with equation (4.6), the 
differences being that all vectors, new and previous, are included in the B matrix. 
Cholesky/QR orthogonalization may be understood as a algorithm which, first, 
locating the non-orthogonal vectors by the B T B multiplication, and second, determines 
a transformation which will shift those vectors to an orthogonal space. 

That the R transformation matrix can be determined by Cholesky 
decomposition of the B T B subspace is demonstrated in the following equations. First, 
the Cholesky factor is defined as follows: if A is positive definite and symmetric then 
the LU decomposition, 

A =LDU ( 4 - 9 ) 


is equivalent to 


a=ldl t 


(4.10) 
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which may be written as 

A =LD V2 (LD V2 ) T = C T C ( 411 ) 

where C = D I/2 L T , and is called the Cholesky factor. That B T B is symmetric can be 
shown as follows 

B t B = {B t B) t = B t (B t ) t ( 4 - 12 ) 

Since the QR decomposition is defined by B = QR, then the B T B matrix 
product is 

B T B = (QR) t QR =r t q t qr ( 4 13) 

and since, by definition, Q T Q = /, then equation (3.24) becomes 

B T B=R T R = C T C (414) 

and because both R and C are upper triangular matrices then R = C, and therefore R 
is the upper Cholesky factor of B T B. 

The steps in the Cholesky/QR orthogonalization algorithm can be summarized 
as follows: 

1) The matrix product B T B , or using the notation of chapter 2, q h ■ T q h *, is 
calculated, where q t * is the concatenation of the vectors <7 l y ., defined in equation 
(2.11) and the vector block q } defined in equation (2.8), 



52 


2) The matrix product B T B, is decomposed into, R T R, where R is upper triangular. 

3) The equation, R T Q T = B, is solved, with only one forward substitution required 
since R T is lower triangular. 

The Cholesky/QR decomposition algorithm can also be used to orthonormalize 
vectors with respect to the mass matrix. If equation (4.13) is re-written as 

B T MB = (QR) T MQR = R T Q T MQR ( 415 ) 

and since Q T MQ = /, if orthonormalizing with respect to the mass matrix, then 
equation (4.14) can be re-written as 

b t mb = r t r = c t c ( 416 ) 

which would yield a different transformation matrix, R, and a different set of 
orthogonal vectors, Q, than the previous example of orthonormalizing with respect to 
the Euclidean vector norm. As discussed in section 2.4, the vectors must be 
orthonormalized with respect to the mass matrix for the tridiagonal form to be 
achieved, if the Krylov sequence of equation 2.2 is used. The practical aspect of the 
tridiagonal form in boundary flexibility component synthesis is that each new vector 
block theoretically only needs to be orthogonalized with the previous two vector 
blocks. 

It is evident, upon examination of the above algorithm, that Cholesky/QR 
decomposition orthonormalizes a set, or subset, of vectors simultaneously. This set, 
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or subset, would include all the vectors within a specific Krylov block, and these 
vectors are nonorthogonal. However, it is important to note that the previous vectors, 
which have been successfully orthonormalized, would produce a partition of B J B that 
is the identity matrix. Decomposing and solving an identity partition produces a 
transformation which does not alter vectors which formed the mutually orthonormal 
subset. As a result, the theoretical description of the tridiagonal form is applicable to 
vectors orthonormalized using Cholesky/QR decomposition, if, as presented above, 
they have been orthonormalized with respect to the mass matrix. Theoretically then, 
it is possible to orthonormalize, using Cholesky/QR, each new vector, or block, with 
the previous two vectors, or vector blocks, and be orthonormal to all previous vectors, 
in an identical manner as in the Gram-Schmidt algorithm. Practically, the 
Cholesky/QR algorithm has an advantage over Gram-Schmidt in that, when the loss 
of orthogonality occurs within the Lanczos algorithm, the inevitably of which is 
discussed was section 4.2.1, it does not need to be eliminated immediately. All non- 
orthogonal vectors can be reorthogonalized simultaneously at intervals, as required, 
and at the termination of the Krylov sequence. The numerous repetitions of the 
orthogonality checks required in some implementations of the Gram-Schmidt 
algorithm is not required in the Cholesky/QR algorithm. However, a method to 
automatically determine when a reorthogonalization is required, such as presented in 
reference [26] for Gram-Schmidt, would be required if orthogonalization with the mass 
matrix and the previous two blocks only were to be implemented practically. 
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4.3) Issues Concerning the Use of Blocks in the Boundary Flexibility Method: 

In the boundary flexibility method of component synthesis, several issues of 
orthogonalization and usage present themselves, due to the use of blocks in the Krylov 
sequence. Nour-Omid and Clough 20 described how vector blocks may be used with 
Wilson’s algorithm and that the results pertaining to the tridiagonal form and 
orthogonalization with the previous two blocks, described in section 2.3, are 
applicable. A difference resulting from the use of blocks is that orthogonalization 
must be performed with the previous two blocks, not merely the previous two vectors. 
Those equations will not be repeated here. The need for orthogonalization within each 
block was also discussed by Nour-Omid and Clough. However, it was not presented 
that, depending on the choice of the initializing vectors (forces), the vectors within the 
block may be approaching linear dependency. In the boundary flexibility method of 
component synthesis, nearly dependent vectors within a block can and do occur. This 
is not a numerical convergence of the Krylov sequence, as discussed in section 4.2.1. 
The nearly dependent vectors within a block is the natural result of the static solution 
of a structure under generalized loading, such as the mass matrix multiplied by the 
boundary flexibility matrix as in this dissertation. An illustrating example, is 
presented in section 4.3.1. 

As presented in chapter 2, the initializing block of forces in the boundary 
flexibility method is the mass matrix multiplied by the constraint modes, or the 
flexibility modes. The size of these matrices is therefore dependent on the size of the 
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of the boundary set of the structural component under consideration. In typical 
aerospace applications the boundary set of these models can be quite large, larger than 
can be practically used in the block-Krylov sequence. Block size reduction by 
filtering can correct an excessive block size and this addition to the boundary 
flexibility method of component synthesis is presented in section 4.3.2. 

4.3.1) Dependence of Vectors Within the Block: 

In the boundary flexibility method of component synthesis, nearly dependent 
vectors within a block result from the static solution of a structure under generalized 
loading. Chapter 2 details this generalized loading as the mass matrix multiplied by 
the boundary flexibility matrix. The initial set of Ritz vectors is the static response 
the component exhibits for the generalized loading. There is no theoretical basis to 
expect that this set of vectors within the initial, or any other, block should be linearly 
independent. The boundary flexibility algorithm, as presented in reference [1] and 
reviewed in chapter 2, made no orthogonality check of the vectors within the Krylov 
block. Furthermore, if internal block orthogonalization is not performed, subsequent 
Gram-Schmidt steps are ineffective because the blocks are orthogonalized with a set 
of vectors that are not orthogonal. 

Consider the case of a simple beam represented by a finite element model 
consisting of only two nodes and one element. This single element model is used as 
an example to demonstrate that the vectors within the initial Krylov block can be 
almost linearly dependent. The beam element will be processed as a component with 
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a fixed node interface. The interior partition of the component stiffness and mass 
matrix, for the simple beam element model, may be obtained directly from beam 
theory 7 . 
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The inverse of is as follows, 


M El 


l J_ 

3 21 

21 L 2 


(4.20) 


The constraint modes are given by equation (2.5), as 


d> = 

^ic "'ll K'i c 


-1 -L 
0 -1 


(4.21) 
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Because the component interface is statically determinate, the constraint modes are 
identical to the rigid body displacement matrix. The initial block of static Ritz 
vectors, q", as given by equation (2.6), is 


a ** = k 1 m $ = - 
9i "a m ii^ic 


ml? 

420 El 


41 35.6661 
* 49 


(4.22) 


assuming a lumped mass approach. The vectors contained in q" can be examined 
for dependency. For the first vector, q u over q 2I is equal to .732 1L. For the second 
vector, q I2 over q 22 is equal to .7279L. These two vectors are, once normalized, 
almost identical and linearly dependent. 

A large amount of dependency can also occur in larger blocks. Consider a 
structure with a statically indeterminate interface. In many cases the nodes may be 
positioned closely together, or in a symmetric fashion, either of which may result in 
some of the constraint modes being nearly identical. Since the initial set of Ritz 
vectors is the static displacement of the component to the mass matrix times the 
constraint modes, it is obvious that many component models will yield nearly 
dependent vectors, in the initial block. 

As mentioned previously, Nour-Omid and Clough 20 presented the requirement 
for orthogonalization within each block of vectors. They suggested that the Gram- 
Schmidt orthogonalization procedure, as shown in equations (2.9) through (2.12) be 
applied in a two step process. First, the Gram-Schmidt orthonormalizing is applied 
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to each vector within each Krylov block. Then, orthogonalization is performed with 
the previous vector blocks. An alternate solution to the problem of dependent vectors 
within the blocks, was presented by Su and Craig in reference [27]. In this solution 
a singular value decomposition is performed on the q’* T mq** subspace. The obtained 
transformation matrix orthogonalizes the block. 

The recommended method for orthonormalizing the vectors within the Krylov 
block has been discussed previously in section 4.2.5, Cholesky/QR orthogonalization. 
Some advantages of using Cholesky/QR decomposition to obtain an orthonormalizing 
transformation matrix have already been discussed. In relation to vector blocks, this 
methodology can orthogonalize a new block separately, or with the previous two 
blocks, or with all previous blocks, simultaneously. One potential disadvantage is 
that, if the vectors are almost linearly dependent, than the matrix to be decomposed 
is numerically singular. As discussed above, this near dependence is to be expected 
in even simple problems. The potential singularly of the B T B matrix product can be 
solved by block filtering, which will be discussed next. 

4.3.21 Block Filtering: 

Block filtering is a procedure by which nearly identical, or dependent, vectors 
are removed from a vector block. It is based upon a standard orthogonality check. 
The use of block filtering simultaneously solves two problems. First, by filtering the 
vectors, the size of the block created by the boundary flexibility method can be 
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reduced to a manageable level. Second, it can eliminate the potential singularities of 
the B T B matrix product which is used in Cholesky/QR orthonormalization. 

The first function of block filtering is block size reduction which corrects 
excessive block size. The size of the block in the boundary flexibility method is 
determined by the size of the boundary set of the structural component. The initial 
block of forces in the boundary flexibility method is the mass matrix multiplied by the 
constraint modes, or the flexibility modes. For the fixed interface method the equation 
establishing the first block, (2.6), is repeated here. 

«r = v l (" | «®fc +m fc) (26) 

The boundary set of practical structural models, which require the use of component 
mode synthesis, is usually large enough to cause problems in use of the algorithm. 
For instance, a typical Space Shuttle cargo element component model might have a 
boundary set of forty-eight degrees of freedom, eight nodes with six degrees of 
freedom each. This would lead to a block size of forty-eight. In contrast, the default 
block size in the MSC/NASTRAN implementation of the Lanczos eigenvalue 
extraction method 32 is seven, with a maximum of fifteen. 

There are several reasons why a large block size is a disadvantage. First, as 
discussed in the previous section, many of the vectors in the blocks may be nearly 
identical. The information retained after orthogonalizing these nearly dependent 
vectors may not be significant, and in extreme cases may only be the product of 
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numerical roundoff. In these cases the Ritz vectors obtained from the second Krylov 
sequence iteration are more likely to be dynamically significant than the product of 
numerical roundoff. Second, truncation of the Krylov sequence and final resulting set 
of vectors becomes very imprecise. The static Ritz vectors are generated block by 
block and so if the block size is very large many more vectors than desired may be 
obtained. Finally, with large blocks, orthogonalization within each block becomes 
more expensive and can become numerically difficult due to the previously discussed 
dependencies. 

The second function of block filtering is the elimination the singularities in the 
B T B matrix product, which is used in Cholesky/QR orthonormalization. If two 
vectors in the B matrix are nearly identical than, after normalization, the B T B matrix 
will have a unity term on both (lower and upper) off-diagonal positions corresponding 
to the column number of the identical vectors. Since each row of the B r B matrix has 
a unity term in the diagonal position, the two rows corresponding to the two identical 
vectors will be dependent, and Cholesky decomposition becomes problematic. If one 
of the two vectors is eliminated by block filtering than decomposition can be 
accomplished. No information is lost in the block filtering because, by definition, the 
vectors in question are nearly identical, and so the discarded vector is a duplication. 

The block filtering procedure can be summarized as follows: 

1) The vectors within the block, q*, are normalized as follows 

q]' = q* [<diag{q/ T m u q/)>' c l 12 ] 


( 4 . 23 ) 
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or, depending on the orthonormalization selection, 

*;■-»/ [<*«*(«/ '*/)»«"] 


(4.24) 


where rfiag is defined as the diagonal terms of the matrix product and the exponent, 
' m , is applied to each term within the resulting diagonal matrix, not to the complete 
matrix itself. 

2) The cross-orthogonality of the vectors is calculated 


L = a* T m.Q . 

m “j u “ ] 


(4.25) 


or, depending on the orthonormalization selection, 


l e 



(4.26) 


At this point no orthogonalization has occurred and so there is no reason to expect the 
matrix product to result in the identity matrix. In this way, the above equations differ 
from equation (4.6) and (4.7). 

3) The L matrix is partitioned into its lower triangular portion, excluding the diagonal 
terms. 

4) The infinity norm, as defined in equation (4.8), of each column in the resulting 
lower triangular matrix is calculated. These terms are compared to an arbitrary filter 
value, e, determined by practice to be initially set to .995. Vectors with associated 
terms greater than this value are partitioned from the normalized block. 
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5) The size of the revised block is determined and if the block is too small or too 
large then the filter value is raised or lowered, respectively, and step 4) is 
subsequently repeated. The minimum value for a block has been set at six and the 
maximum at eighteen. 

4.4) Summary and Revised Orthogonalization Algorithm: 

The following tables contain revised orthogonalization algorithms for the 
boundary flexibility method of component synthesis, using static Ritz vectors. These 
algorithms are a synthesis of the basic methodology described in section 2.2 and the 
revisions and additions to the method which have been presented in this chapter. 
Table (4.1) presents the algorithm using the mass matrix for orthonormalization. This 
algorithm initially orthonormalizes with the previous two vector blocks, and then at 
intervals, and at the termination of the sequence, full reorthogonalization occurs. 
Table (4.2) presents the algorithm using the Euclidean vector norm for 
orthonormalization. Both algorithms are presented for fixed interface components, 
however, for free interface components the body of algorithms presented are identical. 
The initialization and transformation of the free interface component are different, as 


documented in section 2.2. 
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After the assembly of the component mass, m, and stiffness, k, matrices: 

1) Initialization 

® = — jt -1 jfc [create boundary flexibility matrix, k u 

ic ii ic 

and k ic were defined in equation (2.3)] 
q* - fc -i ( m $ +m . ) [create the initial block] 


(For a free interface component, (If the component has rigid body modes then 
the elastic flexibility matrix, g t , defined by (2.35) through (2.39) is used): 


8 a 


8cc 

Sic 


where g = k~ l 


[boundary flexibility matrix] 


<h =gmg a 


[create the initial block]) 


2) Filtering of Initial Block 

*r = *•[<**&/ 7 *««,•)> j’ va 

^.=*r r *«*r 

\ O' 

L_ - - 

\J 

If (< |£j„> c ) * dependent 


8i 


Qdep 


[normalization] 

[cross orthogonality] 

[partition into lower triangular] 

[infinity norm of each vector] 
[partition out dependent vectors] 


Table (4.1) - Revised Boundary Flexibility Algorithm 
Using Orthonormalization With Respect to the Mass Matrix 
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3) Oithonormalization of the Initial Block 


L * _ ~ *** T *** 

m~9 1 m ii9 1 


[B r MB matrix product] 


L' m -a?R 


[Cholesky factor decomposition] 


j x +** 

R T q x = q x 


[Solve by forward substitution] 


4) For Blocks j = 2,3, .... / 
9j - V lm »«/-i 


q]' * r m »9/) > ccl 1/2 

\ 0 

L - - 

[ L * N 

If (<|L m L> ) * e, fftan dependent 


[Krylov sequence] 
[normalization] 

[cross orthogonality] 

[partition into lower triangular] 
[infinity norm of each vector] 


* * [" * * * 
9j "* [ 9 / 


[partition out dependent vectors] 


5) Orthonormalization (At intervals, orthonormalize with all previous blocks) 

(when j = /, transform system to form component) 


9j.j-l.j-2 m H 9j.j-l.j-2 


[B t MB matrix product] 


L‘ m =R T R 


[Cholesky factor decomposition] 


R t 


T *** 

9j.j-l.j-2 ~ 9j.j-l.j-2 


[Solve by forward substitution] 


Table (4.1) (Continued) - Revised Boundary Flexibility Algorithm 
Using Orthonormalization With Respect to the Mass Matrix 
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After the assembly of the component mass, m, and stiffness, k, matrices: 

1) Initialization 

(h = -ir -t k [create boundary flexibility matrix, k u 

ic ii ic 

and k ic were defined in equation (2.3)] 
q* = + w ic ) [create the initial block] 


[For a free interface component, (If the component has rigid body modes then 
the elastic flexibility matrix, g e , defined by (2.35) through (2.39) is used): 


S a 


8cc 

Sic 


where g=k~ x 


[boundary flexibility matrix] 


9i=gmg a 

2) Filtering of Initial Block 
q* t * =q l *[<diag(q 1 * T q l *) > C cY V2 

L E =qr T *r 

\ O' 

l f ~ - 

L E N 

If (< ||Z £ L> c ) ^ e, than dependent 

♦ * r * * * 

9i [qi Sdep 


[create the initial block]] 

[normalization] 

[cross orthogonality] 

[partition into lower triangular] 

[infinity norm of each vector] 
[partition out dependent vectors] 


Table (4.2) - Revised Boundary Flexibility Algorithm 
Using Orthonormalization With the Euclidean Vector Norm 
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3) Orthonormalization of the Initial Block 

L * * * * T — *** 

£ = 9i 

l e =r t r 

j T *** 

R 9i = 9i 


[B t B matrix product] 

[Cholesky factor decomposition] 
[Solve by forward substitution] 


4) For Blocks j = 2,3, — , l 

9 

t" ‘q;i«Ba g (t; T q/->>„r m 

L _ * * 7* — * * 

E-ij % 

\ o 

l e \ 


[Krylov sequence] 


[normalization] 


[cross orthogonality] 


[partition into lower triangular] 


If (< || L L> ) * e, than dependent [infinity norm of each vector] 


Qj [ 9j Qdep 


[partition out dependent vectors] 


5) Orthonormalization (each block) 

(when j = /, transform system to form component) 


*** T „ *** 


L E=<h j l 9u 


l e =r t r 


nT T *** 

R T 9ij = 9ij 


[B t B matrix product] 

[Cholesky factor decomposition] 
[Solve by forward substitution] 


Table (4.2) (Continued) - Revised Boundary Flexibility Algorithm 
Using Orthonormalization With Respect the Euclidean Vector Norm 
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Figure (4.1) - Normalized Vectors Obtained from a Krylov Sequence of 
Order Nine, with x 2 as the Iteration Matrix 


Chapter 5 

Numerical Examples of Orthonormalization 

5.11 Introduction: 

The following chapter describes the method with which the algorithms 
presented in chapter 4 were implemented. Also included are a number of illustrative 
examples of theoretical concepts, timing comparisons, and descriptions of algorithm 
trials on practical models. The software and hardware, and finite element models used 
in these examples is also described. 

5.2) Tools and Programming: 

The algorithms presented in section 4.3 were implemented in 
MSC/NASTRAN 32 , Version 67. The use of a standard, commercially available 
computer program allows the results of this work to be transferred easily to other 
structural dynamists. Adding the static Ritz vector algorithms to MSC/NASTRAN is 
implemented by the use of the internal programming language called DMAP (Direct 
Matrix Abstraction Programming). The standard solution sequences of 
MSC/NASTRAN are written in DMAP, and the source code of MSC/NASTRAN is 
available at the DMAP level. For example, equations (2.3) through (2.5), (creation 
of structural mass and stiffness matrices, partition, and constraint mode creation) and 
(2.16) through (2.23) (transformation into the modal component) are currently 
contained in the standard MSC/NASTRAN normal modes solution sequence. The 
equations presented in Tables (4.1) and (4.2) were written using DMAP, and were 
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then incorporated into the MSC/NASTRAN normal modes solution sequence. In 
addition, Gram-Schmidt orthonormalization was also implemented, to allow 
comparisons with the presented Cholesky/QR decomposition algorithm. 

A Cray X-MP running under the UNICOS operating system was the resident 
hardware running the version of MSC/NASTRAN utilized. This fact is relevant when 
considering the timing numbers presented in section 5.4. 

5.31 Finite Element Models: 

All finite elements models were created using standard MSC/NASTRAN. 
They required no special processing. The normal modes solution sequence created the 
finite element component mass and stiffness matrices and performed the partitions into 
internal and external degrees of freedom. They were then ready to be processed by 
the boundary flexibility algorithm using Ritz vectors. 

5.3.1) Simple Beam Model: 

The first example case of a simple beam was derived from a finite element 
model of the Space Station Freedom photovoltaic array central mast. The length of 
the beam was 1 179.9 inches. The modulus of elasticity, E, was 10.1 x 10 6 lbs/in 2 and 
the moment of inertia of the cross section was 108.9 in 4 . Its weight per unit length 
was .2296 lbs/in. The simple beam was modeled with eleven nodes and ten beam 
finite elements. Several different boundary conditions, both at the component and 
system level, were imposed upon this beam, yielding cantilevered and free-free 
conditions. The different boundary conditions cases will be described in section 5.4. 
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5.3.2) Space Station Electrical Power System Radiator: 

The next example used was a finite element model of the Space Station 
Electrical Power System (EPS) Radiator (Figure (5.1)). The main contractor 
constructing the EPS Radiator is Loral Vought Systems. Its purpose is to expel excess 
heat created by Space Station Electrical Power System into space. The structure 
weighs approximately 1440 pounds and, when deployed, is approximately 50 feet long 
and 1 2 feet wide. The finite element model representing this structure was created by 
Loral Vought Systems (Figure (5.2)). This finite element model contains 
approximately 4000 degrees of freedom. With the boundary degrees of freedom fixed, 
the EPS Radiator finite element produces eight normal modes below 5 Hz, the first 
three being at .19, .73, and .94 Hz (Figures (5.3), (5.4), and (5.5), respectively). 

5.3.3) Cassini Spacecraft: 

The third example used was a finite element model of the Cassini Spacecraft 
deep space probe (Figure (5.6)). The primary organization responsible for the Cassini 
spacecraft is the Jet Propulsion Laboratory of the California Institute of Technology 
and NASA. This spacecraft will be launched upon a Titan IV launch vehicle and will 
explore the Saturn planetary system. The structure weighs approximately 12,890 
pounds and is approximately 23 feet long and 14 feet wide. The finite element model 
representing this structure was created by the Jet Propulsion Laboratory (Figure (5.7)). 
This finite element model contains approximately 11,100 degrees of freedom. With 
the boundary degrees of freedom fixed, this Cassini finite element model produces 
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sixty-three normal modes below 50 Hz, the first three primary structural modes being 
at 7.36 (bending), 7.70 (bending), and 15.78 Hz (torsion). 

5.4) Numerical Results: 

The following numerical examples can be divided into two groups. The first, 
contained in sections 5.4.1 and 5.4.2, are illustrations of some of the theoretical 
properties, presented in chapter 4, of generalized static Ritz vectors in the boundary 
flexibility method. These examples also serve to validate the correctness of the 
implementation of the presented theory. Sections 5.4.3 and 5.4.4 contain examples 
of a more practical nature. These examples serve to provide a physical understanding 
of the methodology and its benefits in terms of computer cost. 

5.4.1) Numerical Illustrations of Theoretical Properties: 

If a number of Ritz vectors, equal to the number of degrees of freedom in a 
finite element model, are used to form a component, then those vectors do not 
represent a Ritz approximation but are an exact transformation. That a component 
so formed is exact was proven in section 3.2. A demonstration of that proof, and of 
the correct implementation of the boundary flexibility/static Ritz vector methodology 
and algorithms presented in chapter 4, is shown in Table (5.1). A complete set of 
Ritz vectors, equal to the number of degrees of freedom in the beam model, were 
calculated. The interface was assumed to be fixed and so the beam was cantilevered. 
A boundary flexibility component model was subsequently formed. The normal 
eigenvalues of the transformed component were then calculated and compared to the 
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eigenvalues obtained directly from the finite element model of the cantilevered beam. 
This is the comparison contained in Table (5.1). 

The cross orthogonality of the Ritz vectors, as defined in equations (3.6) and 
(3.7), used in the exact transformation discussed in the preceding paragraph, were 
calculated. The six largest pairs of off-diagonal terms resulting from the cross 
orthogonality calculation for two different cases are given in Tables (5.2) and (5.3). 
The vectors used to form the matrix, from which the terms contained in Table (5.2) 
were extracted, were orthonormalized with respect to mass matrix. In order to 
facilitate obtaining the exact transformation, orthonormalization was performed on all 
previous vectors at each step in the Krylov sequence, not with just the previous two 
vector blocks. The vectors used to form the matrix, from which the terms contained 
in Table (5.3) were extracted, were orthonormalized using the Euclidean vector norm. 

The orthogonality properties of the Krylov sequence are illustrated in Tables 
(5.4) and (5.5). As was discussed in section 2.3, theoretically, if Ritz vectors are 
orthonormalized with respect to the mass matrix and the preceding two blocks of 
vectors, then they are orthogonal with all previously calculated vectors. Table (5.4) 
contains the orthogonality check of the first five blocks of vectors produced for the 
cantilevered beam example described above, orthonormalized with respect to the mass 
matrix. Theoretically, all vectors should be mutually orthonormal and this matrix 
should be the identity matrix. The extremely large terms at position (1,30) and (30,1) 
of the matrix are illustrative of the inherent breakdown of orthogonality which was 
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discussed in section 4.2.1, and demonstrated by Theorem 4.1. The rapid creation of 
non-orthogonal vectors in this small idealized example clearly demonstrates the 
requirement for re-orthogonalization. When calculating Ritz vectors for a more 
realistic problem, the requirement for continual re-orthogonalization can outweigh the 
advantages of orthogonalization with just the previous two blocks and it can become 
more efficient to orthogonalize with all previous vector blocks. 

Table (5.5) also contains an orthogonality check with the previous two vector 
blocks. However, in this example the vectors have been orthogonalized using the 
Euclidean vector norm. There is no theoretical reason why these vectors should be 
orthogonal with all previous vectors and they are not. The vectors are orthogonal with 
the previous two blocks, but non-orthogonality is manifest between the other vectors. 
Specifically, the vectors in block four of the example are not orthogonal with those 
in block one and the vectors in block five are not orthogonal with the vectors in 
blocks one and two. As a result, whatever advantages exist for using 
orthonormalization with the Euclidean vector norm must be weighed against the 
requirement for orthogonalization with all previous vectors or more difficult re- 
orthogonalization. 

5.4.2) Numerical Illustrations of Block Issues: 

The Krylov block issues of dependence, size, and filtering, which were 
discussed in section 4.3, can be illustrated by the following example, using the EPS 
Radiator finite element model. This model has a boundary set of six nodes, which is 
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a typical size in practical aerospace applications. The six nodes each have six degrees 
of freedom and so the boundary set in this component model has thirty-six degrees of 
freedom, which in turn leads to thirty-six constraint modes and an initial block size 
of thirty-six. Clearly, as discussed in section 4.3, this is an unwieldy amount of 
vectors to process. The cross-orthogonality of the initial block is presented in Table 
(5.6). As can be seen, this is a large matrix which contains many nearly identical 
vectors, as identified by the many cross-orthogonality terms approaching unity. This 
is a further demonstration of the initial nonorthogonality of vectors within a block. 
The initial, large block of vectors was filtered to produce a new set of eight vectors. 
The cross-orthogonality of the filtered vector block presented in Table (5.7). Since 
only nearly identical vectors were removed, the filtered vector block contains almost 
the identical response information as the original block, is a more convenient size to 
work with, and is numerically cleaner and easier to orthonormalize. 

5.4.3) Simple Beam Numerical Results: 

A component representation of the ten element beam finite element model was 
created using the boundary flexibility method and static Ritz vectors. The fixed 
interface approach, with constraint modes and two Krylov blocks, was used to form 
the component. The interface of the component consisted of one node and six degrees 
of freedom. The number of constraint nodes is equal to the number of interface 
degrees of freedom, and the size of the Krylov block is equal to the number of 
constraint modes. Therefore, each Krylov block contained six Ritz vectors. Since the 
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component model was formed with two Krylov blocks, it contained a total of twelve 
generalized coordinates. 

Plots of the lateral static Ritz vectors, which represented the cantilevered beam, 
are shown in Figures (5.8-11). The unorthogonalized vectors, as produced by 
equations (2.6) and (2.8), or as given in Table (4.1), are shown in Figures (5.8-9)(a) 
and (5. 10-1 l)(a), respectively. The first normalized vector, as produced by equation 
(2.7), is given in Figure (5.8)(b). (The first vector does not need to be 
orthogonalized.) The remaining orthogonalized and normalized vectors, as produced 
by equation (2.12), are given in Figures (5.9-1 l)(b). The first two unorthogonalized 
vectors (given in Figures (5.8-9)(a)), which are in the first Krylov block, appear to be 
nearly identical. The first static Ritz vector produced by the boundary flexibility 
method is similar to the classic first bending normal mode shape of a beam. After 
orthonormalization, the second static Ritz vector has been modified into the classic 
second bending normal mode shape (as shown in Figures (5.8-9)(b)). 

The similarity of the Ritz vectors to the normal modes of a cantilevered beam 
provides an important insight into the numerical difficulties of orthogonalizing vectors 
which result from the Krylov sequence. It was shown in theorem 4.1, section 4.2.1, 
that the Krylov sequence converges to an eigenvector. In other words, when an 
eigenvector is input into the sequence, it produces the same eigenvector. The resulting 
complete vector dependency will cause numerical orthogonalization to be very 
difficult. Even when the Ritz vectors closely resemble the normal modes, but are not 



76 


exactly the normal modes, near duplication of these vectors from the subsequent 
Krylov iteration occurs. This convergence to normal modes is the potentially 
important disadvantage of creating boundary flexibility static Ritz vectors which 
closely resemble normal modes. This convergence also partially explains the high 
degree of nonorthogonality evident in Tables (5.4) and (5.5). 

The eigenvalues of the reduced component subspace were then calculated to 
evaluate the static Ritz vector representation. The first five single plane natural 
frequencies from this reduced system are shown in Table (5.8). For comparison, 
Table (5.8) also includes the frequencies of a reduced system where the component 
was formed using traditional normal modes. This component was also formed with 
a fixed interface, but thirteen normal modes were used for numerical convenience. 
The full, or "exact", finite element eigenvalue solution is also shown. In the case of 
the Ritz vectors, no modal selection of any kind was used. For the case of the normal 
modal component, modal selection by truncation was used. The superior accuracy of 
the normal modal component, in the fourth bending mode, does not necessarily 
represent a limitation of the boundary flexibility methodology, but instead 
demonstrates the requirement for an adequate Krylov sequence truncation criteria. 
This subject will be discussed in chapter 7. 

In addition to the fixed interface example, two free interface examples were 
created. Both were based on the same ten element beam finite element model, but 
with different boundary conditions. The first free component model considered was 
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free-fixed, with the component interface being the free boundary condition, and hence 
it had no rigid-body modes. Equations (2.20) through (2.34) define the formulation 
of this free interface component. This component was also formed with two Krylov 
blocks, and therefore twelve Ritz vectors. Table (5.9) presents the first five single 
plane system natural frequencies of the reduced component subspace, compared to the 
natural frequencies of the full finite element model. The second free component 
model considered consisted of the ten element beam with free-free boundary 
conditions and six rigid body modes. Equations (2.35) through (2.42) define the 
formulation of this component. Table (5.10) contains the six rigid body frequencies 
and the first four single plane elastic frequencies of the reduced component subspace, 
compared to the frequencies of the full finite element model. In the two free interface 
cases, there is no comparison with the normal component mode synthesis. This is 
because standard free-interface MSC/NASTRAN routine does not use the ' Rubin- 
MacNeal" method, and so a direct comparison was not performed. 

5.4.4) Timin 2 Comparisons: 

Static Ritz vector component models were formed using the boundary 
flexibility methodology for the EPS Radiator finite element model and the Cassini 
spacecraft finite element model. A variety of options was used in performing these 
computer runs. Orthonormalization was performed using Cholesky/QR and Gram- 
Schmidt, and with respect to the mass matrix and with the Euclidean vector norm. 
When orthonormalizing with respect to the mass matrix, orthonormalization was 
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performed both with all previous blocks and with the previous two blocks, as allowed 
theoretically. When Gram-Schmidt was used, reorthogonalization was performed with 
all vectors, with selected vectors, and with selected vectors using modified Gram- 
Schmidt. The timing comparisons are presented in Table (5.1 1) for the EPS Radiator, 
and Table (5.12) for the Cassini. In these comparisons a specified, consistent number 
of Ritz vectors was created. The number of vectors specified was sufficient to create 
an accurate component, for an arbitrary frequency cutoff. A Cray X-MP, using the 
UNICOS operating system, was the computer hardware system used to perform these 
timing comparisons. 

The algorithms were implemented using DMAP, in version 67 of 
MSC/NASTRAN. As a result, some caution should be used in interpreting the timing 
comparison data. Each DMAP module calls an independent set of compiled fortran 
routines, and each call takes a certain amount of computer time, which is essentially 
overhead. As a result, a large sequence of DMAP statements, especially a loop which 
will be repeated many times, will not be efficient as programs written in some 
compiled computer languages, such as FORTRAN. The results of the modified Gram- 
Schmidt reorthogonalization option would be particularly misleading, because of the 
large number of separate DMAP calls. On the positive side, the results for static Ritz 
vector, boundary flexibility component synthesis represent a lower bound estimate of 
the likely improvement in computer time. If implemented more efficiently, such as 
in the FORTRAN code of NASTRAN, this methodology should demonstrate a greater 
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time savings in forming the component then the already significant amount shown for 
the Cassini model in Table (5.12). The Cassini spacecraft model was the largest finite 
element model considered. 
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Mode 

Ritz 

Normal Modes 

1 

.5463995 

.5463995 

2 

.5463995 

.5463995 

3 

3.414649 

3.414649 

4 

3.414649 

3.414649 

5 

9.522067 

9.522067 

6 

9.522067 

9.522067 

7 

18.56508 

18.56508 

8 

18.56508 

18.56508 

9 

27.57835 

27.57835 

10 

30.53347 

30.53347 

11 

30.53347 

30.53347 

12 

45.41006 

45.41006 

13 

45.41006 

45.41006 

14 

63.18172 

63.18172 

15 

63.18172 

63.18172 

16 

81.89710 

81.89710 

17 

83.59296 

83.59296 

18 

83.59296 

83.59296 

19 

84.63568 

84.63568 

20 

104.4086 

104.4086 

21 

104.4086 

104.4086 

22 

133.7274 

133.7274 

23 

151.3342 

151.3342 

24 

151.3342 

151.3342 

25 

181.4946 

181.4946 

26 

182.6918 

182.6918 

27 

182.6918 

182.6918 

28 

223.0661 

223.0661 

29 

223.0661 

223.0661 

30 

223.7470 

223.7470 

31 

251.5984 

251.5984 

32 

259.2011 

259.201 1 

33 

271.7271 

271.7271 

34 

271.7271 

271.7271 

35 

286.7794 

286.7794 


Table (5.1) - A Comparison of Cantilevered Beam Frequencies (Hz) 
Obtained From a Static Ritz Component Transformed into Normal Eigenvalues 
and a Direct Normal Eigenvalue Solution 
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Mode 

Ritz 

Normal Modes 

36 

305.6441 

305.6441 

37 

315.2220 

315.2220 

38 

329.7375 

329.7375 

39 

329.7375 

329.7375 

40 

397.8221 

397.8221 

41 

397.8221 

397.8221 

42 

411.6982 

411.6982 

43 

473.4397 

473.4397 

44 

473.4397 

473.4397 

45 

543.8633 

543.8633 

46 

543.8633 

543.8633 

47 

560.5679 

560.5679 

48 

692.9518 

692.9518 

49 

692.9518 

692.9518 

50 

694.1468 

694.1468 

51 

808.7912 

808.7912 

52 

901.3738 

901.3738 

53 

969.3694 

969.3694 

54 

1010.923 

1010.923 


Table (5.1)(Continued) - A Comparison of Cantilevered Beam Frequencies (Hz) 
Obtained From a Static Ritz Component Transformed into Normal Eigenvalues 
and a Direct Normal Eigenvalue Solution 
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Row 1, Selected Columns 

1.1) 1.0000e+00 1,54) 1.9043e-14 

Row 6. Selected Columns 

6,6) 1.0000e+00 6,54) -1.4917e-14 

Row 12, Selected Columns 

12.12) 1.0000e+00 12,54) -1.6464e-14 

Row 13. Selected Columns 

13.13) 1.0000e+00 13,54) 1.3798e-14 

Row 39. Selected Columns 

39.39) 1.0000e+00 39,51) -1.3741e-14 

Row 49. Selected Columns 

49.49) 1.0000e+00 49,54) -2.2577e- 14 

Row 51. Selected Columns 

51.39) -1.3741e-14 51,51) 1.0000e+00 

Row 54. Selected Columns 

54.1) 1.9043e-14 54,6) -1.4917e-14 54,12) -1.6464e-14 54,13) 1.3798e-14 

54.49) -2.2577e-14 54,54) 1.0000e+00 


Table (5.2) - The Largest Off-Diagonal Terms From an Orthogonality 
Check of the Static Ritz Vectors Representing a Cantilevered Beam, 
Using Orthonormalization With Respect to the Mass Matrix 
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Row 39, Selected Columns 

39.39) 1.0000e+00 39,51) 9.4066e-15 

Row 40, Selected Columns 

40.40) 1.0000e+00 40,51) 8.4148e-15 

Row 42, Selected Columns 

42.42) 1.0000e+00 42,54) -1.7574e-14 

Row 43, Selected Columns 

43.43) 1.0000e+00 43,54) 1.1849e-14 

Row 48, Selected Columns 

48.48) 1.0000e+00 48,54) 1.2303e-14 

Row 49, Selected Columns 

49.49) 1.0000e+00 49,54) 2.1308e-14 

Row 51, Selected Columns 


51,39) 9.4066e-15 

51,40) 

OO 

4^ 

E 

00 

0 

1 

U\ 

51,51) 

1.0000e+00 

Row 54, Selected Columns 




54,42) -1.7574e-14 

54,43) 

1.1849e-14 

54,48) 

1.2303e-14 

54,49) 2.1308e-14 

54,54) 

1.0000e+00 




Table (5.3) - The Largest Off-Diagonal Terms From an Orthogonality 
Check of the Static Ritz Vectors Representing a Cantilevered Beam, 
Using Orthonormalization With the Euclidean Vector Norm 
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Block One 

U) 

1.0000e+00 

2,2) 

1 .0000e+00 

3,3) 

1.0000e+00 

4,4) 

1.0000e+00 

5,5) 

1.0000e+00 

6,6) 

1.0000e+00 

Block Two 

7,7) 

1.0000e+00 

8,8) 

1.0000e+00 

9,9) 

1.0000e+00 

10,10) 

1.0000e+00 

11,11) 

1.0000e+00 

12,12) 

1.0000e+00 

Block Three 

13,13) 

1 .0000e+00 

14,14) 

1.0000e+00 

15,15) 

1.0000e+00 

16,16) 

1.0000e+00 

17,17) 

1.0000e+00 

18,18) 

1.0000e+00 

Block Four 

19,19) 

1.0000e+00 

20,20) 

1.0000e+00 

21,21) 

1.0000e+00 

22,22) 

1.0000e+00 

23,23) 

1.0000e+00 

24,24) 

1.0000e+00 

Block Five 

25,1) 

1.6720e-02 

26,26) 

1.0000e+00 

27,27) 

1 .0000e+00 

28,3) 

9.8957e-02 

29,29) 

1.0000e+00 

30,1) 

5.8706e-01 


1,25) 1.6720e-02 
3,28) 9.8957e-02 


25,25) 1.0000e+00 

28,28) 1.0000e+00 
30,30) 1.0000e+00 


1,30) 5.8706e-01 


Table (5.4) - The Cross Orthogonality of Static Ritz Vectors 
Created by Orthonormalizing With the Previous Two Blocks, 
Using Orthonormalization With Respect to the Mass Matrix 
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Block One 


U) 

l.OOOOe+OO 

1,19) -9.9924e-01 

2,2) 

l.OOOOe+OO 

2,20) -9.9135e-01 

3,3) 

l.OOOOe+OO 

3,21) -9.9924e-01 

4,4) 

l.OOOOe+OO 

4,21) -3.9085e-02 

5,5) 

l.OOOOe+OO 


6,6) 

l.OOOOe+OO 

6,19) 3.9085e-02 

Block Two 


7,7) 

l.OOOOe+OO 

7,25) -4.9289e-02 

8,8) 

l.OOOOe+OO 

8,26) -2.5613e-01 

9,9) 

l.OOOOe+OO 

9,27) -4.9289e-02 

10,10) 

l.OOOOe+OO 

10,28) -7.4190e-02 

1UD 

l.OOOOe+OO 


12,12) 

l.OOOOe+OO 

12,30) -7.4190e-02 

Block Three 


13,13) 

l.OOOOe+OO 


14,14) 

l.OOOOe+OO 


15,15) 

l.OOOOe+OO 


16,16) 

l.OOOOe+OO 


17,17) 

l.OOOOe+OO 


18,18) 

l.OOOOe+OO 


Block Four 


19,1) 

-9.9924e-01 

19,6) 3.9085e-02 

20,2) 

-9.9135e-01 

20,20) l.OOOOe+OO 

21,3) 

-9.9924e-01 

21,4) -3.9085e-02 

22,3) 

3.8427e-02 

22,4) -9.8199e-01 

23,23) 

l.OOOOe+OO 


24,1) 

-3.8427e-02 

24,6) -9.8199e-01 

Block Five 


25,6) 

-1.8463e-01 

25,7) -4.9289e-02 

26,2) 

-1.2688e-01 

26,8) -2.5613e-01 

27,4) 

1.8463e-01 

27,9) -4.9289e-02 

28,9) 

-9.960 le-01 

28,10) -7.4190e-02 

29,29) 

l.OOOOe+OO 


30,7) 

9.9601e-01 

30,12) -7.4190e-02 


1,24) 

2,26) 

3.22) 

4.22) 

-3.8427e-02 

-1.2688e-01 

3.8427e-02 

-9.8199e-01 

4,27) 

1.8463e-01 

6,24) 

-9.8199e-01 

6,25) 

-1.8463e-01 

7,30) 

9.960 le-01 



9,28) 

-9.960 le-01 




19,19) 1.0000e+00 

21.21) 1.0000e+00 

22.22) 1.0000e+00 

24.24) 1.0000e+00 

25.25) 1.0000e+00 

26.26) l.OOOOe+OO 

27.27) l.OOOOe+OO 

28.28) l.OOOOe+OO 

30,30) l.OOOOe+OO 


Table (5.5)- The Cross Orthogonality of Static Ritz Vectors 
Created by Orthonormalizing With the Previous Two Blocks, 
Using Orthonormalization With the Euclidean Vector Norm 
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Row h Columns 1 Thru 36 

I) 1.0000e+00 -9.9999e-01 -9.9604e-01 -9.9996e-01 9.5350e-01 

II) 7.5250e-01 -1.6221e-Ol -9.994le-01 -9.9906e-01 -9.9289e-01 

21) -9.9854e-01 -9.9967e-01 9.9860e-01 9.9959e-01 -9.7430e-0I 
30 9.995Ie-01 9.9954e-01 -9.9908e-01 -9.9963e-01 -9.9900e-0l 

Row 2. Columns 1 Thru 36 

0 -9.9999e-01 L.OOOOe+OO 9.960le-01 9.9995e-01 -9.5348e-01 

10 -7.5179e-01 1.5939e-01 9.993le-0I 9.9896e-01 9.9282e-OI 

21) 9.9867e-01 9.9969e-01 -9.9872e-01 -9.9964e-01 9.7493e-01 

30 -9.9937e-0] -9.9941e-01 9.9888e-0l 9.995Ie-01 9.9879e-01 
Row 3. Columns I Thru 36 

1) -9.9604e-01 9.9601e-01 1.0000e+00 9.9615e-0I -9.7650e-01 

0) -6.9905e-01 1.0225e-01 9.9328e-01 9.9217e-01 9.8160e-01 

21) 9.9044o *01 9.9359e-Oi -9.9057e-01 -9.9359e-01 9.84l8e-OI 

31) -9.9442e-01 -9.9452e-01 9.9475e-01 9.9450e-01 9.9456e-01 
Row 4. Columns 1 Thru 36 

0 -9.9996e-01 9.9995e-01 9.9615e-01 l.OOOOe+OO -9.5398e-01 

11) -7.5412e-01 1.6492e-01 9.9947e-01 9.9915e-01 9.9343c-01 

21) 9.9843e*01 9.9965e-0t -9.9848e-01 -9.996ie-01 9.7366e-01 

31) -9.9956e-01 -9.9959e-01 9.9901e-01 9.9963e-01 9.9892e-0I 
Row 5. Columns 1 Thru 36 

1) 9.5350e-01 *9.5348e*01 -9.7650e-01 *9.5398e*OI l.OOOOe+OO 

11) 5.4647e-01 5.0223e-02 -9.4547e-01 -9.4255e-01 -9.2212e-01 

21) -9.3835e-01 -9.4610e-01 9.3865e-01 9.4636e-01 -9.7639e-01 
31) 9.4900e-01 9.4928e-01 -9.5089e-01 -9.4899e-01 -9.5045e-01 

Row 6. Columns 1 Thru 36 

1) -9.9999e-01 l.OOOOe+OO 9.9589e-OI 9.9994e-01 -9.5308e-01 

11) -7.5237e-01 1.5987e-0l 9.9933e-01 9.9898e-01 9.9290e*01 

21) 9.9874e-01 9.9972e*01 -9.9879e-0l *9.9967e*01 9.748 le-01 

30 -9,9937e-01 -9.9940e-01 9.9886e-01 9.9951e-01 9.9878e-01 

Rpw 7„ l , T !u v ,36 

I) -9.7551e-Ol 9.7619e-01 9.8395e-01 9.7488e*0I -9.7285e-0I 

II) -6.071 le-01 -5.8122e-02 9.6749e-0i 9.6537e-OI 9.4835e-01 

21) 9.7127e-01 9.7230e*01 -9.7154e-01 -9.727 le-01 9.9978e-01 

31) -9.6844e-01 -9.6868e-01 9.6838e-OI 9.6942e-01 9.6805e-0l 
Row 8. Columns 1 Thru 36 

I) -9.9976e-01 9.9971e-01 9.9409e-01 9.9975e-01 -9.4745e-01 

II) -7.6648e-01 1.8253e-01 9.9989e-01 9.9975e-01 9.9490e-01 

21) 9.9902e-01 9.9988e-01 -9.9905e-01 -9.9978e-0I 9.6940e-01 

31) -9.9979e-01 -9.9980e-01 9.9922e-01 9.9986e-01 9.9916e-01 
Row 9. Columns 1 Thru 36 

I) 9.9978e-01 -9.9973e-01 -9.9425e-01 -9.9979e-01 9.4793e-0l ■ 

II) 7.6607e-01 -1.8215e-01 -9.9990e-01 -9.9974e-01 -9.9492e-01 

21) -9.9895e-01 -9.9987e-01 9.9898e-01 9.9978e-01 -9.695 le-01 
31) 9.9981e-01 9.9982e-01 -9.9921e-01 -9.9987e-01 -9.991 5e-0 1 

Row 10. Columns 1 Thru 36 

I) -9.9964e-01 9.9959e-01 9.9473e-01 9.9979e-01 -9.4997e-0I 

II) -7.6582e-01 1.8200e-01 9.9975e-0l 9.9960e-0I 9.9536e-01 

21) 9.9856e-01 9.9967e-01 -9.9856e-0L -9.9966e-01 9.6944e-01 

31) -9.9968e-01 -9.9970e-01 9.9885e-01 9.9967e-0I 9.9876e-0I 
Row 1 1 . Columns 1 Thru 36 

I) 7.5250e-0l -7.5179e-01 -6.9905e-01 -7.5412e-01 5.4647e-01 ■ 

II) I 0000e+00 -7.1435e-01 -7.7317e-01 -7.7935e-01 -8.2077e-01 

21) -7.7476e-01 -7.6670e-01 7.738Ie-01 7.6684e-0I -5.9718e-01 
31) 7.6653e-01 7.6609e-01 -7.5886e-01 -7.6493e-0I -7.5924e-01 

Row 12. Columns 1 Thru 36 

I) - 1 .622 Ie-0 1 1.5939e-01 1.0225e-01 1.6492e-01 5.0223e-02 

II) -7.1435e-01 l.OOOOe+OO 1.9580e-01 2.0415e-01 2.5665e*01 

21) i.784Oe-01 1.7672e-01 -1.7734e-Oi -1.7494e-01 -6.4138e-02 

31) -1.9090e-01 -1.9001e-01 1.8765e-01 1.8719e-01 1.8867e-01 


-9.9999e-0 1 -9.7551e-01 -9.9976e-01 9.9978e-01 -9.9964e-0l 
-9.9924e-01 -9.43I2e-0I 9.9905e-01 -9.9947e-01 9.9967e-01 
■9.9972e-01 9.9975e-01 -9.9973e-01 7.83l3e-01 1.3869e-02 
-9.9952e-01 

l.OOOOe+OO 9.7619e-0l 9.9971e-01 -9.9973e-01 9.9959e-01 
9.9916e-0I 9.4310e-01 -9.9895e-01 9.9954e-01 *9.997 le-01 
9.9968e-0 ! -9.9971e-01 9.9968e-01 -7.8530e-01 *1.5068e-02 
9.9938e-0l 

9.9589e-01 9.8395e-01 9.9409e-01 -9.9425e-0I 9.9473e-01 
9.9282e-0 1 9.2030e-01 -9.9210e-01 9.9320e-01 -9.9385e-01 
9.9381e-01 -9.9433e-01 9.9407e-0l -8.1058e-01 -8.9784e-02 
9.9483e-0 1 

9.9994e-01 9.7488e-01 9.9975e-01 -9.9979e-01 9.9979e-01 
9.9937e-01 9.4482e-01 -9.9913e-01 9.9946e-01 -9.9968e-01 
9.9968e-01 -9.998 le-01 9.9975e-0l -7.8143e-01 -U353e-02 
9.9958e-01 

-9.5308e-0I -9.7285e-01 -9.4745e-0I 9.4793e-01 -9.4997e-01 
-9.4442e-01 -8.3595e-01 9.4234e-01 -9.4538e-01 9.4702e-01 
-9 .465 6e -01 9.4828e-01 -9.4744e-01 8.5363e-01 2.7095e-01 
-9. 503 2e -01 

1.0000e+00 9.7610e-01 9.9973e-01 -9.9974e-01 9.9959e-01 
9.9918e-01 9.433 le-01 -9.9897e-01 9.9957e-01 -9.9974e-01 
9.9970e-0 1 -9.9972e-01 9.9969e-01 -7.8513e-01 -1.4196e-02 
9.9937e-01 

9.7610e-01 l.OOOOe+OO 9.7086e-01 -9.7096e-01 9.7087e-01 
9.669 le-01 8.6941e-01 -9.6530e-01 9.7294e-01 -9.7297e-01 
9.7074e-01 -9.7084e-01 9.7077e-01 -8.9989e-01 -2.092 le-01 
9.6889e-01 

9.9973e-0I 9.7086e-01 l.OOOOe+OO -9.9999e-01 9.9979e-01 
9.9980e-0 1 9.4872e-01 -9.9975e-01 9.9967e-01 -9.9982e-01 
9.9999e-01 -9.9997e-01 9.9996e-01 -7.71 lle-0 1 0.0000e+00 
9.9973e-0I 

-9.9974e-01 -9.7096e-01 -9.9999e-0I l.OOOOe+OO -9.9985e-01 
-9.998 le-01 -9.4877e-01 9.9974e-0I -9.9966e-01 9.9982e-01 
-9.999 8e -01 9.9999e-01 -9.9998e-01 7.7128e-Ol 0.0000e+00 
-9.997 5e -01 

9.9959e-0 1 9.7087e-0I 9.9979e-01 -9.9985e-01 l.OOOOe+OO 
9.9980e-01 9.5065e-0i -9.9958e-01 9.9952e-01 -9.9969e-01 
9. 997 2e -01 -9.999!e-01 9.9987e-01 -7.7124e-0l O.OOOOe+OO 
9.9965e-01 

-7.5237e-01 -6.071 le-01 -7.6648e-01 7.6607e-01 -7.6582e-01 
-7.7680e-01 -9.0551e-0I 7.7955e-01 -7.6740e-0l 7.655 le-01 
-7.6694e-0i 7.6661e-01 -7.6683e-01 2.5430e-01 -6.478 le-01 
-7.642 le-01 

1.5987e-01 -5.8122e-02 1.8253e-01 -1.8215e-0I 1.8200e-01 
1.9830e-01 3.95I2e-01 -2.0440e-01 1.737Ie-0I -1.739!e-0I 
1.8295e-01 -1.8267e-0l 1.8287e-01 4.7562e-01 9.0808e-01 
1.8869e-01 


Table (5.6) - The Cross Orthogonality of the First Krylov Block of the EPS Radiator 
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Row 13. Columns 1 Thru 36 

I) -9.994 ie-01 9.9931e-01 9.9328e-01 9.9947e-01 -9.4547e-01 9.9933e-01 9.6749e-01 9.9989e-01 -9.9990e-01 9.9975e-01 

II) -7.7317e-01 1.9580e-01 1.0000e+00 9.9995e-01 9.9574e-01 9.9993e-OJ 9.5135e-01 -9.9995e-01 9.9935e-01 -9.9955e-01 

21) 9.9864e-0I 9.9964e-01 -9.9866e-01 -9.9951e-01 9.6604e-01 9.9987e-01 -9.9989e-01 9.9987e-01 -7.6210e-01 1 .85 10e-02 

31) -9.9989e-0 1 -9.9989e-01 9.9931e-01 9.9989e-01 9.9927e-01 9.9982e-01 

Row 14, Columns 1 Thru 36 

1) -9.9906e-01 9.9896e-0l 9.9217e-01 9.9915e-01 -9.4255e-01 9.9898e-01 9.6537e-01 9.9975e-01 -9.9974e-0l 9.9960e-01 

11) -7.7935e-01 2.0415e-01 9.9995e-01 l.OOOOe+OO 9.9651e-01 9.9995e-0t 9.5402e-01 -l.OOOOe+OO 9.9925e-01 -9.994 le-01 
21) 9.9867e-01 9.9951e-0l -9.9867e-01 -9.9939e-01 9.6377e-01 9.9974e-01 -9.9975e-01 9.9973e-01 -7.5718e-01 2.8265e-02 

31) -9.9974e-01 -9.9973e-01 9.9903e-01 9.997!e-01 9.9899e-01 9.9963e-0l 

Row 15. Columns 1 Thru 36 

I) -9.9289e-01 9.9282e-01 9.8160e-0I 9.9343e-0l -9.2212e-01 9.9290e-0) 9.4835e-01 9.9490e-01 -9.9492e-01 9.9536e-01 

II) -8.2077e-01 2.5665e-0I 9.9574e-0i 9.9651e-01 l.OOOOe+OO 9.9660e4)l 9.7569e-01 -9.9650e-01 9.9525e-01 -9.9505e-01 

21) 9.9533e-01 9.9503e-01 -9.9515e-01 -9.9526e-01 9.4543e-01 9.9477e-0l -9.9518e-01 9.9500e-01 -7.2508e-01 9.7575e-02 

31) -9. 9459c -01 -9.9458e-01 9.9220e-01 9.9431e-01 9.9210e-01 9.9424c-01 

Row 16. Columns 1 Thru 36 

I) -9.9924e-01 9.9916e-01 9.9282e-01 9.9937e-01 -9.4442e-01 9.99 1 8e-0 1 9.6691e-01 9.9980e-01 -9.9981e-01 9.9980e-01 

11) -7.7680e-01 1.9830e-01 9.9993e-01 9.9995e-01 9.9660e-01 l.OOOOe+OO 9.5435e-0l -9.9994e-01 9.9946e-01 -9.9960e-01 

21) 9.9880e-01 9.9964e-01 -9.9879e-01 -9.9959e-01 9.6528e-01 9.9975e-01 -9.9986e-01 9.9981e-01 -7.6135e-01 2.4039e-02 

31) -9.9970e-01 -9.9971e-01 9.9884e-01 9.9967e-01 9.9877e-01 9.9961e-01 

Row 17. Columns 1 Thru 36 

I ) -9,4312e-01 9.4310e-01 9.2030e-01 9.4482e-01 -8.3595c-01 9.4331e-01 8.6941e-01 9.4872c-01 -9.4877e-01 9.5065e-01 

II) -9.055 le-01 3.9512e-01 9.5135e-01 9.5402e-01 9.7569e-01 9.5435e-01 l.OOOOe+OO -9.5400e-01 9.5099e-01 -9.4993e-01 

21) 9.5281e-01 9.4963e-01 -9.5217e-01 -9.5072e-01 8.6322e-01 9.4831e-01 -9.4967e-01 9.4907e-01 -6.0968e-01 2.8219e-01 

31) -9.4782e-01 -9.4779e-0l 9.4127e-01 9.4685e-01 9.4106e-0I 9.4686e-0t 

Row 1 8. Columns 1 Thru 36 

I) 9.9905e-0 1 -9.9895e-01 -9.9210e-01 -9.9913e-01 9.4234e-01 -9.9897e-01 -9.6530e-01 -9.9975e-0I 9.9974e-01 -9.9958e-01 

II) 7. 795 5e -01 -2.0440e-01 -9.9995e-01 -l.OOOOe+OO -9.9650e-01 -9.9994e-0l -9.5400e-01 l.OOOOe+OO -9. 9925e-01 9.9940e-01 

21) -9.9868e-01 -9.9951e-01 9.9869e-01 9.9939e-01 -9.6370e-01 -9.9974e-01 9.9974e-01 -9.9973e-01 7.5704e-01 -2.8583e-02 

31) 9.9973e-01 9.9972e-01 -9.9902e-01 -9.9970e-01 -9.9899e-01 -9.9962e-0l 

Row 19, Columns 1 Thru 36 

!) .9 9947e-01 9.9954e-01 9.9320e-01 9.9946e-01 -9.4538e-OI 9.9957e-01 9.7294e-0I 9.9967e-01 -9.9966e-01 9.9952e-01 

11) -7.6740e-01 1.7371e-01 9.9935e-01 9.9925e-01 9.9525e-01 9.9946e-01 9.5099e-01 -9.9925e-01 l.OOOOe+OO -9. 9997e-01 

21) 9.9970e-01 9.9994e-01 -9.9971e-0l -9.9998e-01 9.7105e-01 9.9966e-01 -9.9967e-01 9.9966e-01 -7.7909e-01 0.0000e+00 

31) -9.9897e-01 -9.9900e-01 9.9790e-01 9.991 le-01 9.9781e-01 9.9885e-01 

Row 20. Columns 1 Thru 36 

I) 9.9967e-01 -9.9971e-01 -9.9385e-01 -9.9968e-01 9.4702e-01 -9.9974e-01 -9.7297e-01 -9.9982e-01 9.9982e-01 -9.9969e-01 

II) 7.655 le-01 -1.7391e-01 -9.9955e-01 -9.9941e-01 -9.9505e-01 -9.9960e-01 -9.4993e-01 9.9940e-01 -9.9997e-01 l.OOOOe+OO 

21) -9.995 1 e-0 1 -9.9998e-01 9.9953e-01 1 .0000e+00 -9.7124e-Ol -9.9980e-01 9.9983e-01 -9.998 le-01 7.7823e-01 0.0000e+00 
31) 9.9927e-01 9.9930e-01 -9.9834e-01 -9.9939e-01 -9.9826e-01 -9.9918e-01 

Row 21. Columns 1 Thru 36 

I) -9.9854e-01 9.9867e-01 9.9044e-01 9.9843e-01 -9 3835e-01 9.9874e-01 9.7127e-01 9.9902e-01 -9.9895e-01 9.9856e-01 

II) -7.7476e-01 1.7840e-01 9.9864e-01 9.9867e-0I 9.9533e-01 9.9880e-0l 9.528 le-01 -9.9868e-01 9.9970e-01 -9.995 le-01 

21) l.OOOOe+OO 9.9953e-01 -l.OOOOe+OO -9.9957e-0l 9.6897e-01 9.9908e-01 -9.9892e-01 9.9896e-01 -7.7756e-01 2.0543e-02 

31) -9.9797e-01 -9.9799e-01 9.9670e-01 9.9818e-01 9.9663e-01 9.9774e-0I 

Row 22. Columns 1 Thru 36 

I) -9.9967e-0 1 9.9969e-0l 9.9359e-01 9.9965e-01 -9.4610e-01 9.9972e-0l 9.7230e-01 9.9988e-01 -9.9987e-01 9.9967e-01 

II) -7.6670e-01 1.7672e-01 9.9964e-01 9.9951e-01 9.9503e-01 9.9964e-01 9.4963e-0l -9.9951e-01 9.9994e-01 -9.9998e-0I 

21) 9.9953e-0I l.OOOOe+OO -9.9955e-0l -9.9997e-01 9.7060e-01 9.9988e-Ol -9.9986e-0l 9.9986e-01 -7.7626e-01 O.OOOOe+OO 

31) -9.9937e-0 1 -9.9939e-01 9.9854e-01 9.9950e-01 9.9846e-01 9.9927e-01 

Row 23. Columns 1 Thru 36 

I) 9.9860e-01 -9.9872e-01 -9.9057e-01 -9.9848e-01 9.3865e-01 -9.9879e-01 -9.7154e-01 -9.9905e-01 9.9898e-01 -9.9856e-01 

II) 7.738 le-01 -1.7734e-01 -9.9866e-01 -9.9867e-01 -9.9515e-01 -9.9879e-01 -9.5217e-01 9.9869e-01 -9.9971e-01 9.9953e-01 

21) -l.OOOOe+OO -9.9955e -01 l.OOOOe+OO 9.9958e-01 -9.6927e-01 -9.991 2e -01 9.9894e-01 -9.9899e-01 7.78l2e-01 -1.9052e-02 
31) 9.980 le-01 9.9802e-01 -9.9679e-01 -9.9823e-01 -9.9672e-01 -9.9779e-0I 

Row 24. Columns 1 Thru 36 

I) 9.9959e-01 -9.9964e-01 -9.9359e-0l -9.9961e-01 9.4636e-01 -9.9967e-01 -9.7271e-01 -9.9978e-01 9.9978e-01 -9.9966e-01 

II) 7.6684e-01 -1.7494e-01 -9.9951e-01 -9.9939e-01 -9.9526e-01 -9.9959e-01 -9.5072e-01 9.9939e-01 -9.9998e-01 l.OOOOe+OO 

21) -9.9957e-01 -9.9997e-0I 9.9958e-01 l.OOOOe+OO -9.709 le-01 -9.9976e-0l 9.9979e-01 -9.9977e-01 7.7785e-01 O.OOOOe+OO 

31 ) 9 99i9 e -01 9.9922e-01 -9.9819e-01 -9.9931e-01 -9.9810e-01 -9.9909e-01 

Table (5.6)(Cont.) - The Cross Orthogonality of the First Krylov Block of the EPS Radiator 
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Row 25. Columns I Thru 36 

I) -9.7430e-01 9.7493e-01 9.8418e-01 9.7366e-0l -9.7639e-01 9.7481e-01 9.9978e-01 9.6940e-01 -9.695 le-01 9.6944e-0l 

II) -5.9718e-01 -6.4138e-02 9.6604e-01 9.6377e-01 9.4543e-0l 9.6528e-01 8.6322e-01 -9.6370e-01 9.7105e-01 -9.7124e-01 

21) 9.6897e-01 9.7060e-01 -9.6927e-01 -9.709 lc-01 i.0000e+00 9.6926e-01 -9.6937c-01 9.6930e-01 -9.0022e-01 -2.2244e-01 

31) -9.6730e-01 -9.6753e-01 9.6765e-01 9.6825e-01 9.6734e-0l 9.6783e-01 

Row 26. Columns 1 Thru 36 

I) -9.9972e-0l 9.9968e-01 9.9381e-01 9.9968e-01 -9.4656e-01 9.9970e-0I 9.7074c-01 9.9999c-0I -9.9998e-01 9.9972e-0l 

II) -7.6694e-01 1.8295e-01 9.9987e-01 9.9974e-01 9.9477e-01 9.9975e-01 9.483Ie-01 -9.9974e-0I 9.9966e-01 -9.9980e-01 

21) 9.9908e-01 9.9988e-01 -9.991 2e-01 -9.9976e-01 9.6926e-01 I.0000e+00 -9.9994e-01 9.9996e-01 -7.709 le-01 O.OOOOe-fOO 

31) -9.9976e-01 -9.9976e-0) 9.9922e-01 9.9984e-0) 9.99l7e-01 9.9968e-01 


Row 27. Columns 1 ThraJ36 

I) 9.9975e-0! -9.997 le-01 -9.9433e-01 -9.9981 e-01 9.4828e-01 -9.9972e-01 -9.7084e-01 -9.9997e-01 9.9999e-01 -9.999 lc-01 

II) 7.666 le-01 -1.8267e-01 -9.9989e-01 -9.9975e-01 -9.9518e-01 -9.9986e-01 -9.4967e-01 9.9974c-0l -9.9967e-0I 9.9983e-01 

21) -9.9892e-01 -9.9986e-01 9.9894e-01 9.9979e-01 -9.6937e-01 -9.9994e-01 I.OOOOmOO -9.9997e-01 7.7105e-01 0.0000e+00 

31) 9.9979e-01 9.9981e-01 -9.991 le-01 -9.9983e-01 -9.9904e-01 -9.9974e-01 

Row 28. Columns I Thru 36 

I) -9.9973e-01 9.9968e-01 9.9407e-01 9.9975e-01 -9.4744e-01 9.9969c-01 9.7077e-0I 9.9996e-0I -9.9998e-01 9.9987e-01 

II) -7.6683e-01 1.8287e-01 9.9987e-01 9.9973e-01 9.9500e-01 9.9981e-01 9.4907e-01 -9.9973c-01 9.9966e-01 -9.998 le-01 

21) 9.9896e-01 9.9986e-01 -9.9899e-01 -9.9977e-01 9.6930e-01 9.9996e-01 -9.9997e-01 1.0000c+00 -7.7095e-01 0.0000e+00 

31) -9.9977e-01 -9.9977e-01 9.9912e-01 9.9983e-01 9.9906e-01 9.9970e-01 

Row 29. Columns 1 Thru 36 

I) 7.8313e-01 -7.8530e-01 -8.1058e-01 -7.8143e-01 8.5363e-01 -7.8513e-01 -8.9989e-01 -7.711 le-01 7.7128e-01 -7.7124e-01 

II) 2.5430e-01 4.7562e-01 -7.6210e-01 -7.5718e-01 -7.2508e-01 -7.6I35e-01 -6.0968e-0l 7.5704e-01 -7.7909e-01 7.7823e-01 

21) -7.7756e-01 -7.7626e-01 7.78I2e-01 7.7785e-01 -9.0022e-01 -7.709 1 e-01 7.7105e-01 -7.7095e-01 1.0000e+00 5.2942e-01 

31) 7.636 le-01 7.6423e-01 -7.6291e-01 -7.6626e-01 -7.6212e-01 -7.6457e-01 

Row 30, . C o l umn s i . i hmJP 

I) 1.3869e-02 -I.5068e-02 -8.9784e-02 -1.1353e-02 2.7095e-01 -1.4l96e-02 -2.0921e-01 0.0000e+00 0.0000e+00 0.0000e+00 

II) -6.478 le-01 9.0808e-01 1.8510e-02 2.8265e-02 9.7575e-02 2.4039e-02 2.8219e-01 -2.8583e-02 0.0000e+00 O.OOOOe-KX) 

21) 2.0543e-02 0.0000e+00 -1.9052e-02 0.0000e+00 -2.2244e-01 0.0000e+00 0.0000e+00 O.OOOOe+OO 5.2942e-0l ].0000e+00 

Row 31. Columns 1 Thru 36 

I) 9.995 le-01 -9.9937e-01 -9.9442e-01 -9.9956e-01 9.4900e-01 -9.9937e-01 -9.6844e-01 -9.9979e-01 9.9981e-01 -9.9968e-01 

II) 7.6653e-01 -1.9090e-01 -9.9989e-01 -9.9974e-01 -9.9459e-01 -9.9970e-01 -9.4782e-01 9.9973e-01 -9.9897e-01 9.9927e-01 

21) -9.9797e-01 -9.9937e-01 9.9801e-01 9.9919e-01 -9.6730e-01 -9.9976e-01 9.9979e-01 -9.9977e-01 7.6361e-01 0.0000e+00 

31) 1.0000e+00 1.0000e+00 -9.9970e-01 -9.9998e-01 -9.9966e-01 -9.9999e-01 

Row 32. Columns I Thru 36 

I) 9.9954e-0! -9.9941e-01 -9.9452e-01 -9.9959e-01 9.4928e-01 -9.9940e-01 -9.6868e-0I -9.9980e-01 9.9982e-01 -9.9970e-0I 

II) 7.6609e-01 -1.9001e-01 -9.9989e-01 -9.9973e-01 -9.9458e-0I -9.997 le-01 -9.4779e-01 9.9972e-01 -9.9900e-01 9.9930e-01 

21) -9.9799e-01 -9.9939e-01 9.9802e-01 9.9922e-01 -9.6753e-01 -9.9976e-0I 9.9981e-01 -9.9977e-01 7.6423e-01 0.0000e+00 

31) 1.0000e+00 1.0000e+00 -9.9969e-01 -9.9998e-0l -9.9964e-01 -9.9999e-01 

Row 33. Columns 1 Thru 36 

1) -9. 990 8e -01 9.9888e-01 9.9475e-01 9.9901e-01 -9.5089e-0I 9.9886e-01 9.6838e-0I 9.9922e-01 -9.9921e-0L 9.9885e-01 

11) -7.5886e-01 1.8765e-01 9.9931e-01 9.9903e-01 9.9220e-01 9.9884e-01 9.4127e-01 -9.9902e-01 9.9790e-0I -9.9834e-01 

21) 9.9670e-01 9.9854e-01 -9.9679e-01 -9.981 9e-01 9.6765e-01 9.9922e-0 1 -9.991 le-01 9.99 12e-0 1 -7.629 le-01 0.0000e+00 

31) -9.9970e-01 -9.9969e-01 1.0000e+00 9.9971e-01 1.0000e+00 9.9975e-01 
Row 34. Columns 1 Thru 36 

I) -9.9963e-0l 9.995Ie-01 9.9450e-01 9.9963e-01 -9.4899e-01 9.995le-01 9.6942e-01 9.9986e-01 -9.9987e-01 9.9967e-01 

II) -7.6493e-01 1.8719e-01 9.9989e-01 9.9971e-01 9.9431e-01 9.9967e-01 9.4685e-01 -9.9970e-01 9.991 lc-01 -9.9939e-0I 

21) 9.9818e-01 9.9950e-01 -9.9823e-01 -9.9931e-01 9.6825e-01 9.9984e-01 -9.9983e-01 9.9983e-01 -7.6626e-01 0.0000e+00 

31) -9.9998e-01 -9.9998e-01 9.9971e-01 1.0000e+00 9.9967e-01 9.9997e-01 

Row 35. Columns 1 Thru 36 

I) -9.9900e-0I 9.9879e-0I 9.9456e-01 9.9892e-01 -9.5045e-0I 9.9878e-01 9.6805e-01 9.9916e-01 -9.9915e-0t 9.9876e-01 

II) -7.5924e-01 1.8867e-01 9.9927e-01 9.9899e-01 9.9210e-01 9.9877e-0l 9.4106e-01 -9.9899e-0I 9.978 1 e-0 1 -9.9 826e -01 

21) 9.9663e-0l 9.9846e-0I -9.9672e-01 -9.9810e-01 9.6734e-01 9.9917e-01 -9.9904e-01 9.9906e-01 -7.621 2e-0 1 0.0000e+00 

31) -9.9966e-OI -9.9964e-01 l.0000e+00 9.9967e-01 1.0000e+00 9.997 le-01 

Row 36. Columns 1 Thru 36 

I) -9.9952e-01 9.9938e-01 9.9483e-01 9.9958e-01 -9.5032e-01 9.9937e-01 9.6889e-01 9.9973e-01 -9.9975e4)l 9.9965e-01 

II) -7.642 le-01 1.8869e-01 9.9982e-01 9.9963e-01 9.9424e-01 9.9961e-01 9.4686e-01 -9.9962e-01 9.9885e-01 -9.991 8e -01 

21) 9.9774e-01 9.9927e-01 -9.9779e-01 -9.9909e-01 9.6783e-01 9.9968e4)l -9.9974e-01 9.9970e-01 -7.6457e-01 0.0000e->00 

31) -9.9999e-0l -9.9999e-01 9.9975e-01 9.9997e-01 9.997 le-01 l.OOOOe+OO 


Table (5.6)(Cont.) - The Cross Orthogonality of the First Krylov Block of the EPS Radiator 
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Row 1 , Columns I Thru 8 

1.0000e+00 9.5350e-01 -9.755 le-01 7.5250e-01 -1.6221e-0I -9.4312e-01 7.83I3e-01 1.3869e-02 
Row 2. Columns 1 Thru 8 

9.5350e-01 1 0000e+00 -9.7285e-01 5.4647e-01 5.0223e-02 -8.3595e-01 8.5363e-01 2.7095e-01 
Row 3. Columns 1 Thru 8 

-9.755 le-01 -9.7285e-01 1.0000e+00 -6.071 le-01 -5.8122e-02 8.6941e-01 -8.9989e-01 -2.0921e-01 


Row 4, Columns 1 Thru 8 

7.5250e-01 5. 4647e-01 -6.071 le-01 1.0000e+00 -7.1 435e-01 -9.055 le-01 2.5430e-01 -6.478 le-01 


Row 5. Columns 1 Thru 8 

-1.622 le-01 5.0223e-02 -5.81 22e-02 -7.1 435e-01 1.0000e+00 3.9512e-01 4.7562e-01 9.0808e-01 


Row 6. Columns 1 Thru 8 

-9.4312e-01 -8.3595e-01 8.6941e-01 -9.0551e-01 3.9512e-01 1.0000e+00 -6.0968e-01 2.8219e-01 


Row 7, Columns 1 Thru 8 

7.8313e-01 8.5363e-01 -8.9989e-01 2.5430e-01 4.7562e-01 -6.0968e-01 1.0000e+00 5.2942e-01 


Row 8. Columns 1 Thru 8 

1.3869e-02 2.7095e-01 -2.0921e-01 -6.4781e-01 9.0808e-01 2.8219e-01 5.2942e-01 1.0000e+00 


Table (5.7) - The Cross Orthogonality of the First Krylov Block 
of the EPS Radiator, After Block Filtering 
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Modes 

Finite 

Element 

Ritz 

Vectors 

Percent 

Difference 

Normal 

Modes 

Percent 

Difference 

1st Bend 

.5464 Hz 

.5464 Hz 

0 . % 

.5464 Hz 

0 . % 

2nd Bend 

3.415 Hz 

3.415 Hz 

0 . % 

3.415 Hz 

0 . % 

3rd Bend 

9.522 Hz 

9.610 Hz 

.92 % 

9.522 Hz 

0 . % 

4th Bend 

18.57 Hz 

24.24 Hz 

31. % 

18.57 Hz 

0 . % 

1st Tors 

27.58 Hz 

27.58 Hz 

0 . % 

27.58 Hz 

0 . % 


Table (5.8) - Comparison Between the Frequencies of a 
Cantilevered Beam Obtained From Fixed Interface Component 
Normal Modes, Ritz Vectors, and Finite Elements 
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Modes 

Finite 

Element 

Ritz 

Vectors 

Percent 

Difference 

1st Bend 

.5464 Hz 

.5464 Hz 

0. % 

2nd Bend 

3.415 Hz 

3.415 Hz 

0. % 

3rd Bend 

9.522 Hz 

9.524 Hz 

.02 % 

4th Bend 

18.57 Hz 

19.23 Hz 

3.6 % 

1st Tors 

27.58 Hz 

27.58 Hz 

0. % 


Table (5.9) - Comparison Between the Frequencies of a 
Cantilevered Beam Obtained From Free Interface Component 
Ritz Vectors, and Finite Elements 
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Modes 

Finite 

Element 

Ritz 

Vectors 

Percent 

Difference 

1st Rigid 

0. Hz 

0. Hz 

0 . % 

2nd Rigid 

0. Hz 

0. Hz 

0 . % 

3rd Rigid 

0. Hz 

0. Hz 

0 . % 

4th Rigid 

0. Hz 

0. Hz 

0 . % 

5th Rigid 

0. Hz 

0. Hz 

0 . % 

6th Rigid 

0. Hz 



1st Bend 

3.474 Hz 

3.474 Hz 

0 . % 

2nd Bend 

9.549 Hz 

9.549 Hz 

0 . % 

3rd Bend 

18.65 Hz 

18.69 Hz 

.21 % 

4th Bend 

30.70 Hz 

33.15 Hz 

8. % 


Table (5.10) - Comparison Between the Frequencies of a 
Free-Free Beam Obtained From Free Interface Component 
Ritz Vectors, and Finite Elements 















Various Normalization Methods 


Various 

Orthogonalization 

Methods 

Orthonormal ized 
wrt the Euclidean 
Norm 

Orthonormalized 
wrt the Mass 
Matrix, and all 
previous blocks 

Orthonormalized 
wrt the Mass 
Matrix, previous 
two blocks only 

Cholesky/QR 

Decomposition 

97.0 

140.9 

112.4 

Gram-Schmidt, 
Reorthogonalize 
with all vectors 

103.9 

379.2 

385.6 

Gram-Schmidt, 
Reorthogonalize 
with selected 

116.0 

409.1 

468.9 

Gram-Schmidt, 
Reorthogonalize 
using modified 
Gram-Schmidt 

461.0 

3747.0 

4196.9 


Various Eigenvalue Extraction Methods 


Normal 

Lanczos 

Modified Givens 

Inverse Power 

Eigenvalues 

115.5 

5758.8 

6455.5 


Table (5.11) - Comparisons of the Computer Time (seconds) Required 
to Form a Component of the EPS Radiator Finite Element Model 
(4000 DOF), Using Selected Orthonormalization Options 























Various Normalization Methods 


Various 

Orthogonalization 

Methods 

Orthonormalized 
wrt the Euclidean 
Norm 

Orthonormalized 
wrt the Mass 
Matrix, and all 
previous blocks 

Orthonormalized 
wrt the Mass 
Matrix, previous 
two blocks only 

Cholesky/QR 

Decomposition 

323.2 

399.5 

376.1 

Gram-Schmidt, 
Reorthogonalize 
with all vectors 

341.6 

534.3 

542.1 

Gram-Schmidt, 
Reorthogonalize 
with selected 

348.5 

558.3 

565.2 

Gram-Schmidt, 
Reorthogonalize 
using modified 
Gram-Schmidt 

559.5 

3062.5 

3095.3 


Various Eigenvalue Extraction Methods 


Normal 

Lanczos 

Modified Givens 

Inverse Power 

Eigenvalues 

743.7 




Table (5.12) - Comparisons of the Computer Time Required 
to Form a Component of the Cassini Spacecraft Finite Element 
Model (11,100 DOF), Using Selected Orthonormalization Options 
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Figure (5.1) - Space Station Electrical Power System Deployed Radiator 
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Figure (5.2) - Space Station EPS Radiator Undeformed Finite Element Model 
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Figure (5.3) - EPS Radiator First Normal Mode, Fixed Interface 
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Figure (5.4) - EPS Radiator Second Normal Mode, Fixed Interface 
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Figure (5.5) - EPS Radiator Third Normal Mode, Fixed Interface 
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Figure (5.6) - Cassini Spacecraft 
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Figure (5.7) - Cassini Spacecraft Finite Element Model 
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(a) Initial Calculation 



(b) After Gram-Schmidt Orthogonalization 


Figure (5.8) - Fixed Interface Ritz Vector, 
Block One, Vector One 
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(a) Initial Calculation 



(b) After Gram-Schmidt Orthogonalization 


Figure (5.9) - Fixed Interface Ritz Vector, 
Block One, Vector Two 


(a) Initial Calculation 



(b) After Gram-Schmidt Orthogonalization 


Figure (5.10) - Fixed Interface Ritz Vector, 
Block Two, Vector One 
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(a) Initial Calculation 



(b) After Gram-Schmidt Orthogonalization 


Figure (5.11) - Fixed Interface Ritz Vector, 
Block Two, Vector Two 



Chapter 6 
Targeted Shifting 

6.1) Introduction: 

Inverse operators have been widely used in a variety of eigenvalue extraction 
routines. The use of inverse operators is commonly called sequence shifting. By 
shifting the dominant frequency sought by an algorithm, an eigenvalue extraction 
routine can locate missing eigenvalues within a specific frequency range. In Lanczos 
eigenvalue extraction, the inverse operator shifts the dominant frequency of the Krylov 
sequence to a range of frequency where missing eigenvalues are located. By doing 
so, the number of eigenvalues beyond the range of interest is minimized while 
insuring that interesting modes are calculated. 

The utility made of the ability of the inverse operator to shift the dominant 
frequency of the Krylov sequence is not limited to searching for missing eigenvalues. 
This chapter presents shifting that is targeted to the frequency of the dominant applied 
dynamic load vector, rather than to missing eigenvalues. By shifting the dominant 
frequency sought by the Krylov sequence, static Ritz vectors which can contribute to 
an accurate dynamic response prediction are generated. The use of the inverse 
operator in this manner has been identified in this work as targeted shifting. 
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6.21 Theory: 

As discussed in section 2.5, an inverse operator can direct a solution to 
particular frequency range. The problem discussed in section 2.5 was the computation 
of selected eigenpairs of the generalized eigenvalue problem, 

(K-XM)x = 0 ( 61 ) 

with X being the eigenvalues and jc the eigenvectors, the matrix A of equation (2.46) 
being re-written to the common structural dynamics usage of K. A system is created 
using the operator 


(x-oiifr'M (6 - 2) 

which has the same eigenvectors as (6.1). The shifted system eigenvalues are 
transformed to l/(k - a), with a being defined as the shift frequency. When the 
inverse operator is applied the eigenvalue nearest a becomes the dominant one and 
the sequence will converge to the corresponding eigenvector. The cost for performing 
the shift is the factorization of (K - oM ) ‘. 

If the selected shift frequency, o, is coincidentally an exact eigenvalue, then 
(K - oM)' 1 is, by definition, singular. To obtain a non-singular matrix the shift value 
is merely altered by a small constant value, such as .1, to move the shift away from 
the eigenvalue frequency. 
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To add an inverse operator to the static Ritz vector calculation the Krylov 
sequence presented in chapter 2 must be modified. Several different inverse operators 
could be used. Some inverse operators, such as that defined by equation (6.2), 
transform the system equations so that they produce transformed eigenvalues, but the 
produced eigenvectors are not transformed. Others inverse operators transform the 
system equations such that both the eigenvalues and eigenvectors are transformed. (In 
eigenvalue extraction routines this transformation necessitates a back transformation 
to obtain a final solution.) The operator defined by equation (6.2) was selected for 
this work because this operator does not transform vectors resulting from the Krylov 
sequence when shifting is applied. 

When this operator is applied to the Krylov sequence presented in Table (4.1) 
and chapter 2, equation (2.8), 





m 


it Qj - 1 


( 2 . 8 ) 


is revised to become 


q'j =(k ir om u y 1 m ii q j _ 1 


(6.3) 


Note that if a shift frequency is defined as zero, then equation (6.3) becomes equation 
(2.8). The standard Krylov sequence can be considered to be a shifted sequence with 
the shift frequency permanently defined as zero. 
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If the inverse operator is applied to the equation used to generate the initial 
block of Ritz vectors, equation (2.6) 


9 1 




( 2 . 6 ) 


would be revised to become 

«r = ( k u - ° w s ) _l m ti( m u*ic +m ic) 


(6.4) 


6.3) Illustrative Example: 

A short numerical example is given to illustrate the effectiveness of shifting. 
Consider a system with the following defining matrices. 


and 



1 0 
2 1 
1 2 


(6.5) 


2 0 0 

M = 0 10 (66) 

0 0 1 

In this example, a load, rather than the boundary flexibility matrix, will be used to 
initiate the sequence. The initial load dependent static Ritz vector is defined by the 
following equation 


q ;'=K l P 


(6.7) 
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If p is defined as 

p = [ 11 If ( 6 - 8 > 

then the normalized initial Ritz vector, q„ is 

q x = [.408 .408 .816f (6 9) 

The Rayleigh quotient of a vector derived from a two matrix, or general 
system, can be defined as 


Pi = 


tfMq. 


( 6 . 10 ) 


When the vector used with equation (6.10) is an eigenvector, then the Rayleigh 
quotient is the eigenvalue associated with that eigenvector. The square root of a 
vector’s Rayleigh quotient implies, as a eigenvalue does, a frequency at which the 
vector is most likely to respond. The initial static Ritz vector, q,, has a Rayleigh 
quotient of p = 2.857. 

The eigenvalues of the system defined in equations (6.5) and (6.6) are .775, 
2., and 3.225. The system eigenvectors, <j>, are 


* = 


-.309 .707 .312 
.756 0.000 .764 
-.617 -.707 .624 


( 6 . 11 ) 
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The cross orthogonality of the third eigenvector and the first Ritz vector, <j> 3 T </ , , is .948. 
The initial static Ritz vector, produced by equation (6.8), is most numerically similar 
to the third eigenvector, both by mode and by frequency. 

The shifted initial load dependent static Ritz vector is defined by the following 
equation 

q” = (K-oMY 1 Mp ( 612 ) 

Using the same force vector, p, defined by equation (6.8), and a shift with an arbitrary 
frequency of 1.0, (an assumption being made that the force vector,/?, has a dominant 
frequency component of 1.0) the normalized initial Ritz vector, q,, is 

^=[0.000 .894 - 447] r ( 613 ) 

The Rayleigh quotient of this vector is p = 1.199, relatively close to the shift value 
of 1.0. The cross orthogonality of the first eigenvector and the first Ritz vector, 
is .952. The shifted initial static Ritz vector’s frequency content is closer to a target 
frequency than that of the vector (equation 6.9) produced without the shift. Rather 
than resembling the third eigenvector of the system, the shifted initial Ritz vector 
resembles the first eigenvector. 

6.4) Targeted Shifting: 

After considering the example presented in section 6.3, it is reasonable to 
consider using the frequency content of an applied dynamic loading in the shift 
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methodology. This can be accomplished by determining the frequency content of the 
applied dynamic load and using the dominant frequencies as shift values. The 
resulting shifted static Ritz vectors will be more likely to accurately predict system 
responses than static Ritz vectors which have a frequency content far removed from 
that of the applied dynamic load. In addition, it is also reasonable to augment the 
static Ritz vectors derived from the boundary flexibility matrices with static load 
dependent Ritz vectors. 

Each separate targeted shift that is performed entails the computation cost of 
decomposing (k u - am,,)' 1 into the upper and lower triangular factors V and L. In a 
large finite element model the computational cost to perform this factorization can be 
significant and over-aggressive shifting can make the computational cost of this 
procedure exceed that of the orthonormalization and make the static Ritz vector 
calculation inefficient. As a result it is beneficial to minimize the amount of shifts in 
the overall targeted shifting strategy. The overall targeted shifting strategy will 
discussed in the next section. 

Since relatively few frequencies can be targeted in a computationally efficient 
shift strategy, the method chosen to select frequencies important to the dynamic 
response should seek, at most, several dominant frequencies. The absolute magnitude 
of any particular frequency is unimportant since the relative magnitude of a frequency 
is used to identify the dominant frequencies. The choice of which method is used to 
determine the frequency content of the time domain dynamic loading can be made 
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somewhat subjectively. Most analytical procedures which transform a time domain 
response into the frequency domain, such as calculating a power spectral density, or 
a frequency response function, etc. would be sufficient for identifying several 
dominant frequencies. Alternately, if frequency information is available describing the 
dynamic environment from non-analytical experience, then that information can also 
easily be used. A shift can also easily be targeted to test derived frequency 
information describing the structure itself. There is, of course, a possibility in some 
applications, due to the nature of the dynamic excitation, that no dominant frequencies 
can be identified. In those circumstances, the unshifted Krylov sequence is sufficient. 

Static Ritz vectors derived from the spatial distribution of the applied dynamic 
load, in conjunction with targeted shifting, can augment the static Ritz vectors derived 
from the boundary flexibility matrices. A load dependent static Ritz vector may be 
calculated from a dynamic load by using the spatial distribution of the dynamic load 
at a single time step to create a representative static load. Selecting a time step where 
the applied dynamic load, at a particular node, is at a maximum, is one possible 
criteria. Consideration should be given, in the time step selection, to the dominant 
frequency chosen for shifting. Another possible method for creating a representative 
static load derived from the applied dynamic load would be to, determine the 
maximum applied load over all time steps for each node, and from these maximums 
synthesize a single static force vector. This representative static vector may have no 
physical relation to the dynamic spatial characteristics of the applied loading, and the 
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resulting static Ritz vector may not resemble a response to dynamic excitation. Once 
a representative static load has been created, static Ritz vectors may be calculated by 
using equation (6.12) with the selected dominant frequency. This load dependent 
static Ritz vector is appended to the vectors calculated by the algorithm given in 
Tables (4.1) and (4.2). 

6.51 Shift Strategy: 

A definite theory does not exist which describes an optimal ordering and value 
of shifts in Krylov sequence methodologies. In Lanczos eigenvalue extraction a fairly 
complicated heuristic approach, which is successful in extracting eigenvalues, has been 
developed 13 . That particular strategy has limited applicability to the creation of 
boundary flexibility method static Ritz vectors. As a result, an alternate heuristic 
strategy has been developed specifically to create accurate static Ritz vectors as early 
in the Krylov sequence as possible. 

The new shift strategy for static Ritz vectors is introduced as follows. No shift 
(or equivalently, a shift of zero) is applied to the initial block of vectors. The 
representative static load, derived from applied the dynamic loading, is also not 
applied to calculate the initial block. Therefore, equation (2.6), as previously 
presented, is used to generate the initial block of Ritz vectors 


?r=v 1 ( m «*ic +m «c) 


( 2 . 6 ) 
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The non-zero shift equation, (6.4), is not used in the computer implementation of this 
work. The initial zero shift is performed because a basic set of vectors targeted to the 
lowest frequencies is desired. These vectors are desired because, with only a limited 
amount of shifts targeting the several selected dominant frequencies, they describe 
other important dynamic characteristics not in the frequency spectrum of the targeted 
shifts. 

Two selected dominant frequencies are allowed, and a cutoff frequency is also 
used as a shift frequency, for a total of three non-zero targeted shifts. When 
approximately one quarter of the total required vectors have been calculated (the next 
chapter presents a method to determine this value) a shift targeted to the first selected 
frequency is performed. At this shift the representative static load, derived from the 
applied dynamic loading, is also used by appending it to the previous vector block. 
The resulting modified Krylov sequence equation (6.3) becomes 

q'j =(k ii -o k m u r 1 m ii [ p k ] (614) 

with the addition of the representative static load, p k , and with the subscript k 
indicating the shift number, initially k = 1. The vector block, q*, will have one more 
column in it than the vector block, q H , due to the augmentation of q^, with the 
representative static load, p k . Subsequent to the block increment, j, at which the shift 
is applied, equation (6.3) is used to calculate additional vector blocks. When 
approximately one half of the total required vectors have been calculated a second 
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shift, k = 2, is applied with the second selected target frequency. The third shift is 
applied when approximately three quarters of the vectors required have been 
generated. The frequency used for the shift is a user defined cutoff frequency which 
will be discussed in the next chapter. 

This strategy was developed through trial and error and is as a result heuristic. 
There may be other strategies which would be successful at creating a set of static 
Ritz vectors which accurately represent a substructure. The effectiveness of the 
presented strategy will not be demonstrated until chapter 8, after a sequence 
termination methodology has also been presented. 

6.61 Applied Dynamic Loading of Example Structural Models: 

Several examples applied to the structural examples given in chapter 5 are 
presented. For the beam example the time domain loading is shown in Figure (6.1). 
The response spectrum of this time domain load is shown in Figure (6.2). The 
dominant frequency selected is the frequency where the response spectrum reaches it’s 
peak, at .546 HZ. Because the time domain dynamic load is a rectangular impulse, 
this frequency is equal to the inverse of twice the impulse width. The dynamic 
loading is applied to a single node, the end point. As a result, the maximum spatial 
load for the structure is the maximum applied dynamic load at that node and a 
representative static load is easily created. 

For the EPS Radiator the spatial distribution of the dynamic loading is the 
same at every time step that the load is applied. As a result, all choices for an 
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appropriate time step from which to derive a representative static load are identical. 
The load is applied in a series of rectangular impulses as shown in Figure (6.3). The 
response spectrum of this applied dynamic load is shown in Figure (6.4). The peak 
of the response spectrum, at a frequency of approximately 5 HZ, is a result of the 
pulse width of .08 sec of the applied dynamic load. This applied dynamic loading is 
derived from a plume impingement loading event. Plume impingement loading occurs 
when the space shuttle approaches the space station, directing it’s reaction control 
system jets toward the station. It is the critical loading event for the space station 
structure in it’s on-orbit configuration. The specific loading case shown in Figure 
(6.3) is a result of the shuttle translating in yaw along the main station axis. 

All Cassini spacecraft loading comes through the spacecraft/launch vehicle 
interface. As a result, the boundary flexibility vectors are the only available loading 
vectors. (A separate applied dynamic load vector does not exist.) Figures (6.5), (6.6), 
and (6.7) show the interface acceleration for a liftoff loading event in the three 
translational directions. Response spectra of this inertia loading are shown in Figures 
(6.8), (6.9), and (6.10). Dominant frequencies of 5 HZ and 13 HZ were selected for 
the two targeted shifts. The 13 HZ frequency was derived from the pitch and yaw 
excitation response spectrum and the 5 HZ frequency was selected from the 
longitudinal excitation response spectrum. No representative static load was appended, 
because of the lack of a spatial distribution for the applied dynamic load. 
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Figure (6.1) - Time Domain Beam Model Applied Dynamic Loading 
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Figure (6.2) - Response Spectrum of the Beam Model Applied Dynamic Loading 
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Figure (6.3) - Plume Impingement Dynamic Loading on the EPS Radiator 
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Figure (6.4) - Response Spectrum of the Plume Impingement Dynamic Loading 
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Figure (6.5) - Cassini Spacecraft Pitch Translation Interface Acceleration 
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Figure (6.6) - Cassini Spacecraft Yaw Translation Interface Acceleration 



MtO>W H> WO^O'fl 


124 



TIME (SECS) 

ASSINI LIFTOFF 

S0U7TICN 69, FORCE RESPONSE SPECTRA 


Figure (6.7) - Cassini Spacecraft Longitudinal Interface Acceleration 
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Figure (6.8) - Response Spectrum of the Pitch Interface Acceleration 
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Figure (6.9) - Response Spectrum of the Yaw Interface Acceleration 




MCn>Cd H> OWHt 4 *tjro> 


127 



CASSINI LIFTOFF 

SOLUTION 69, FORCE RESPONSE SPECTRA 


Figure (6.10) - Response Spectrum of the Longitudinal Interface Acceleration 



Chapter 7 

Krylov Sequence Termination Techniques 

7.D Introduction: 

Determining when a sufficient number of static Ritz vectors have been 
calculated is a problem which all static Ritz vector methods, including the boundary 
flexibility method, share. When enough static Ritz vectors have been calculated to 
accurately represent the dynamics of the given finite element model the Krylov 
sequence can be terminated. A heuristic method for sequence termination has been 
developed which is based, not on rigorous mathematics, but on the observed properties 
of load dependent and boundary flexibility static Ritz vector creation. 

A mathematically rigorous basis for judging when the sequence can be 
terminated was extensively sought for, both in the literature and by analytical 
investigation. It is possible that such a solution exists, however, a rigorous solution 
was neither located nor could be created here. It is also possible that a mathematically 
rigorous sequence termination methodology, which would be independent of physical 
knowledge of typical structural dynamic systems (this knowledge was used to create 
the heuristic method presented in this chapter), cannot be created. 

The mathematically rigorous solution sought was analogous to modal 
truncation. Modal truncation, based upon an eigenvalue cutoff, is the most popular 
basis for determining if a sufficient number of normal eigenvectors have been 
calculated. Static Ritz vectors do not have an eigenvalue with which to associate a 
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truncation limit, but they do have a Rayleigh value, as presented in chapter 6. Using 
Rayleigh values as a partial substitute for eigenvalues was explored and was found to 
not be practical. 

7.2) Error Criteria and Effective Mass: 

Several error criteria have been proposed which truncate the Krylov sequence 
when a somewhat arbitrary variable reaches a arbitrary value. These proposed criteria 
were briefly described in chapter 2. Most of the proposed termination criteria are 
applicable only to load dependent static Ritz vectors, with a goal of obtaining an 
accurate solution to the given static problem. Error criteria which are completely 
aimed toward the solution of static problems have a limited direct utility in their 
application to dynamic substructuring and boundary flexibility static Ritz vectors. 

The single proposed error criteria for Krylov sequence termination which is 
directed to dynamic problems is that proposed by Yiu and Landress 30 . This criteria 
is based on a parameter, effective mass, commonly used in structural dynamics for the 
identification of globally important normal eigenvectors of a structure. Effective mass 
is a measure of the amount of the total structural mass represented in each individual 
eigenvector. Effective mass is calculated as follows, beginning with an intermediate 
matrix, M ER , being defined as 

M ER = * I M<b liB 


(7.1) 
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where O is the matrix of eigenvectors, and is the rigid body transformation matrix 
of the structure. Each term within the resulting matrix of equation (7.1) is squared to 
form the effective mass of each eigenvector. 

Effective Mass = < M m > 2 (7.2) 

Note that eigenvalues are not used in the effective mass calculation, allowing static 
Ritz vectors to directly replace the eigenvectors in equation (7.1). Each eigenvector 
has a total of six effective mass values, three for the translations and three for the 
rotations. The individual effective mass terms of the eigenvectors can, separately in 
each translation and rotation coordinate, be summed and the sum can be compared to 
the rigid body mass of the structure. The percentage of mass represented by 
eigenvectors, in each direction, is referred to as the total effective mass. 

Yiu and Landress proposed that the effective mass calculation can be made 
using static Ritz vectors, and that when the total effective mass reaches an arbitrary 
percentage, the sequence can be terminated. The recommended arbitrary percentage 
was a minimum of 90% total effective mass, which is consistent with standard 
aerospace practice when using normal eigenvectors. In standard aerospace practice 
this cutoff value is used to determine if the important dynamics of a structure has been 
test verified, and in some cases even a greater percentage than 90% is specified. 

The primary assumption made when using total effective mass as an error 
criteria is that, when enough vectors have been calculated to achieve the selected 
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arbitrary percentage of represented mass, an accurate solution is guaranteed. 
Engineering experience has demonstrated that if this percentage is set at a high enough 
value an accurate dynamic solution will often be possible. However, usually an 
accurate solution is possible with far fewer vectors than the amount required to reach 
the cutoff value. The effective mass of the normal eigenvectors, from the beam 
example model presented in chapter 5, is shown in Table (7.1). The effective mass 
of the normal eigenvectors, from the EPS Radiator and the Cassini finite element 
models, is shown in appendix A. These are typical effective mass results for typical 
finite element models. The total effective mass of these models does not converge 
monotonically, nor does it reach any particular percentage within the frequency range 
of interest. Table (7.1) demonstrates that, for the torsional rotation (R2) of the beam 
model, eigenvectors with eigenvalues up to 150 HZ yield a total effective mass of 
76.7%. The contribution of the 22nd mode at 133.73 HZ was required to bring the 
axial translation (T2) total effective mass above 90%. Eigenvectors with eigenvalues 
in this frequency range are not required to allow an accurate dynamic solution for the 
excitation described in chapter 6. The rectangular impulse excitation has a relatively 
broad frequency spectrum, and excitation at a particular high frequency would require 
eigenvectors in that range for an accurate solution, but in any case a required total 
effective mass of 90% is arbitrary. 

For the EPS Radiator model, there are 175 modes below 140 HZ. These 
modes yield a total effective mass of 78.1% for axial translation and 75.6% for out- 
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of-plane translation. The in-plane rotational total effective mass for this set of modes 
is 75.7%. If the desired total effective mass was set at a value of 90%, hundreds of 
vectors would be generated, but most are not required for an accurate dynamic 
solution, at frequencies of interest in loads calculation. The Cassini model total 
effective mass reaches 90% much quicker than the example beam and the EPS 
Radiator. Even here, if the desired value is set higher than 90%, many unnecessarily 
vectors might be calculated. In conclusion, effective mass, while a useful tool, was 
investigated but was not implemented in this work. Other similar error criteria were 
also found to not exhibit rapid, monotonic or well behaved convergence. 

7.31 Modal Density Truncation: 

A new method has been developed that is based upon the density of the modal 
space in a finite element model and the observation that the most dynamically 
significant static Ritz vectors are generated near the beginning of the Krylov sequence. 
This method was developed from the observation that almost all of the boundary 
flexibility static Ritz vectors calculated early in the Krylov sequence are required to 
obtain an accurate solution to the dynamic response problem. This includes static Ritz 
vectors with Rayleigh values which are high relative to the frequency range of interest, 
which might erroneously imply relatively low importance to dynamic response 
prediction. It was also observed, in general, for the test examples utilized and from 
the current literature, that accurate dynamic response predictions are possible with 
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fewer load or boundary flexibility derived static Ritz vectors than with normal 
eigenvectors. 

Krylov sequence termination using modal density truncation can be described 
in detail as follows. First, a cutoff frequency is defined by the user. This frequency 
should be based upon the dynamic range of interest and the frequency spectrum of the 
forcing function. Next, the number of normal eigenvectors with eigenvalues below 
the cutoff frequency is determined, which is referred to here as the modal density. 
The actual eigenvectors and eigenvalues of the model are not required, just their 
number. The Krylov sequence is then initiated and when the number of static Ritz 
vectors created is equal to a fraction of the modal density, the sequence is terminated 
and the component is formed. A value for this fraction has been determined by trial 
and error. The effectiveness of this sequence termination strategy will be 
demonstrated in the next chapter. 

Modal density truncation is not mathematically rigorous, but it is practical. 
It is effective because the order of the system, within the frequency range of interest, 
is related to the number of static Ritz which are needed to represent this model. This 
algorithm is simple, and does not pursue unobtainable arbitrary parameter values, 
which can be the case when using error functions. The computational cost of this 
method is the factorization of the mass and stiffness matrices, shifted to the cutoff 
frequency. The factorization for determining the modal density can be reused as a 
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shift in the Krylov sequence, targeting static Ritz vectors near the cutoff frequency 
at no additional computational cost. 

7.41 Determination of Modal Density and Truncation: 

A method, called spectrum slicing, for determining the number of eigenvalues 
in a system, below a certain value, was presented in section 2.6. Equation (2.48) 
from that section, has been rewritten below in a form consistent with the equations in 
Tables (4.1) and (4.2). The triangular factorization of the matrices k u and m u , shifted 
to a cutoff value a c , is calculated 

n u )‘LDL T < 7 - 3 > 

where a c is calculated from the user defined cutoff frequency, by 

o e -(2« lf c f < 7 - 4 > 

Equation (2.49) can then rewritten as 

v(A -a c I) = v(k a -a c m a )=v(D) (7-5) 

where v is number of negative eigenvalues and A = diag (X,, Xi , ... , X n ). The 
number of negative elements of D is equal to the number of eigenvalues, n c , of the 
matrices k u and m u , which are less than C c . 

The Krylov sequence is continued until the total number of calculated static 
Ritz vectors, here defined as n r is greater than the number of eigenvalues below o c . 



135 


n c , multiplied by the fraction \|/. The fraction \\f was established by trial and error. 
Writing the termination logic in the form of the sub-indices from the equations of 
Tables (4.1) and (4.2) yields 

ifn r > ^fn c , thenj = l (7.6) 

All calculated static Ritz vectors, n„ are retained for subsequent dynamic analysis. 
The boundary flexibility algorithm, with targeted shifting and modal density truncation 
included, and numerical examples will be presented in the next chapter. 
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Mode 

No. 

Freq. 

(Hz) 

T1 

Effective Masses 
T2 T3 

(%) 

R1 

R2 

R3 

1 

0.55 

0.0 

0.0 

61.3 

97.0 

0.0 

0.0 

2 

0.55 

61.3 

0.0 

0.0 

0.0 

0.0 

97.0 

3 

3.41 

18.1 

0.0 

0.8 

0.1 

0.0 

2.4 

4 

3.41 

0.8 

0.0 

18.1 

2.4 

0.0 

0.1 

5 

9.52 

2.8 

0.0 

3.7 

0.2 

0.0 

0.1 

6 

9.52 

3.7 

0.0 

2.8 

0.1 

0.0 

0.2 

7 

18.57 

3.2 

0.0 

0.1 

0.0 

0.0 

0.1 

8 

18.57 

0.1 

0.0 

3.2 

0.1 

0.0 

0.0 

9 

27.58 

0.0 

80.6 

0.0 

0.0 

0.0 

0.0 

10 

30.53 

2.0 

0.0 

0.0 

0.0 

0.0 

0.0 

11 

30.53 

0.0 

0.0 

2.0 

0.0 

0.0 

0.0 

12 

45.41 

1.3 

0.0 

0.1 

0.0 

0.0 

0.0 

13 

45.41 

0.1 

0.0 

1.3 

0.0 

0.0 

0.0 

14 

63.18 

0.7 

0.0 

0.3 

0.0 

0.0 

0.0 

15 

63.18 

0.3 

0.0 

0.7 

0.0 

0.0 

0.0 

16 

81.90 

0.0 

8.6 

0.0 

0.0 

0.0 

0.0 

17 

83.59 

0.4 

0.0 

0.3 

0.0 

0.0 

0.0 

18 

83.59 

0.3 

0.0 

0.4 

0.0 

0.0 

0.0 

19 

84.64 

0.0 

0.0 

0.0 

0.0 

76.7 

0.0 

20 

104.41 

0.3 

0.0 

0.0 

0.0 

0.0 

0.0 

21 

104.41 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

22 

133.73 

0.0 

2.8 

0.0 

0.0 

0.0 

0.0 


T1 

Total Effective Mass (%) 
T2 T3 R1 

R2 

R3 

95.4 

92.1 95.4 100.0 

76.7 

100.0 


Table (7.1) - Beam Finite Element Model Effective Masses 
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Numerical Examples of Targeted Shifting and Modal Density Truncation 

8.1) Algorithm with Targeted Shifting and Modal Density Truncation: 

The complete static Ritz vector, boundary flexibility algorithm, with targeted 
shifting and modal density truncation included, is shown in Table (8.1). This 
algorithm was extended from that presented in chapter 4, which included block 
filtering and Cholesky/QR orthonormalization. Only the revised boundary flexibility 
algorithm for euclidean vector orthonormalization is presented in this chapter, because 
there is no fundamental difference in the application of targeted shifting and modal 
density truncation to the various orthonormalization options presented in chapter 4. 

8.2) Numerical Results and the Determination of the Fraction 

Three sets of time response problems were performed for the three example 
models described in chapter 5 using the applied loading described in chapter 6. These 
direct transient response problems were all performed using dynamically reduced 
models, a structural damping ration of 2%, with the beam example and the EPS 
Radiator having a fixed interface. The Cassini model was excited using enforced 
acceleration on a seismic mass. Two comparison cases for each model were generated 
using eigenvectors, the first using a modal truncation frequency greater than the 
dynamic range of interest, and the second using a frequency cutoff at the frequency 
range of interest. The beam example cutoff was 250 HZ for the high frequencies and 
100 HZ for the representative frequencies. The EPS Radiator and the Cassini cutoff 
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frequencies were 150 HZ for the high frequencies and 60 HZ for the representative 
frequencies. Three time response problems, using three different values for the fraction 
\j/ = 1.0, .5 and .25, for each example model created with static Ritz vectors, were 
performed. 

A selected set of physical accelerations and loads for each of these models was 
recovered. Table (8.2) shows the ratios of the minimum and maximum peak 
responses of the beam example, formulated using the options discussed above, to the 
prediction using the high frequency cutoff eigenvalue model. Tables (8.3) and (8.4) 
present similar tables for the EPS Radiator and the Cassini model transient solutions. 
There is no guarantee that the high frequency cutoff eigenvalue model represents a 
completely converged solution. The cutoff values were selected by multiplying the 
frequency at the cutoff of the dynamic range of interest by 2.5, and, considering the 
large size of the Cassini and the EPS Radiator models, this cutoff produced a size 
dynamic model for which a transient solution could still conveniently be obtained. 

Several observations can be made of the beam example ratios presented in 
Table (8.2). First, in the case of the axial tip acceleration, which is dependent on high 
frequency beam dynamics, the \|/ = 1.0 Ritz vector response predicted a more 
complete solution than did the eigenvector response. Figures (8.1) and (8.2) show 
axial tip acceleration time history plots predicted by the 250 Hz cutoff eigenvector and 
the \y = 1.0 Ritz vector models. The ratios for the vj/ = .50 and the \|/ = .25 beam 
examples are identical. This is because the truncation criteria algorithm produced 
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identical representations. There are so few eigenvectors in the frequency range of 
interest that, after the initial two vector blocks are calculated, there are already enough 
vectors to satisfy the truncation criteria. As a result, no more vectors were calculated 
for the \|/ = .50 beam example than were calculated for the \|t = .25 beam example, 
and targeted shifting and the representative static load is not utilized for either of these 
two cases. Not utilizing targeted shifting and a representative static load vector 
partially explain the relatively low response predictions of the vy = .50 and the \j/ = 
.25 cases. The lateral tip acceleration time history for the \|f = 1.0 case is illustrated 
in Figure (8.3), and the base bending moment for the same case is illustrated in Figure 
(8.4). The time histories of the other cases are not included because they are 
essentially identical. 

Table (8.3) presents the time response ratios generated using the EPS Radiator 
model. All the results from the various cases are reasonably consistent with the 150 
Hz frequency cutoff case. There begins to be some divergence from the 150 Hz case 
in the responses generated using the \|/ = .25 model. This would suggest that a 
fraction value of V|/ = .50 might be appropriate for accurate, but low computational 
cost, response predictions. In general, the responses predicted using the y = 1.0 Ritz 
vector component is a marginally closer match to the high frequency cutoff case than 
the 60 Hz eigenvector model. 

Figures (8.5), (8.6), (8.7), and (8.8) present the time response of the X 
acceleration, the Z acceleration, and the two bending moments, for the \\i = 1.0 case. 
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As in the beam example, the time histories plots of the EPS Radiator from the various 
cases are virtually identical. Examining the X acceleration time history explains the 
somewhat lower ratios obtained for the \|/ = .50 and the \\f = .25 maximum values 
comparisons. The minimum value is much greater than the maximum value, which 
represents a relatively small overshoot, and the ratio of the maximums represent a 
comparison of relatively small values. 

Table (8.4) presents the ratios of the time responses of the Cassini model using 
various y’s, compared to the 150 Hz cutoff eigenvector component. There is more 
divergence in the results of the various cases using the Cassini model than was in the 
previously discussed models. It was more difficult for the static Ritz vectors to match 
appendage accelerations than other responses. Once again, the results suggest that 
fraction value of approximately Vj / = .50 might be the appropriate choice for 
reasonably accurate response predictions. The High Gain Antenna Strut response 
predictions of the 60 Hz cutoff eigenvector component indicate that, for this element, 
the 60 Hz cutoff was too low. The Ritz vector models, whose vector size is based 
upon the number of vectors in the 60 Hz model, were also unable to accurately predict 
these strut loads. However, for the element forces in the main body stringer, the Ritz 
vectors were able to obtain more accurate response predictions than the 60 Hz 
eigenvector component. In conclusion, the use of a fraction value of vj/ = .50 may 
not provide enough Ritz vectors for accurate response predictions if the frequency 
cutoff of interest selected excludes important dynamic effects. 
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Figures (8.9), (8.10), (8.1 1), and (8.12) present the time response of the X, Y, 
and Z acceleration, and the strut axial load, for the \j/ = 1.0 case. The ratios of the 
maximum values of the accelerations in the Z direction were not presented in Table 
(8.4). An examination of Figure (8.11), the Z direction acceleration, will show that 
the maximum value occurs near the initial time step and it’s value is near zero. Ratios 
of these small response values was not meaningful and were not calculated. Appendix 
B contains the minimum and maximum responses of all the models and cases, and it 
was from this data that Tables (8.2), (8.3), and (8.4) were generated. 
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After the assembly of the component mass, m, and stiffness, fc, matrices: 
1) Modal Density Calculation 


° c = ( 2*/ c ) 2 

{k a -o c m a )=L c D e Lj = L c U c 

o c = o c + .l 

v(D) = n c 


[shift value is calculated from user 
defined cutoff frequency, f c ] 

[decompose shifted matrices] 

[if and only if (k u - cf is singular] 
[count the eigenvalues below a c ] 


2) Initialization 

<j> = -k r 1 k [create boundary flexibility matrix, kjj 

and ki C were defined in equation (2.3)] 

n =k ~ l (m +m I [create the initial block] 

HI v ii ic ic ' 


[For a free interface component, (If the component has rigid body modes then 
the elastic flexibility matrix, g e , defined by (2.35) through (2.39) is used): 


S a 


See 

Sic 


where g = fc' 1 


[boundary flexibility matrix] 


ql = smg a 


[create the initial block]] 


3a) Filtering of Initial Block 

=9 1 *[«fiag(?i* r «i*) > ec ]' 1/2 

^= 9 r * v * 


[normalization] 
[cross orthogonality] 


Table (8.1) - Revised Boundary Flexibility Algorithm 
With Targeted Shifting and Modal Density Truncation 
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3a) (Continued) Filtering of Initial Block 


\ 0 

L E ~ - 

L e \J 

[partition into lower triangular] 

If (< ||Z £ L> c ) * e, than dependent 

[infinity norm of each vector] 

* * r * * * 

01 "*[01 0DEP. 

[partition out dependent vectors] 

3b) Orthonormalization of the Initial Block 

r* ~ *** T _ *** 

l e = 01 01 

[B t B matrix product] 

l\=r t r 

[Cholesky factor decomposition] 

j t *** 

* 0i =0i 

[Solve by forward substitution] 


4) For Blocks j = 2,3, .... /, until n r > \\f n c 
4a) Set Shift Variables, Based on Count of n r 

If n r < ^ then k = 0 , o t = 0 [ n0 shift] 

If n > — — then k = 1 , a k = a T1 , a k = ^o c [ use °r;, ^ defined] 
r 4 3 

If n r > — — then k = 2 , o t = °* = ~°c [use cr r2 , if defined] 

If n r > ^-^ £ , then k = 3 , a k = a c [reuse L f and £7 C from modal 

density calculation] 

Table (8.1) (Continued) - Revised Boundary Flexibility Algorithm 
With Targeted Shifting and Modal Density Truncation 
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4b) Vector Calculation and Filtering 



[if k = 3, then L s = L c and U s = U c ] 

[if and only if (k u - Gjn u ) is singular] 

[create new vector block, unless k = 1] 
[create new vector block, if k = 1 ] 
[normalization] 

[cross orthogonality] 

[partition into lower triangular] 

[infinity norm of each vector] 

[partition out dependent vectors] 

[B t B matrix product] 

[Cholesky factor decomposition] 

[solve by forward substitution] 

[count total calculated vectors] 


5) When n r > \\i n c , Transform System to Form Component 


Table (8.1) (Continued) - Revised Boundary Flexibility Algorithm 
With Targeted Shifting and Modal Density Truncation 
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Eigen - 100 Hz 

Ritz - v 

= 1.0 

Ritz - \ 

= .50 

Ritz - i 

/= .25 

Min 

Max 

Min 

Max 

Min 

Max 

Min 

Max 

| Tip Acceleration 

Lat 

1.00 

1.00 

1.00 

1.00 

.88 

.85 

.88 

.85 

Axial 

0.87 

0.89 

1.11 

1.13 

.85 

.85 

.85 

.85 

| Element Forces at the Base of the Beam 

M 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

V 

1.00 

1.00 

1.00 

1.00 

.99 

1.00 

.99 

1.00 


Table (8.2) - Ratios of Responses of the Beam Example Model Using Various 
Representations to the 250 HZ Cutoff Eigenvector Representation Responses 



Eigen - 60 Hz 

Ritz - 1 

= 1.0 

Ritz - \ 

= .50 

Ritz - \j/ = .25 

Min 

Max 

Min 

Max 

Min 

Max 

Min 

Max 1 

| Tip Acceleration jj 

X 

1.00 

1.01 

1.00 

.98 

.97 

.88 

.97 

.87 

Y 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

Z 

1.00 

1.01 

.99 

.99 

.98 

1.02 

.97 

1.00 

1 Element Forces at the Base of the Scissors Beam jj 

Ml 

1.02 

1.01 

1.01 

1.01 

1.00 

1.00 

1.01 

1.01 

M2 

1.02 

1.02 

1.00 

-J 

1.01 

1.02 

1.04 

1.07 

1.11 

V 

1.02 

1.02 

1.01 

1.00 

1.04 

1.02 

1.11 

1.07 

P 

1.02 

1.00 

1.00 

1.00 

1.02 

1.00 

1.12 

1.10 

T 

1.02 

1.01 

1.01 

1.01 

1.00 

1.00 

1.01 

1.01 


Table (8.3) - Ratios of Responses of the EPS Radiator Model Using Various 
Representations to the 150 HZ Cutoff Eigenvector Representation Responses 
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Eigen - 60 Hz 

Ritz - \ 

f= 1.0 

Ritz - \ 

r = .50 

Ritz - \ 

= .25 

Min 

Max 

Min 

Max 

Min 

Max 

Min 

Max 

| Acceleration at the Oxidizer Tank eg 

X 

1.00 

1.00 

1.07 

1.03 

1.06 

1.20 

.77 

.77 

Y 

1.00 

1.00 

.97 

1.02 

.97 

.99 

.96 

1.09 

Z 

1.00 


1.02 


.93 


.92 


Acceleration at the Probe eg 

X 

1.00 

1.00 

1.00 

.98 

1.16 

1.07 

1.03 

.96 

Y 

1.00 

1.00 

.98 

1.00 

.97 

1.01 

.97 

.92 

Z 

1.00 


1.00 


1.02 


1.01 


High Gain Antenna Strut Beam Element Forces || 

Ml 

1.00 

1.00 

1.00 

.96 

1.08 

.93 

.87 

.94 

M2 

.98 

1.02 

1.00 

.95 

.92 

.88 

.84 

.89 

P 

1.05 

.90 

.98 

.94 

.89 

.82 

.90 

.84 

|| Main Body Stringer Beam Element Forces |j 

Ml 

1.06 

.94 

.99 

1.01 

.99 

.99 

.97 

1.02 

M2 

.99 

1.00 

.99 

1.03 

.93 

1.00 

.99 

.99 

P 

1.02 

.99 

.99 

1.01 

1.00 

1.01 

.96 

i 

1.03 

|| Lower Equipment Module Skin Membrane Forces j] 

Fx 

1.00 

1.00 

1.00 

1.00 

1.01 

1.04 

.99 

.99 

Fy 

.99 

1.01 

1.00 

1.00 

.99 

1.00 

1.00 

1.00 

Fxy 

.99 

1.00 

.99 

1.00 

1.00 

1.01 

.92 

1.07 


Table (8.4) - Ratios of Responses of the Cassini Model Using Various 
Representations to the 150 HZ Cutoff Eigenvector Representation Responses 
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TIME I SECS) 


FIXED INTERFACE BEAM - VERSION 3.3 1 
LANCZOS EIGENVECTORS, 250 HZ 


SUPTOELQJENT 10 
SUBCASE 10 


Figure (8.1) - Beam Example Tip Axial Acceleration 
Generated Using a 250 Hz Eigenvector Cutoff 
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Figure (8.2) - Beam Example Tip Axial Acceleration 
Generated Using a vj/ = 1.0 Ritz Vector Representation 
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Figure (8.3) - Beam Example Tip Lateral Acceleration 
Generated Using a \|/ = 1.0 Ritz Vector Representation 
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FIXED INTERFACE BEAM - VERSION 3.3 1 
FRBO = 100HZ, MULTIPLE^ 1 . 0 , TAKGETl=.S46 
SOLUTION 69, SAMPLE TRANSIENT 


SUBCASE 10 


10 


Figure (8.4) - Beam Example Base Bending Moment 
Generated Using a \j/ = 1.0 Ritz Vector Representation 
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Figure (8.5) - EPS Radiator X Tip Acceleration 
Generated Using a \j / = 1.0 Ritz Vector Representation 
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Figure (8.6) - EPS Radiator Z Tip Acceleration 
Generated Using a \\f = 1.0 Ritz Vector Representation 
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EPS RADIATOR - VERSION 3.3 1 
FREQ = 60 HZ , MULTIPLES., TARGET 1 = 5. 
SOLUTION 69 , SAMPLE TRANSIENT 


SUPERELEMENT 422 
SUBCASE 10 


Figure (8.7) - EPS Radiator Scissor Beam Bending Moment 1 
Generated Using a \j/ = 1.0 Ritz Vector Representation 
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Figure (8.8) - EPS Radiator Scissor Beam Bending Moment 2 
Generated Using a \|/ = 1.0 Ritz Vector Representation 
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Figure (8.11) - Cassini Oxidizer Tank Z Acceleration 
Generated Using a \\f = 1.0 Ritz Vector Representation 
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CASSINI MODEL TEST - VERSION 3.3 1 
ORTHOPT = 1, ORTHNFM = 1 FREQ - 6 0 HZ 
MULTIPLE = i . , TARGET 1 * 5 HZ, TARGET2 = 13 HZ 


SUPERELUQ3TT 100 
SUBCASE 100 


Figure (8.12) - High Gain Antenna Strut Axial Force 
Generated Using a \ji = 1.0 Ritz Vector Representation 



Chapter 9 


Summary and Conclusions 

Using the boundary flexibility method to initialize the block-Krylov recurrence 
algorithm provides an efficient and simple method for generating static Ritz vectors. 
Static Ritz vectors so generated accurately can represent the dynamics of a 
substructure. Because this methodology does not require the solution of the 
component eigenvalue problem, a component can be formed with a significant 
decrease in computational cost. The component formed using static Ritz vectors can 
include fewer vectors than the comparable eigensolution, for similar accuracy, and the 
computational cost of the transient solution is then also reduced. 

This dissertation presents new developments in several areas related to static 
Ritz vector calculation using the boundary flexibility method. It has been shown that 
the loss of orthogonality, discussed in the literature, is directly related to convergence 
to an eigenvector in a power extraction method. Orthonormalization using the 
euclidean norm rather than the mass matrix has been demonstrated. The replacement 
of Gram-Schmidt with Cholesky/QR orthonormalization has also been demonstrated. 
These two modifications to the orthonormalization algorithm were developed to 
decrease the primary computational cost of the block-Krylov sequence. 

The Krylov blocks produced by the boundary flexibility method are initially 
to large for efficient handling and computation. As a result, a method of block 
filtering has been developed which retains the physically significant information 
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contained in the vector block, while eliminating the redundant information. Block 
filtering reduces the block size to that commonly used in commercial Lanczos 
eigenvalue extraction routines. 

The use of inverse operators, or shifting, commonly used in Lanczos 
eigenvalue extraction routines to search for missing eigenvectors, has been applied to 
the boundary flexibility method. Shifting alters the Krylov sequence so that vectors 
near a selected frequency are created. Targeted shifting utilizes the frequency 
spectrum of the applied dynamic loading to select a targeted frequency. The spacial 
distribution of the applied dynamic loading can be used in conjunction with targeted 
shifting to further refine the static Ritz vector creation. 

Determining when a sufficient number of static Ritz vectors have been 
calculated to accurately represent a component has been difficult. Truncation of the 
Krylov sequence based upon the modal density of the given finite element model has 
been developed and presented. This heuristic method is based on the observation that 
dynamically significant static Ritz vectors are calculated early in the Krylov sequence, 
and that fewer static Ritz vectors than eigenvectors are necessary to accurately 
represent a component. 

Potential future work on the boundary flexibility method of static Ritz vector 
creation could include the use of synthetic load vectors to supplement the spatial 
distribution of the applied dynamic load. A synthetic load could be created which 
would cause a particular element, or sets of elements, to deflect. Applying this load 
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to the Krylov sequence would create static Ritz vectors which should guarantee 
accurate response predictions for specific, critical data recovery items. This technique 
could be investigated with no additional theoretical development. 

The shift strategy developed, and the value of the fraction y, could be tested 
with a larger number of example models. This might allow for a either a 
simplification, or a greater sophistication of the shift strategy. In particular, an 
alteration of the shift strategy, to allow the use of the applied load and targeted 
shifting, on models with a number of eigenvectors in the frequency range of interest 
smaller than the block size, should be investigated. 
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Appendix A1 - EPS Radiator Finite Element Model Effective Masses 


Mode 

Frequency 


Effective Masses (%) 



No. 

(Hz) 

T1 

T2 

T3 

R1 

R2 

R3 

1 

0.19 

0.0 

0.0 

46.8 

93.7 

35.8 

0.0 

2 

0.70 

0.1 

0.0 

0.0 

0.0 

11.2 

0.1 

3 

0.91 

43.4 

0.0 

0.0 

0.0 

1.0 

85.0 

4 

1.20 

0.0 

35.7 

0.1 

0.0 

0.1 

2.2 

5 

1.30 

0.0 

0.3 

12.3 

3.6 

9.4 

0.0 

6 

2.59 

0.4 

0.0 

0.0 

0.0 

1.6 

0.2 

7 

3.32 

19.2 

0.0 

0.0 

0.0 

0.6 

6.1 

8 

3.65 

0.0 

0.1 

3.7 

0.7 

2.8 

0.0 

9 

5.15 

0.0 

9.9 

0.0 

0.0 

0.0 

0.6 

10 

5.78 

0.2 

0.0 

0.0 

0.0 

0.6 

0.0 

11 

6.73 

0.0 

0.5 

1.7 

0.3 

1.3 

0.0 

12 

6.88 

4.3 

0.0 

0.0 

0.0 

0.2 

0.4 

13 

8.92 

0.0 

0.5 

0.4 

0.0 

0.3 

0.0 

14 

9.72 

0.9 

0.0 

0.0 

0.0 

0.0 

0.1 

15 

9.96 

0.0 

1.0 

1.2 

0.2 

0.9 

0.1 

16 

10.33 

0.8 

0.0 

0.0 

0.0 

0.5 

0.1 

17 

10.72 

0.0 

0.2 

5.3 

0.5 

4.0 

0.0 

18 

12.08 

1.1 

0.0 

0.0 

0.0 

0.1 

0.1 

19 

13.15 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

20 

14.51 

0.8 

0.0 

0.0 

0.0 

0.0 

0.1 

21 

15.05 

0.0 

0.0 

0.0 

0.0 

0.2 

0.0 

22 

16.34 

0.7 

0.0 

0.0 

0.0 

0.0 

0.0 

23 

16.34 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

24 

16.91 

0.0 

0.0 

0.0 

0.0 

0.1 

0.0 

25 

17.83 

0.7 

0.0 

0.0 

0.0 

0.0 

0.0 

26 

18.71 

0.7 

0.0 

0.0 

0.0 

0.0 

0.0 

27 

19.15 

0.0 

1.2 

0.1 

0.0 

0.1 

0.1 

28 

20.25 

0.0 

0.3 

0.0 

0.0 

0.0 

0.0 

29 

21.03 

0.1 

0.0 

0.0 

0.0 

1.2 

0.0 

30 

22.80 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

31 

22.94 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

32 

23.13 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

33 

24.03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

34 

24.04 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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35 

24.80 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

36 

25.52 

0.0 

0.5 

0.1 

0.0 

0.1 

0.0 

37 

25.64 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

38 

25.75 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

39 

26.23 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

40 

26.34 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

41 

27.01 

0.0 

1.1 

0.5 

0.0 

0.4 

0.1 

42 

27.43 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

43 

27.69 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

44 

28.19 

0.0 

0.0 

0.1 

0.0 

0.1 

0.0 

45 

28.46 

0.0 

0.3 

0.0 

0.0 

0.0 

0.0 

46 

28.91 

0.0 

1.7 

0.8 

0.1 

0.6 

0.1 

47 

29.61 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

48 

30.23 

0.0 

0.0 

0.5 

0.0 

0.4 

0.0 

49 

30.88 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

50 

30.89 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

51 

31.13 

0.0 

2.2 

0.0 

0.0 

0.0 

0.1 

52 

33.69 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

53 

33.79 

0.0 

0.4 

0.0 

0.0 

0.0 

0.0 

54 

34.35 

0.0 

2.1 

0.0 

0.0 

0.0 

0.1 

55 

34.71 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

56 

36.57 

0.0 

0.3 

0.0 

0.0 

0.0 

0.0 

57 

37.37 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

58 

37.51 

0.0 

1.6 

0.0 

0.0 

0.0 

0.1 

59 

37.53 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

60 

38.26 

0.0 

0.4 

0.0 

0.0 

0.0 

0.0 

61 

38.61 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

62 

39.79 

0.0 

0.0 

0.0 

0.0 

0.1 

0.1 

63 

39.83 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

64 

40.88 

0.0 

9.5 

0.0 

0.0 

0.0 

0.6 

65 

41.08 

0.0 

0.2 

0.0 

0.0 

0.0 

0.0 

66 

42.45 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

67 

44.29 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

68 

44.52 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

69 

45.73 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

70 

46.20 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

71 

46.53 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

72 

46.74 

0.0 

2.6 

0.1 

0.0 

0.0 

0.2 

73 

48.69 

0.0 

0.3 

0.0 

0.0 

0.0 

0.0 

74 

49.37 

0.0 

1.8 

0.0 

0.0 

0.0 

0.1 

75 

49.70 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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76 

50.80 

0.0 

0.2 

0.0 

0.0 

0.0 

0.0 

77 

51.21 

0.0 

0.6 

0.0 

0.0 

0.0 

0.0 

78 

51.81 

0.0 

0.0 

0.0 

0.0 

0.0 

0.2 

79 

52.04 

0.0 

0.5 

0.0 

0.0 

0.0 

0.0 

80 

53.69 

0.0 

1.7 

0.1 

0.0 

0.0 

0.1 

81 

55.50 

0.1 

0.0 

0.0 

0.0 

0.0 

0.1 

82 

56.10 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

83 

56.65 

0.0 

4.7 

0.0 

0.0 

0.0 

0.3 

84 

56.97 

0.1 

0.1 

0.0 

0.0 

0.0 

0.0 

85 

57.69 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

86 

58.28 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

87 

59.33 

0.0 

0.5 

0.0 

0.0 

0.0 

0.0 

88 

60.24 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

89 

61.50 

0.0 

0.6 

0.0 

0.0 

0.0 

0.0 

90 

62.01 

0.0 

0.8 

0.1 

0.0 

0.1 

0.0 

91 

62.45 

0.0 

0.9 

0.0 

0.0 

0.0 

0.1 

92 

63.12 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

93 

63.86 

0.0 

0.2 

0.0 

0.0 

0.0 

0.0 

94 

64.84 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

95 

67.15 

0.0 

0.5 

0.0 

0.0 

0.0 

0.0 

96 

70.52 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

97 

70.53 

0.0 

0.0 

0.4 

0.0 

0.3 

0.0 

98 

71.34 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

99 

73.06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

100 

74.40 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

101 

75.66 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

102 

76.92 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

103 

79.49 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

104 

79.89 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

105 

80.46 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

106 

82.72 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

107 

83.36 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

108 

84.53 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

109 

85.43 

0.0 

1.5 

0.0 

0.0 

0.0 

0.1 

110 

86.57 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

111 

87.00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

112 

88.86 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

113 

89.31 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

114 

89.90 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

115 

92.33 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

116 

92.75 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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117 

94.21 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

118 

96.19 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

119 

96.53 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

120 

98.07 

1.0 

0.1 

0.0 

0.0 

0.1 

0.1 

121 

98.77 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

122 

99.72 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

123 

100.06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

124 

100.19 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

125 

100.46 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

126 

101.95 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

127 

102.02 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

128 

102.69 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

129 

103.97 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

130 

104.88 

0.0 

0.1 

0.5 

0.0 

0.4 

0.0 

131 

106.50 

0.0 

0.2 

0.1 

0.0 

0.1 

0.0 

132 

107.69 

0.0 

0.7 

0.0 

0.0 

0.0 

0.1 

133 

107.93 

0.1 

0.6 

0.0 

0.0 

0.0 

0.0 

134 

108.61 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

135 

108.78 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

136 

109.35 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

137 

110.03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

138 

111.76 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

139 

111.93 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

140 

113.42 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

141 

113.43 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

142 

113.70 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

143 

113.70 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

144 

114.31 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

145 

114.41 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

146 

115.00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

147 

115.16 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

148 

116.00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

149 

117.72 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

150 

117.95 

0.2 

0.0 

0.0 

0.0 

0.0 

0.0 

151 

118.99 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

152 

121.38 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

153 

122.25 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

154 

123.42 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

155 

123.68 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

156 

124.66 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

157 

125.34 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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158 

125.37 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

159 

126.61 

0.0 

0.0 

0.0 

0.0 

.0.0 

0.0 

160 

127.27 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

161 

128.22 

1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

162 

128.68 

0.6 

0.0 

0.0 

0.0 

0.0 

0.0 

163 

129.60 

0.2 

0.0 

0.0 

0.0 

0.0 

0.0 

164 

131.13 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

165 

131.41 

0.2 

0.0 

0.0 

0.0 

0.0 

0.0 

166 

133.15 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

167 

133.57 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

168 

133.86 

0.2 

0.0 

0.0 

0.0 

0.0 

0.0 

169 

134.39 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

170 

135.99 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

171 

136.67 

0.0 

0.0 

0.3 

0.0 

0.2 

0.0 

172 

137.31 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

173 

137.41 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

174 

138.89 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

175 

139.59 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


Total Effective Mass (%) 

T2 T3 R1 R2 R3 

90.3 75.6 99.4 75.7 98.8 


T1 

78.1 



Appendix A2 - Cassini Finite Element Model Effective Masses 


Mode Frequency Effective Masses (%) 


No. 

(Hz) 

T1 

T2 

T3 

R1 

R2 

R3 

1 

7.36 

49.2 

8.6 

0.0 

13.2 

78.7 

0.0 

2 

7.70 

7.7 

50.9 

0.0 

78.8 

12.5 

0.0 

3 

12.19 

0.0 

0.0 

0.0 

0.0 

0.0 

0.2 

4 

14.89 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

14.91 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

6 

15.50 

0.0 

0.0 

0.0 

0.0 

0.0 

13.9 

7 

15.75 

0.0 

0.2 

0.0 

0.0 

0.0 

60.5 

8 

15.88 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9 

15.90 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

10 

17.85 

0.3 

0.0 

5.2 

0.0 

0.1 

0.0 

11 

18.28 

0.0 

0.7 

0.0 

0.1 

0.0 

1.0 

12 

19.01 

3.0 

0.2 

0.4 

0.0 

0.7 

2.2 

13 

19.08 

4.3 

0.0 

0.5 

0.0 

0.9 

0.8 

14 

19.21 

0.1 

0.9 

0.0 

0.2 

0.0 

0.0 

15 

19.54 

1.7 

0.1 

5.2 

0.0 

0.6 

0.2 

16 

19.81 

0.5 

0.0 

0.5 

0.0 

0.1 

0.4 

17 

19.96 

0.0 

0.2 

0.0 

0.1 

0.0 

0.0 

18 

20.04 

0.0 

1.7 

0.0 

0.6 

0.0 

0.5 

19 

20.34 

0.4 

2.0 

4.1 

0.3 

0.1 

0.9 

20 

20.47 

1.4 

1.8 

3.5 

0.2 

0.3 

1.2 

21 

20.63 

1.2 

0.3 

0.0 

0.0 

0.2 

0.1 

22 

21.39 

5.4 

0.0 

1.1 

0.0 

1.3 

0.0 

23 

23.44 

0.0 

0.0 

0.0 

0.0 

0.0 

0.7 

24 

23.73 

1.5 

0.0 

0.8 

0.0 

0.3 

0.0 

25 

25.55 

0.0 

9.0 

0.4 

2.5 

0.0 

3.4 

26 

25.94 

1.0 

0.0 

0.1 

0.0 

0.1 

0.0 

27 

26.54 

0.0 

0.6 

0.1 

0.2 

0.0 

2.5 

28 

27.14 

1.0 

0.0 

0.0 

0.0 

0.2 

0.0 

29 

27.23 

0.0 

0.0 

0.1 

0.0 

0.0 

0.1 

30 

27.45 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

31 

27.52 

0.0 

0.0 

0.2 

0.0 

0.0 

0.0 

32 

28.33 

1.7 

0.1 

14.2 

0.1 

0.5 

0.1 

33 

29.60 

0.3 

0.6 

2.8 

0.1 

0.1 

0.0 

34 

29.70 

0.5 

1.1 

2.3 

0.2 

0.1 

0.0 
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35 

31.27 

1.9 

2.9 

4.6 

0.5 

0.4 

0.1 

36 

31.51 

0.0 

0.1 

2.1 

0.0 

0.0 

0.0 

37 

33.53 

0.4 

0.0 

1.9 

0.0 

0.1 

0.0 

38 

34.49 

2.6 

0.0 

0.0 

0.0 

0.5 

0.0 

39 

34.99 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

40 

35.23 

0.3 

0.0 

0.3 

0.0 

0.1 

0.0 

41 

36.45 

1.2 

3.7 

0.5 

0.6 

0.2 

0.0 

42 

36.73 

2.7 

1.4 

0.1 

0.3 

0.5 

0.1 

43 

37.53 

0.8 

2.1 

15.6 

0.4 

0.2 

0.1 

44 

38.04 

0.0 

1.1 

1.6 

0.2 

0.0 

0.0 

45 

38.29 

0.0 

0.1 

19.6 

0.0 

0.0 

0.0 

46 

38.53 

0.1 

0.0 

0.1 

0.0 

0.0 

0.1 

47 

39.26 

0.2 

0.7 

0.0 

0.1 

0.0 

0.2 

48 

39.39 

1.2 

0.0 

0.5 

0.0 

0.3 

0.1 

49 

40.65 

0.0 

3.0 

0.0 

0.5 

0.0 

0.0 

50 

41.70 

0.4 

0.2 

0.2 

0.0 

0.1 

0.0 

51 

41.97 

0.2 

0.1 

0.4 

0.0 

0.0 

0.0 

52 

42.69 

0.2 

0.1 

0.1 

0.0 

0.0 

0.1 

53 

43.38 

1.6 

0.0 

0.1 

0.0 

0.2 

0.1 

54 

44.23 

0.1 

0.0 

0.5 

0.0 

0.0 

0.0 

55 

44.58 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

56 

44.98 

0.1 

0.0 

1.3 

0.0 

0.0 

0.0 

57 

45.11 

0.2 

0.0 

0.9 

0.0 

0.0 

0.0 

58 

46.03 

0.0 

1.0 

0.1 

0.2 

0.0 

0.0 

59 

47.12 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3 

60 

47.71 

0.0 

0.0 

0.2 

0.0 

0.0 

0.0 

61 

48.40 

0.1 

0.1 

1.3 

0.0 

0.0 

0.0 

62 

48.89 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

63 

49.13 

0.2 

0.0 

0.1 

0.0 

0.0 

0.7 

64 

50.11 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

65 

50.11 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

66 

50.12 

0.2 

0.0 

0.0 

0.0 

0.0 

0.0 

67 

50.93 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

68 

52.72 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

69 

53.08 

0.2 

0.1 

0.0 

0.0 

0.0 

0.0 

70 

53.31 

0.2 

0.0 

0.0 

0.0 

0.0 

0.0 

71 

54.10 

0.0 

0.2 

0.0 

0.0 

0.0 

0.0 

72 

55.12 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

73 

56.42 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

74 

57.31 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

75 

58.35 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 
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76 

58.75 

0.0 

0.1 

0.1 

0.0 

0.0 

0.0 

77 

59.10 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

78 

59.12 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

79 

59.12 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

80 

59.13 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

81 

59.30 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

82 

59.72 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

83 

59.93 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

84 

60.32 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

85 

61.40 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

86 

62.01 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

87 

62.77 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

88 

63.21 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

89 

64.52 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

90 

65.18 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

91 

67.86 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

92 

67.88 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

93 

67.88 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

94 

67.88 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

95 

67.91 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

96 

67.94 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

97 

68.75 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

98 

69.50 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

99 

69.86 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

100 

71.31 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

101 

72.05 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

102 

73.13 

0.0 

0.0 

0.8 

0.0 

0.0 

0.0 

103 

75.57 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

104 

75.57 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

105 

75.59 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

106 

75.59 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

107 

75.64 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

108 

75.68 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

109 

75.84 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

110 

76.07 

0.0 

0.0 

0.0 

0.0 

0.0 

0.2 

111 

76.43 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

112 

76.81 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

113 

78.81 

0.0 

0.0 

0.0 

0.0 

0.0 

0.5 

114 

79.64 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

115 

80.36 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3 

116 

82.57 

0.1 

0.0 

0.0 

0.0 

0.0 

0.3 
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117 

83.03 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

118 

83.48 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

119 

84.86 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

120 

85.93 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

121 

86.70 

0.0 

0.2 

0.0 

0.0 

0.0 

0.0 

122 

87.09 

0.0 

0.0 

0.6 

0.0 

0.0 

0.0 

123 

87.80 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

124 

88.35 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

125 

88.41 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

126 

88.81 

0.2 

0.0 

0.3 

0.0 

0.0 

0.0 

127 

90.55 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

128 

90.56 

0.1 

0.0 

0.1 

0.0 

0.0 

0.0 

129 

90.70 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

130 

91.26 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

131 

91.73 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

132 

92.47 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

133 

92.88 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

134 

93.82 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

135 

94.24 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

136 

94.37 

0.1 

0.1 

0.0 

0.0 

0.0 

0.0 

137 

94.99 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

138 

95.03 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

139 

96.00 

0.1 

0.0 

0.0 

0.0 

0.0 

0.1 

140 

96.25 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

141 

97.55 

0.0 

0.2 

0.1 

0.0 

0.0 

0.0 

142 

99.16 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

143 

100.37 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

144 

101.77 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

145 

101.93 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

146 

103.05 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

147 

104.03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

148 

104.08 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

149 

104.62 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

150 

106.66 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

151 

106.66 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

152 

106.69 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

153 

107.17 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

154 

107.17 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

155 

107.85 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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156 

109.30 

0.0 

0.0 

0.0 

0.0 

0.0 

0.5 

157 

109.53 

0.0 

0.0 

0.0 

0.0 

0.0 

0.2 




Total Effective Mass (%) 




T1 

T2 

T3 

R1 

R2 

R3 



97.8 

97.8 

97.8 

99.8 

99.8 

93.8 



Appendix B 




M 

I N / M A X 

S U M M A 

R Y 


STRUCTURE : 

beam 

PARAM2 

: eigen PARAM3 : 250 

HZ RESPONSE 

: ALL 


ID 


VALUE 

OCCURRENCE 

VALUE 

OCCURRENCE 

ACCE 

10 

3 

-2 . 2539 OE+Ol 

.9300 

1 . 79300E+01 

.0100 

ACCE 

10 

4 

-7.88126E-01 

.9200 

8 . 04331E- 01 

. 0000 

ACCE 

10 

5 

- 2 . 25390E+01 

.9300 

1.79300E+01 

. 0100 

EL FOR 

1 

4 

- 2 . 61165E+03 

1.8250 

2 . 54222E+03 

1.0250 

EL FOR 

1 

5 

-2 - 54222E+03 

1.0250 

2 . 61165E+03 

1.8250 

EL FOR 

1 

6 

-3 . 1983 8E + 00 

1.8200 

3 . 31442E+00 

1.0400 

EL FOR 

1 

7 

-3 . 31442E+00 

1.0400 

3 . 19838E+00 

1.8200 

EL FOR 

1 

8 

- 2 . 29287E- 01 

. 9450 

1.27743E+00 

. 0250 

EL FOR 

1 

9 

-4 . 65317E- 10 

.9450 

2 . 67230E- 08 

. 0250 
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MIN/MAX SUMMARY 

STRUCTURE: beam PARAM2 : eigen PARAM3 : 100 HZ RESPONSE: ALL 

ID VALUE OCCURRENCE VALUE OCCURRENCE 


ACCE 

10 

3 

ACCE 

10 

4 

ACCE 

10 

5 

EL FOR 

1 

4 

EL FOR 

1 

5 

EL FOR 

1 

6 

EL FOR 

1 

7 

EL FOR 

1 

8 

EL FOR 

1 

9 


-2 . 24629E+01 .9300 
-6 . 83512E-01 .9200 
-2 . 24629E+01 .9300 
-2.61165E+03 1.8250 
-2 . 54222E+03 1.0250 
-3 . 19838E+00 1.8200 
-3 . 31442E+00 1.0400 
-2.28109E-01 .9450 
-2.13987E-11 .1400 


1 . 78509E+01 .0100 
6 . 97566E-01 .0000 
1 . 78509E+01 .0100 
2 . 54222E+03 1.0250 
2.61165E+03 1.8250 
3 . 31442E+00 1.0400 
3 . 19838E+00 1.8200 
1 . 08525E+00 .0250 
1 . 09748E- 13 1.0800 
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MIN /MAX SUMMARY 


STRUCTURE : 

beam 

ID 

P ARAM 2 

: ritz PARAM3: fracl 

VALUE OCCURRENCE 

. 0 RESPONSE 

VALUE 

: ALL 

OCCURRENCE 

ACCE 

10 

3 

-2 . 24504E+01 

.9300 

1.78359E+01 

.0100 

ACCE 

10 

4 

- 8 . 72 979E- 01 

. 9200 

8 . 90928E- 01 

.0000 

ACCE 

10 

5 

-2.24 5 04E+01 

.9300 

1 . 78359E+0 1 

.0100 

EL FOR 

1 

4 

-2 . 61165E+03 

1.8250 

2 . 54221E+03 

1.0250 

EL FOR 

1 

5 

-2 . 54221E+03 

1.0250 

2 . 61165E+03 

1.8250 

EL FOR 

1 

6 

-3 . 19838E+0 0 

1.8200 

3 . 31401E+00 

1.0400 

EL FOR 

1 

7 

-3 . 314 01E+0 0 

1.0400 

3 . 19838E+00 

1.8200 

EL FOR 

1 

8 

-2 . 17777E- 01 

.9450 

1.12482E+00 

.0250 

EL FOR 

1 

9 

0 .OOOOOE+OO 

.0000 

0 . OOOOOE+OO 

. 0000 
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STRUCTURE: beam 
ID 


ACCE 

10 

ACCE 

10 

ACCE 

10 

EL FOR 

1 

EL FOR 

1 

EL FOR 

1 

EL FOR 

1 

EL FOR 

1 

EL FOR 

1 


MIN /MAX SUMMARY 


PARAM2 : ritz 

VALUE 


PARAM3: frac.5 
OCCURRENCE 


RESPONSE: ALL 
VALUE OCCURRENCE 


3 

-1 . 98163E+01 

.9300 

4 

- 6 . 70055E- 01 

.0200 

5 

-1 . 98163E+01 

.9300 

4 

-2 . 61174E+03 

1.8250 

5 

-2 . 53853E+03 

1.0250 

6 

- 3 . 20066E+00 

1.8200 

7 

-3 . 27943E+Q0 

1.0450 

8 

- 2 . 27150E- 01 

.9450 

9 

0 . 000 00E+00 

. 0000 


1.52155E+01 

. 0100 

6.68161E-01 

.9400 

1 .52155E+01 

. 0100 

2 . 53853E+03 

1.0250 

2 . 61174E+03 

1.8250 

3 . 27943E+00 

1.0450 

3 . 20 066E+00 

1.8200 

1 . 01833E+Q0 

.0250 

0 .00000E+00 

.0000 
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MIN /MAX SUMMARY 


STRUCTURE : 

beam 

ID 

PARAM2 

: ritz ] 

VALUE 

PARAM3 : f rac . 
OCCURRENCE 

.25 RESPONSE 

VALUE 

: ALL 

OCCURRENCE 

ACCE 

10 

3 

-1.98163E+01 

.9300 

1 . 52155E+01 

.0100 

ACCE 

10 

4 

-6 . 70 055E- 01 

. 0200 

6 . 68161E-01 

. 9400 

ACCE 

10 

5 

-1 . 98 163E+01 

. 9300 

1 . 52155E+0 1 

. 0100 

EL FOR 

1 

4 

-2 . 61174E+03 

1.8250 

2 . 53853E+03 

1.0250 

EL FOR 

1 

5 

- 2 . 53 8 53E+03 

1.0250 

2 . 61174E+03 

1 . 8250 

EL FOR 

1 

6 

- 3 . 20066E+00 

1.8200 

3 . 27943E+00 

1.0450 

EL FOR 

1 

7 

-3 . 27943E+00 

1.0450 

3 . 20066E+00 

1 . 8200 

EL FOR 

1 

8 

-2 . 27150E- 01 

.9450 

1 . 01833E+00 

.0250 

EL FOR 

1 

9 

0 . 000 0 OE + 00 

.0000 

0 . 00000E+00 

. 0000 
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MIN/MAX SUMMARY 


STRUCTURE: radiator PARAM2 : eigen 
ID VALUE 


PARAM3: 150 HZ RESPONSE: ALL 
OCCURRENCE VALUE OCCURRENCE 


ACCE 440827 3 
ACCE 440827 4 
ACCE 440827 5 
EL FOR 440155 4 
EL FOR 440155 5 
EL FOR 440155 6 
EL FOR 440155 7 
EL FOR 440155 8 
EL FOR 440155 9 


-9 . 96496E-01 31.0078 
-5 . 71353E+00 .1150 
-1 . 67945E+00 17.0102 
-2 . 91732E-01 21.6194 
-1 . 54566E+00 16.2353 
-5 . 46363E-02 18.1800 
-4 . 14298E-01 21.6094 
-9 . 38668E-01 16.2503 
-2 . 75216E- 01 21.6194 


2 . 54486E-01 17.0102 
1 . 71285E+00 17.1851 
2 . 06374E+00 21.3244 
9.41954E-01 .1250 
6 . 89794E- 01 21.6094 
7 . 09203E- 02 16.9851 
9 . 28340E- 01 16.2353 
6 . 09971E-01 18.7948 
9 . 15571E- 01 17.0051 
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MIN/MAX SUMMARY 

STRUCTURE: radiator PARAM2 : eigen PARAM3 : 60 HZ RESPONSE: ALL 

ID VALUE OCCURRENCE VALUE OCCURRENCE 


ACCE 

440827 

3 

ACCE 

440827 

4 

ACCE 

440827 

5 

EL 

FOR 

440155 

4 

EL 

FOR 

440155 

5 

EL 

FOR 

440155 

6 

EL 

FOR 

440155 

7 

EL 

FOR 

440155 

8 

EL 

FOR 

440155 

9 


-I . 00129E+00 31.0078 
-5 . 71337E+00 .1150 
- 1 . 6 9103E+00 17.0102 
-2.96674E-01 21.6194 
-1 . 57609E+00 16.2353 
- 5 . 56583E-02 18.1800 
-4 . 24259E-01 21.6094 
-9 . 60995E-01 16.2503 
-2 . 79499E-01 21.6194 


2 . 57273E-01 17.0102 
1 . 71283E+00 17.1851 
2 . 07328E+00 21.3244 
9 . 55594E-01 .1250 
7 . 06379E- 01 21.6094 
7 . 91698E- 02 16.9902 
9 . 46616E- 01 16.2353 
6 . 08715E- 01 18.7948 
9 . 28780E- 01 17.0051 
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MIN /MAX SUMMARY 

STRUCTURE: radiator PARAM2 : ritz PARAM3 : fracl.O RESPONSE: ALL 

ID VALUE OCCURRENCE VALUE OCCURRENCE 


ACCE 440827 3 
ACCE 440827 4 
ACCE 440827 5 
EL FOR 440155 4 
EL FOR 440155 5 
EL FOR 440155 6 
EL FOR 440155 7 
EL FOR 440155 8 
EL FOR 440155 9 


- 9 . 924 73E-01 31.0078 
-5 . 71352E+00 .1150 
-1 . 65486E+00 17.0102 
-2 . 93704E-01 21.6194 
- 1 . 54285E+00 16.2353 
-5 . 56114E-02 18.1800 
-4 . 16776E-01 21.6094 
-9 . 40207E-01 16.2503 
-2 . 76952E- 01 21.6194 


2 . 50247E-01 17.0102 
1 . 71291E+00 17.1851 
2 . 03754E+00 21.3244 
9 . 49249E-01 .1250 
6 . 93919E-01 21.6094 
7 . 54589E-02 16.9902 
9 . 26655E- 01 16.2353 
6 . 07573E-01 18.7948 
9 . 22073E-01 17.0051 
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MIN/MAX SUMMARY 

STRUCTURE : radiator P ARAM 2 : ritz P ARAM 3 : f rac . 5 RESPONSE : ALL 

ID VALUE OCCURRENCE VALUE OCCURRENCE 


ACCE 

440827 

3 

ACCE 

440827 

4 

ACCE 

440827 

5 

EL 

FOR 

440155 

4 

EL 

FOR 

440155 

5 

EL 

FOR 

440155 

6 

EL 

FOR 

440155 

7 

EL 

FOR 

440155 

8 

EL 

FOR 

440155 

9 


- 9 . 684 25E-01 31.0078 
-5 . 72016E+00 .1150 
-1 . 64152E+00 17.0102 
-2 . 9202 IE-01 21.6144 
- 1 . 57293E+00 16.2353 
-5 . 54656E-02 18.1800 
-4 . 30009E-01 21.6094 
-9 . 55818E- 01 16.2553 
-2.76342E-01 21.6194 


2 . 24475E- 01 17.0501 
1 . 71146E+00 17.1851 
2 . 09572E+00 21.3244 
9 . 44502E- 01 .1250 
7 . 15953E- 01 21.6094 
7 . 52191E- 02 16.9902 
9 . 44717E- 01 16.2353 
6 . 10545E- 01 18.7948 
9 . 15235E- 01 17.0051 
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MIN/MAX SUMMARY 


STRUCTURE: radiator PARAM2 : ritz 

ID VALUE 


PARAM3: frac.25 RESPONSE: ALL 
OCCURRENCE VALUE OCCURRENCE 


ACCE 440827 3 
ACCE 440827 4 
ACCE 440827 5 
EL FOR 440155 4 
EL FOR 440155 5 
EL FOR 440155 6 
EL FOR 440155 7 
EL FOR 440155 8 
EL FOR 440155 9 


-9 . 65517E- 01 31.0078 
-5 . 723 86E+00 .1150 
-1 . 62194E+00 17.0102 
-2 . 95189E-01 21.6144 
-1 . 65746E+00 16.2453 
-5 . 35842E- 02 18.1800 
-4 . 61048E-01 21.6094 
-1 . 04982E+00 16.2553 
-2 . 78004E-01 21.6194 


2 . 21647E- 01 17.0551 
1 . 71038E+00 17.1851 
2 . 05417E+00 21.3244 
9 . 55356E-01 .1250 
7 . 67631E- 01 21.6094 
7 . 54796E- 02 16.9851 
9 . 95492E-01 16.2453 
6 . 70256E-01 17.9000 
9 . 23280E-01 17.0051 
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MIN/MAX SUMMARY 


STRUCTURE : 

cassini 

ID 

PARAM2 

: eigen P ARAMS : 150 

VALUE OCCURRENCE 

HZ RESPONSE 

VALUE 

: ALL 

OCCURRENCE 

ACCE 

10004 

3 

-5 . 88261E+02 

1.3670 

4 . 88726E+02 

.7680 

ACCE 

10004 

4 

-4 . 04763E+02 

1 . 9490 

5 . 37344E+02 

.6300 

ACCE 

10004 

5 

-1 .75641E+03 

.5140 

3 . 23883E+02 

.6850 

ACCE 

701 

3 

-3 . 12625E+02 

1.2130 

3 . 79308E+02 

1.1050 

ACCE 

701 

4 

- 2 . 75651E+02 

1.9610 

5. 07647E+02 

.6280 

ACCE 

701 

5 

-1.14259E+03 

.7740 

-6.93264E+00 

.0000 

EL FOR 

11601 

4 

-1 . 04208E+02 

.5960 

5.69385E+01 

.8620 

EL FOR 

11601 

5 

-3 . 49150E+01 

1 . 1870 

4 . 77298E+01 

1.1300 

EL FOR 

11601 

8 

- 6 . 7676 5E+02 

.8570 

6 . 52884E + 02 

1 .1320 

EL FOR 

1262 

4 

-2.22637E+02 

1.2840 

1 . 08921E+02 

.8560 

EL FOR 

1262 

5 

-4 . 89656E+01 

.7670 

2 . 43728E+01 

1.9350 

EL FOR 

1262 

8 

-1.02894E+02 

1.2850 

5 . 93157E+01 

.8550 

EL FOR 

16505 

2 

-9 . 22882E+01 

1.3640 

8 . 67772E+01 

.6080 

EL FOR 

16505 

3 

-4 . 114 28E+02 

. 8510 

3 . 11839E+02 

.5970 

EL FOR 

16505 

4 

- 9 . 4 7152E+01 

1 . 1040 

9 . 85466E+0 1 

.6490 
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STRUCTURE : 


ACCE 
ACCE 
ACCE 
ACCE 
ACCE 
ACCE 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 


cassini 

ID 


MIN/MAX SUMMARY 


PARAM2 : eigen 

VALUE 


PARAM3: 60 HZ 
OCCURRENCE 


RESPONSE: ALL 
VALUE OCCURRENCE 


10004 

3 

-5 . 88162E+02 

10004 

4 

-4 . 04593E+02 

10004 

5 

-1 . 75667E+03 

701 

3 

-3 . 12834E+02 

701 

4 

-2 . 75536E+02 

701 

5 

-1 . 14237E+03 

11601 

4 

-1 . 04183E+02 

11601 

5 

-3 . 40644E+01 

11601 

8 

-7 . 10632E+02 

1262 

4 

-2 . 36321E+02 

1262 

5 

-4 . 85338E+01 

1262 

8 

-1 . 04599E+02 

16505 

2 

-9 . 22948E+01 

16505 

3 

-4 . 08437E+02 

16505 

4 

-9 . 41314E+01 


1.3670 

4 . 89098E+02 

.7680 

1.9490 

5 . 36760E+02 

.6300 

.5140 

3 . 23813E + 02 

.6850 

1.2130 

3 . 78818E+02 

1.1050 

1.9610 

5 . 07715E+02 

.6280 

.7740 

-3 . 90032E+01 

.0000 

.5960 

5 . 70208E+01 

.8620 

1.1870 

4 . 88520E+01 

1.1300 

.8570 

6 . 11217E+02 

1.1320 

1.2840 

1 . 02161E+02 

.8560 

.7670 

2 . 44500E+01 

1.9350 

1.2850 

5 . 8518 0E + 0 1 

.8550 

1.3640 

8.67314E+01 

.6080 

.8510 

3 . 14572E+02 

.5970 

1.1040 

9 . 87876E+01 

.6490 
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MIN /MAX SUMMARY 


STRUCTURE : 

cassini 

ID 

P ARAM 2 

: ritz PARAM3 : f racl 

VALUE OCCURRENCE 

. 0 RESPONSE 

VALUE 

: ALL 

OCCURRENCE 

ACCE 

10004 

3 

-6 . 31701E+02 

1.3650 

5 . 02735E+02 

.7660 

ACCE 

10004 

4 

-3 . 91711E+02 

1 . 9500 

5 . 49198E+02 

.6310 

ACCE 

10004 

5 

-1 . 80007E+03 

.5140 

3 . 46303E+02 

.3990 

ACCE 

701 

3 

-3 . 14088E+02 

1.3620 

3 . 71505E+Q2 

1.1050 

ACCE 

701 

4 

-2 . 69071E+02 

1.9620 

5 . 09647E+02 

.6280 

ACCE 

701 

5 

- 1 . 146 87E + 03 

1 . 0300 

- 1 . 60672E+01 

.0000 

EL FOR 

11601 

4 

- 1 . 03659E+02 

.5970 

5 . 43767E+01 

.8610 

EL FOR 

11601 

5 

-3 . 48669E+01 

.6730 

4 . 51540E+01 

1.1320 

EL FOR 

11601 

8 

- 6 . 65697E+02 

.8570 

6 . 37747E+02 

1 . 1330 

EL FOR 

1262 

4 

- 2 . 214 76E+02 

1.2840 

1. 09841E+02 

.8560 

EL FOR 

1262 

5 

-4 . 853 02E + 01 

.7670 

2 . 4 9850E+01 

1 . 9370 

EL FOR 

1262 

8 

-1 . 0193 8E + 02 

1.2850 

5. 99338E+01 

.8550 

EL FOR 

16505 

2 

- 9 . 21733E + 01 

1.3640 

8 . 67974E+01 

. 6080 

EL FOR 

16505 

3 

-4 . 10596E+02 

.8510 

3 . 12612E+02 

.5970 

EL FOR 

16505 

4 

-9 . 37584E+01 

1.1040 

9 . 82706E+0 1 

.6490 
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STRUCTURE : 


ACCE 
ACCE 
ACCE 
ACCE 
ACCE 
ACCE 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 
EL FOR 


cassini 

ID 


MIN /MAX SUMMARY 


PARAM2: ritz 

VALUE 


PARAM3: frac.5 
OCCURRENCE 


RESPONSE: ALL 
VALUE OCCURRENCE 


10004 

3 

-6 . 21011E+02 

10004 

4 

-3 . 92906E+02 

10004 

5 

-1 . 62924E+03 

701 

3 

-3 . 62351E+02 

701 

4 

-2 . 683 99E+02 

701 

5 

-1.16865E+03 

11601 

4 

-1 . 12753E+02 

11601 

5 

-3 . 22312E+01 

11601 

8 

-6 . 03145E+02 

1262 

4 

-2 . 20645E+02 

1262 

5 

-4 . 57440E+01 

1262 

8 

-1 . 02720E+02 

16505 

2 

- 9 . 30055E+01 

16505 

3 

-4 . 08634E+02 

16505 

4 

-9 . 48909E+01 


1.2210 

5 . 85796E+02 

.7680 

1.1080 

5 . 29591E+02 

.6350 

.5080 

1 . 67419E+02 

.3470 

1.3630 

4 . 04263E+02 

1.1050 

1.0970 

5 . 1453 9E+02 

.6290 

1.0300 

-3.07396E+01 

.0000 

.5950 

5 . 31308E+01 

.8620 

1.1860 

4 . 20591E + 01 

1.1310 

.8530 

5 . 55450E+02 

1.1330 

1.2860 

1 . 08038E + 02 

.8550 

.7680 

2 . 43067E+01 

1.9390 

1.2870 

5 . 99723E+01 

.8550 

1.3620 

9 . 00367E+01 

.6070 

.8510 

3 . 10 949E+Q2 

.5980 

1.1030 

9 . 97088E+0 1 

.6450 
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MIN/MAX SUMMARY 

STRUCTURE: Cassini PARAM2 : ritz PARAM3 : frac.25 RESPONSE: ALL 

ID VALUE OCCURRENCE VALUE OCCURRENCE 


ACCE 

10004 

3 

ACCE 

10004 

4 

ACCE 

10004 

5 

ACCE 

701 

3 

ACCE 

701 

4 

ACCE 

701 

5 

EL FOR 

11601 

4 

EL FOR 

11601 

5 

EL FOR 

11601 

8 

EL FOR 

1262 

4 

EL FOR 

1262 

5 

EL FOR 

1262 

8 

EL FOR 

16505 

2 

EL FOR 

16505 

3 

EL FOR 

16505 

4 


-4 . 52083E+02 1.3610 
-3 . 88050E+02 1.9500 
-1 . 62037E+03 1.2890 
-3 . 23387E+02 1.3680 
-2 . 68010E+02 .5160 
- 1 . 15915E+03 .7770 
- 9 . 03325E+01 .5930 
-3 . 09406E+01 .6650 
-6 . 11457E+02 .8520 
-2 . 15326E+02 1.2890 
-4 . 86087E+01 .7680 
-9 . 91673E+01 1.2890 
-9 . 16432E+01 .8510 
-4 . 12496E+02 .8510 
-8 . 68751E+01 1.1050 


3.76871E+02 .9310 

5 . 85286E+02 .6280 

2 . 13292E+02 .6900 

3 . 62839E+02 1.1070 

4 . 66561E+02 .6310 

-4 . 01343E+01 .0000 

5 . 33065E+01 .8610 

4 . 00999E+01 1.1330 

5 . 70471E+02 1.1370 

1 . 11135E+02 .8530 

2 . 40371E+01 1.9370 

6 . 12986E+01 .8530 

8 . 55358E+01 .6020 

3 . 10529E+02 .5990 

1 . 05739E+02 .6460 
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